Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,076,841 G→C 5.1% R72P (CGC→CCC)  mreC → cell shape‑determining protein MreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,076,8410GC 5.1% 316.7 / 4.8 97R72P (CGC→CCC) mreCcell shape‑determining protein MreC
Reads supporting (aligned to +/- strand):  ref base G (46/46);  new base C (2/3);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 104 total aligned reads displayed.

GATGGGCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCATTCTGATCAACAAGG  >  NC_002947/1076695‑1076990
                                                                                                                                                  |                                                                                                                                                     
gATGGGCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCt                                                                                                                                                    <  1:1155282/150‑1 (MQ=255)
   gggCCGGGTGCTGACGGAGTAGTACTCGCTTACCGACTGGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCa                                                                                                                                                 <  2:21886/150‑1 (MQ=255)
     gCCTTGTGCTCATGGCGTCGCACTGCACCCCCGATTCGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCCGTTCGGTAGCCCAACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGa                                                                                                                                               <  2:1109762/150‑1 (MQ=255)
      ccTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGaa                                                                                                                                              <  1:1162040/150‑1 (MQ=255)
      ccTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGaa                                                                                                                                              <  2:1074611/150‑1 (MQ=255)
      ccTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGCGCCCCTGCAgca                                                                                                                                              <  1:1315475/150‑3 (MQ=255)
          gTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGGAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTg                                                                                                                                          >  2:401199/1‑150 (MQ=255)
          gTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCAGCTTCTGCAGGCGCCCCTGCagcagcag                                                                                                                                          <  1:962554/150‑3 (MQ=255)
             cTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTggcg                                                                                                                                       <  1:2397936/150‑1 (MQ=255)
                aTGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGcc                                                                                                                                    >  1:2296510/1‑150 (MQ=255)
                aTGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGcc                                                                                                                                    >  2:411361/1‑150 (MQ=255)
                  ggAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCt                                                                                                                                  <  2:1447392/150‑1 (MQ=255)
                  ggAGTCGTACTGGAGCACCGATTTGCCAGAGCGTTCCTGGCAAGGTGTTGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCt                                                                                                                                  <  2:2108981/150‑1 (MQ=255)
                      tCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGAcc                                                                                                                              >  2:1979294/1‑150 (MQ=255)
                         tACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTCACCGAg                                                                                                                           <  2:50554/150‑1 (MQ=255)
                         tACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCAGCTTCTGCAGGCGCCCCTGCAGCAGCAgg                                                                                                                                         <  1:997389/136‑1 (MQ=255)
                         tACTGGATCACCGATTCGCCACAGCGCGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCCCCGAGCTGATTGACGAAAACGACAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGACCCCCCTCCGATGCAGGCCGCCCCGACCgtg                                                                                                                           >  2:185916/1‑148 (MQ=255)
                            tGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAg                                                                                                                        >  1:1446444/1‑150 (MQ=255)
                              gATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGaa                                                                                                                      >  1:1387487/1‑150 (MQ=255)
                                 cACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCAGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGt                                                                                                                   >  1:1376660/1‑150 (MQ=255)
                                  aCCGATTTGCCACAGGGTGCCTGGCAAGCTGTCGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAACACCGAAGCCCTCCTACTGCAGCGGCGCCTGCAGAAGCTGGCGGGCCTGACCGACGAGGACGTg                                                                                                                  >  2:1210834/1‑150 (MQ=255)
                                  aCCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTg                                                                                                                  >  1:894296/1‑150 (MQ=255)
                                  aCCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCCGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTg                                                                                                                  >  2:1713413/1‑150 (MQ=255)
                                   ccGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGc                                                                                                                 <  1:1879346/150‑1 (MQ=255)
                                     gATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGGGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCt                                                                                                                                                    >  1:184438/1‑113 (MQ=255)
                                     gATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGGGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCt                                                                                                                                                    <  2:184438/113‑1 (MQ=255)
                                            cacaGCGTGCCTGGCAAGGTGTGGCCGGCCGGTTCGGTAGGCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGAAAGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGcgcg                                                                                                        >  1:1855912/1‑150 (MQ=255)
                                            cacaGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGcgcg                                                                                                        >  1:516424/1‑150 (MQ=255)
                                            cacaGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGcgcg                                                                                                        >  2:747610/1‑150 (MQ=255)
                                             acaGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGaa                                                                                                                                              <  1:365693/111‑1 (MQ=255)
                                             acaGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGaa                                                                                                                                              >  2:365693/1‑111 (MQ=255)
                                             acaGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGa                                                                                                                                >  1:2048335/1‑125 (MQ=255)
                                             acaGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGa                                                                                                                                <  2:2048335/125‑1 (MQ=255)
                                             acaGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGa                                                                                                       >  2:2326044/1‑150 (MQ=255)
                                              caGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCCGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAg                                                                                                      <  1:1713413/150‑1 (MQ=255)
                                                gCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGtt                                                                                                    <  1:2105183/150‑1 (MQ=255)
                                                gCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCAGCTGCCCGAGtt                                                                                                    <  2:532339/150‑1 (MQ=255)
                                                  gTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGc                                                                                                  <  1:1590490/150‑1 (MQ=255)
                                                   tGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCt                                                                                                 <  1:747610/150‑1 (MQ=255)
                                                      cTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCaa                                                                                              <  1:411361/150‑1 (MQ=255)
                                                           aaGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCtt                                                                                         <  2:1392112/150‑1 (MQ=255)
                                                             ggTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCg                                                                                       >  2:747724/1‑150 (MQ=255)
                                                               tgtgGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGgc                                                                                     <  1:1166809/150‑1 (MQ=255)
                                                                     cggccAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTggg                                                                               >  2:1415450/1‑149 (MQ=255)
                                                                            ttCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGAcc                                                                                                                              <  2:1929292/96‑1 (MQ=255)
                                                                            ttCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGAcc                                                                                                                              >  1:1929292/1‑96 (MQ=255)
                                                                              cGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGaaaa                                                                      <  2:1446444/150‑1 (MQ=255)
                                                                                 tAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGa                                                                                                       <  1:316733/114‑1 (MQ=255)
                                                                                 tAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGa                                                                                                       >  2:316733/1‑114 (MQ=255)
                                                                                 tAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAAACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGt                                                                   <  1:135926/150‑1 (MQ=255)
                                                                                      gCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGt                                                                   <  2:70581/145‑1 (MQ=255)
                                                                                      gCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGt                                                                   >  1:70581/1‑145 (MQ=255)
                                                                                      gCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGATAACGTCTTTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAg                                                              >  2:694673/1‑150 (MQ=255)
                                                                                        aCCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTg                                                            >  1:1826505/1‑150 (MQ=255)
                                                                                          cgggCTGATTGCCGCAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACTGCGCAGAACGCGAGGCTGCGCGCGTTGCTCCACGCTTCGCCGCTGGTCAACGAAAAGGTTGACGTGGc                                                          <  2:89692/147‑1 (MQ=255)
                                                                                             gCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGa                                                       <  2:516424/150‑1 (MQ=255)
                                                                                              cTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAg                                                      >  2:1292758/1‑150 (MQ=255)
                                                                                                gATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAgg                                                                    <  1:286448/134‑1 (MQ=255)
                                                                                                gATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAgg                                                                    >  2:286448/1‑134 (MQ=255)
                                                                                                  ttGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAAGTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGa                                                  <  2:645169/150‑1 (MQ=255)
                                                                                                   tGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGAt                                                 <  2:894296/150‑1 (MQ=255)
                                                                                                        aaaaCGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGgtg                                            <  1:185916/150‑1 (MQ=255)
                                                                                                            cGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCaa                                                                            >  2:2284243/1‑114 (MQ=255)
                                                                                                            cGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCaa                                                                            <  1:2284243/114‑1 (MQ=255)
                                                                                                             gagaAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGc                                                     <  1:2266613/136‑1 (MQ=255)
                                                                                                             gagaAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGc                                                     >  2:2266613/1‑136 (MQ=255)
                                                                                                               gaAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGGCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGAccc                                     <  1:2326044/150‑1 (MQ=255)
                                                                                                                     cAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAAccc                               >  2:95460/1‑150 (MQ=255)
                                                                                                                     cAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAAccc                               >  2:117554/1‑150 (MQ=255)
                                                                                                                      aaGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTg                                                                                 <  1:1121782/99‑1 (MQ=255)
                                                                                                                      aaGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTg                                                                                 >  2:1121782/1‑99 (MQ=255)
                                                                                                                      aaGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCg                              >  2:330019/1‑150 (MQ=255)
                                                                                                                       aGACCGAAGCCCTGATGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGt                             >  2:215960/1‑150 (MQ=255)
                                                                                                                        gACCGATGCCTGGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGtt                            <  2:1068428/150‑1 (MQ=255)
                                                                                                                        gACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGtt                            >  2:2184818/1‑150 (MQ=255)
                                                                                                                        gACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGtt                            >  1:833019/1‑150 (MQ=255)
                                                                                                                        gACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGAGGCTGCGCGAGTTGCTCAACTCtt                                                                                         <  1:1874519/89‑1 (MQ=255)
                                                                                                                        gACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGAGGCTGCGCGAGTTGCTCAACTCtt                                                                                         >  2:1874519/1‑89 (MQ=255)
                                                                                                                            gAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTAcc                        >  2:456106/1‑150 (MQ=255)
                                                                                                                              gccccTGCAGCAGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCg                              >  1:1293894/3‑142 (MQ=255)
                                                                                                                              gccccTGCAGCAGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCg                              <  2:1293894/140‑1 (MQ=255)
                                                                                                                              gccccTGCAGCAGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCa                      >  2:1595840/3‑150 (MQ=255)
                                                                                                                              gccccTGCAGCAGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCa                      >  2:320585/3‑150 (MQ=255)
                                                                                                                              aGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCa                      >  1:915024/1‑150 (MQ=255)
                                                                                                                                  ctgctgctgCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGt                             <  2:2282844/139‑1 (MQ=255)
                                                                                                                                  ctgctgctgCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGt                             >  1:2282844/1‑139 (MQ=255)
                                                                                                                                    gctgctgcAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCAt                >  1:627046/1‑150 (MQ=255)
                                                                                                                                        ctgcAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCACCCACTCGACGCTGGTCAACGCATGTGTTCTAGTGGTCTAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCATTCTg            <  1:296494/150‑1 (MQ=255)
                                                                                                                                         tgcAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAAc                                                                           <  2:1457320/86‑1 (MQ=255)
                                                                                                                                         tgcAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAAc                                                                           >  1:1457320/1‑86 (MQ=255)
                                                                                                                                          gcAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCg                                               >  2:1380499/1‑113 (MQ=255)
                                                                                                                                          gcAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCg                                               <  1:1380499/113‑1 (MQ=255)
                                                                                                                                             ggggCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCTCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCATTCTGATcaa       <  2:1581625/150‑1 (MQ=255)
                                                                                                                                                  gcCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCAcc                    <  1:341385/132‑1 (MQ=255)
                                                                                                                                                  gcCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCAcc                    >  2:341385/1‑132 (MQ=255)
                                                                                                                                                  gcCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCATTCTGATCAACAAgg  <  2:833019/150‑1 (MQ=255)
                                                                                                                                                  gcCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCATTCTGATCAACAAgg  >  2:2001134/1‑150 (MQ=255)
                                                                                                                                                  gcCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCCTATCGGTGTCGACCCCAACCCGTTTACCCACCGCATTCTGATCAACAAgg  <  1:1979294/150‑1 (MQ=255)
                                                                                                                                                  gcCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGAGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCATTCTGATCAACAAgg  <  2:1657096/150‑1 (MQ=255)
                                                                                                                                                  gcCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGAGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCATTCTGATCAACAAgg  >  1:1657096/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GATGGGCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCATTCTGATCAACAAGG  >  NC_002947/1076695‑1076990

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: