Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,606,281 +T 1.3% coding (13/1083 nt) mraY → phospho‑N‑acetylmuramoyl‑pentapeptide‑ transferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,606,2811.→T 1.3% 4.8 ‑3.8 ‑0.1234coding (13/1083 nt)mraYphospho‑N‑acetylmuramoyl‑pentapeptide‑ transferase
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (63/168):  total (66/168)
Fisher's exact test strand distribution p-value = 2.17e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.07e-01

Only 100 of 267 total aligned reads displayed.

CGAAAGGCTCGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTGGGGCCTTCTTCACCGCGCTGATCTTCGGCTTCCT  >  NC_014034/1606181‑1606381
                                                                                                     |                                                                                                     
cGAAAGGCTCGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGc                                                                                                        >  2:4654313/1‑101 (MQ=255)
cGAAAGGCTCGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGc                                                                                                        >  2:3934325/1‑101 (MQ=255)
 gAAAGGCTCGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.t                                                                                                      >  1:5284916/1‑101 (MQ=255)
 gAAAGGCTCGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.t                                                                                                      >  1:3229111/1‑101 (MQ=255)
  aaaGGCTCGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.tt                                                                                                     <  2:2390403/101‑1 (MQ=255)
  aaaGGCTCGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.tt                                                                                                     >  1:1936623/1‑101 (MQ=255)
   aaGGCTCGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.ttt                                                                                                    >  1:3906826/1‑100 (MQ=255)
   aaGGCTCGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.ttt                                                                                                    >  1:3638858/1‑100 (MQ=255)
   aaGGCTCGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTg                                                                                                    <  2:2886860/101‑1 (MQ=255)
    aGGCTCGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGc                                                                                                   <  2:5119935/101‑1 (MQ=255)
     tcttCGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCg                                                                                                  <  1:2815202/98‑1 (MQ=255)
     ggCTCGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.ttttt                                                                                                  >  2:2757711/1‑98 (MQ=255)
      gCTCGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGa                                                                                                 <  1:4654589/101‑1 (MQ=255)
        tCGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTggcaga                                                                                               >  2:2815041/1‑96 (MQ=255)
        tCGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAAc                                                                                               <  1:2509915/101‑1 (MQ=255)
        tCGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAAc                                                                                               <  2:925944/101‑1 (MQ=255)
         cGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAAcc                                                                                              <  2:4014647/101‑1 (MQ=255)
          gAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCTGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCt                                                                                             <  1:4898349/101‑1 (MQ=255)
            aGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCttt                                                                                           <  2:5784424/101‑1 (MQ=255)
             ggTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTc                                                                                          <  2:4387209/101‑1 (MQ=255)
             ggTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTc                                                                                          <  2:1565500/101‑1 (MQ=255)
               tttCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.gg                                                                                       <  1:4595635/101‑1 (MQ=255)
                ttCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGa                                                                                      <  1:4456672/101‑1 (MQ=255)
                 tCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGAc                                                                                     <  2:2287003/101‑1 (MQ=255)
                  cGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACg                                                                                    <  1:4044177/101‑1 (MQ=255)
                   gCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACgg                                                                                   <  2:4750748/101‑1 (MQ=255)
                   gCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACgg                                                                                   <  2:4076219/101‑1 (MQ=255)
                     tCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGGAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACgggg                                                                                 <  1:986826/101‑1 (MQ=255)
                     tCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACgggg                                                                                 <  1:1068750/101‑1 (MQ=255)
                        tagatGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCg                                                                              <  2:4882399/97‑1 (MQ=255)
                        tGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTTTGTACCCTTC.TGTCggggggg                                                                              >  1:1435159/1‑99 (MQ=255)
                        tGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCg                                                                              <  2:5497925/101‑1 (MQ=255)
                          gTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGAt                                                                            <  1:735299/101‑1 (MQ=255)
                          gTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGAt                                                                            <  1:318062/101‑1 (MQ=255)
                            tGACGCGGTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGTGAttt                                                                          <  1:89154/101‑1 (MQ=255)
                            tGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTTTTAAACTTTCTGGAGggggggggtt                                                                           >  2:573921/1‑96 (MQ=255)
                                gcgcTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTTTTGACCTTTT.TGACgggggggggtttttt                                                                      >  2:4490469/1‑91 (MQ=255)
                                gcgcTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATttcttc                                                                      <  1:1557118/101‑1 (MQ=255)
                                 cgcTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCa                                                                     <  2:3343441/101‑1 (MQ=255)
                                  gcTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCaa                                                                    <  2:3070612/101‑1 (MQ=255)
                                    tCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAAcc                                                                  <  1:1316375/101‑1 (MQ=255)
                                       gAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACctct                                                               <  1:3511075/101‑1 (MQ=255)
                                        aaaaCTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCtt                                                              <  1:2940636/101‑1 (MQ=255)
                                          aaCTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTcc                                                            <  1:2003183/101‑1 (MQ=255)
                                          aaCTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTcc                                                            <  2:479853/101‑1 (MQ=255)
                                          aaCTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTcc                                                            <  1:449307/101‑1 (MQ=255)
                                             tCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCt                                                         <  1:2249010/101‑1 (MQ=255)
                                             tCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCt                                                         <  1:241448/101‑1 (MQ=255)
                                              cGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTTTGAAACTTTC.GGGGGGGGGGGGGTTTTTTAACCTCGTCCTCtc                                                        >  2:1673126/1‑100 (MQ=255)
                                               ggTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTAc                                                       <  2:427201/101‑1 (MQ=255)
                                                 tCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACAt                                                     <  2:5295131/101‑1 (MQ=255)
                                                   aGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTTTGGACCTTTC.TGGCGGGGGCGGTGTCTTCTAACCCTTCCCCttcctcc                                                   >  2:3694994/1‑95 (MQ=255)
                                                   aGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTTTGAACCTTTC.GGGCGGGGGGGGGTTCTTCTACCTCTTCCCCCCCATCa                                                   >  1:1007735/1‑101 (MQ=255)
                                                     cccGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTTTGGACCTTTT.GGACGGGGGGGGGTTTTTTAAACTCTTCCCCTTCATCCcc                                                 >  1:4299407/1‑101 (MQ=255)
                                                      ccGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCt                                                <  2:3436718/101‑1 (MQ=255)
                                                       cGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCtt                                               <  2:3175743/101‑1 (MQ=255)
                                                       cGGCGCGGCAGAGGCCGTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCtt                                               <  2:639244/101‑1 (MQ=255)
                                                        ggCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTc                                              <  2:5655598/101‑1 (MQ=255)
                                                         gcgcGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTcc                                             <  1:1683782/101‑1 (MQ=255)
                                                           gcgGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCcgc                                           <  1:5648597/101‑1 (MQ=255)
                                                             ggCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTTTGGAACTTTT.TGGCGGGGGGGGGTTTTTTTACCTCTTTCCCTACATCACCTTcccccc                                         >  1:4991168/1‑97 (MQ=255)
                                                                 gagGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCGCCTTCCGCGCCGGt                                     <  1:4780211/101‑1 (MQ=255)
                                                                 gagGCCTACAAGGGAAAAGAATAATGCTTTACTGGCTTTGGGACCCCTTC.GG.CGGGGGCGGGTTCTTTAACCCCCTCCGCCCCCACACCCTCCGCGcccgc                                     >  1:621087/1‑98 (MQ=255)
                                                                  agGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTTTGGAACTTTT.GGACGGGGGGGGGTTTTTCTCTCTCTTCCGCCACACCACCCTCCCCGCCgggg                                    >  1:136743/1‑99 (MQ=255)
                                                                  agGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTTCGGACCTTTT.TGGCGGGGGGGGTTTTTTTTAACCCTTCCCCTTCATCCCCCTCCCCGcccggg                                    >  1:2981369/1‑97 (MQ=255)
                                                                  agGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTg                                    <  1:4725080/101‑1 (MQ=255)
                                                                  agGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTg                                    <  1:4602061/101‑1 (MQ=255)
                                                                  agGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTg                                    <  1:3107104/101‑1 (MQ=255)
                                                                   ggCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTTCGAACCCTTC.TGACGGGGGGGGGTTTTTCAACCTCTTCCCCCACATCACCCTccccccccgggg                                   >  1:1148709/1‑91 (MQ=255)
                                                                    gCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTTCGAACCTTTT.CGGCGGGGGGGGGTTTTTTAACCTCTTCCCCCACATCACCCTCCCCCCCGGGggg                                  >  1:2590049/1‑101 (MQ=255)
                                                                     ccTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTTCGGACCCTTC.GGGCGGGGGGGGTTTTTTTCAACTCTTCCCCCACATCACCCTCCTCGCCCGCgggg                                 >  2:1972014/1‑101 (MQ=255)
                                                                     ccTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTgggg                                 <  1:1494653/101‑1 (MQ=255)
                                                                       ttCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTGGGGcc                               <  2:4232168/101‑1 (MQ=255)
                                                                        tCAAGGGAAAAGAATAATGCTTTACTGGC.TTTTCTACCTTTT.TGTCGGCGGGGGTTTTTTTAATCTTTTTCGCTTCATCACCTTTTTCGCCCGTGGGGcgt                              >  2:1342943/1‑99 (MQ=255)
                                                                         cAAGGGAAAAGAATAATGCTTTACTGGC.TTTCGAACCTTTT.GGGCGGGGGGGGTTTTTTTTACCTCTTCCCCTACATCACCTTCCCCGCCCGGGGGGGctt                             >  1:5072150/1‑101 (MQ=255)
                                                                          aaGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTGGGGCcttc                            <  2:215493/101‑1 (MQ=255)
                                                                           aGGGAAAAGAATAATGCTTTACTGGC.TTTCGAACCTTTC.GGACGGGGGGGGGTTTTTCAACCTCTTCCGCTACATCACCTTCCCCGCCCGGGGGGCcttct                           >  1:1235530/1‑101 (MQ=255)
                                                                           aGGGAAAAGAATAATGCTTTACTGGC.TTTCGAAACCTTT.GGACGGGGGGGGGTTCTTTAACCTCTTCCGCTACATACCCTTCCGCCCCCCTGGGGGcctct                           >  2:266211/1‑101 (MQ=255)
                                                                           aGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTGGGGCcttct                           <  1:3678952/101‑1 (MQ=255)
                                                                             ggAAAAGAATAATGCTTTACTGGC.TTTCGAACCCTTC.GGACGGGGGGGGTTTTTTTAAACCCTTCCCCTACATCCCCCTCCCCCCCCGCGGGGGCttcttc                         >  2:5352288/1‑101 (MQ=255)
                                                                               aaaaGAATAATGCTTTACTGGC.TTTTGGACCTTTT.CGATGGGGGGGGTTTTTTCCACCCTTCCCCCTACCCCACCTTCCCCCCCCGTcggggcgttttttc                       >  1:5682496/1‑91 (MQ=255)
                                                                               aaaaGAATAATGCTTTACTGGC.TTTCGAACCTTTT.GGACGGGGGGGGTTTCTTTAACCTCTTCCGCCACATTACCTTCCGCGCCCGCGGGGCCctcttctc                       >  1:5609031/1‑99 (MQ=255)
                                                                                aaaGAATAATGCTTTACTGGC.TTTCGAACCTTTT.GGGCGGGGGGGATTTTTTTAACCTCTTCCGCTACATCACCTTCCGCGCCGGGGGGGGCTTCTTCAcc                      >  2:5152830/1‑101 (MQ=255)
                                                                                aaaGAATAATGCTTTACTGGC.TTTCGAACCTTTC.GGACGGGGGGGGGTTTTTCAACCTCTTCCCCCACATCACCTTCCCCCCCCGGGGGGCCCTCTTCTcc                      >  1:5219499/1‑101 (MQ=255)
                                                                                 aaGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTGGGGCCTTCTTCACcg                     <  2:141032/101‑1 (MQ=255)
                                                                                  aGAATAATGCTTTACTGGC.TTTCGAACCCTTC.GGGCGGGGGCGGTTTCTTCAACCTCTTTCCCCACATCACCCTCCCCGCCCGCGGGGCcttctttccaga                    >  2:5625191/1‑95 (MQ=255)
                                                                                    aataatGCTTTACTGGC.TTTCGAACCCTTC.GGACGGGGGCGATTTCTTCAAACTCCTCCGCTCCATCACCTTCCCCGCCGGTGGGGGCTTCCTCCcccccc                  >  2:859472/1‑97 (MQ=255)
                                                                                        atGCTTTACTGGC.TTGGGACCCTTTC.GGAAGGGGGGGGTTTTTTTTAACTCTTCCCCTTCAACACCCTCCGCCCCGGGGGGGCCTTGTTCTCCGCGCTGag              >  2:1619374/1‑100 (MQ=255)
                                                                                          gCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTGGGGCCTTCTTCACCGCGCTGATCt            <  2:2460548/101‑1 (MQ=255)
                                                                                           cTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTGGGGCCTTCTTCACCGCGCTGATCtt           <  2:3575187/101‑1 (MQ=255)
                                                                                             ttACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTGGGGCCTTCTTCACCGCGCTGATCTTCg         <  1:2163666/101‑1 (MQ=255)
                                                                                             ttACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTGGGGCCTTCTTCACCGCGCTGATCTTCg         <  2:5391020/101‑1 (MQ=255)
                                                                                                cTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTGGGGCCTTCTTCACCGCGCTGATCTTCGGCt      >  1:4749121/1‑101 (MQ=255)
                                                                                                 tGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTGGGGCCTTCTTCACCGCGCTGATCTTCGGCtt     >  1:513138/1‑101 (MQ=255)
                                                                                                  ggC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTGGGGCCTTCTTCACCGCGCTGATCTTCGGCTTc    >  2:2207183/1‑101 (MQ=255)
                                                                                                  ggC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTGGGGCCTTCTTCACCGCGCTGATCTTCGGCTTc    >  2:2151196/1‑101 (MQ=255)
                                                                                                  ggC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTGGGGCCTTCTTCACCGCGCTGATCTTCGGCTTc    >  1:5449092/1‑101 (MQ=255)
                                                                                                    c.TTGCTAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTGGGGCCTTCTTCACCGCGCTGATCTTCGGCTTCCt  <  1:4388908/101‑1 (MQ=255)
                                                                                                    c.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTGGGGCCTTCTTCACCGCGCTGATCTTCGGCTTCCt  <  1:3952518/101‑1 (MQ=255)
                                                                                                    c.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTGGGGCCTTCTTCACCGCGCTGATCTTCGGCTTCCt  <  2:1148628/101‑1 (MQ=255)
                                                                                                     |                                                                                                     
CGAAAGGCTCGAAGGTTTCGCTCGTGGTTGACGCGCTCAGAAAACTCGGTCAGCCCGGCGCGGCAGAGGCCTTCAAGGGAAAAGAATAATGCTTTACTGGC.TTGCGAACCTTTC.GGACGGGGGCGATTTCTTCAACCTCTTCCGCTACATCACCTTCCGCGCCGGTGGGGCCTTCTTCACCGCGCTGATCTTCGGCTTCCT  >  NC_014034/1606181‑1606381

Base quality scores:  ATCG  < 3 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG