Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,085,498 G→C 5.4% intergenic (+164/‑431) PP_1864 → / → PP_1865 hypothetical protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,085,4980GC 5.4% 222.7 / 4.4 74intergenic (+164/‑431)PP_1864/PP_1865hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base G (37/33);  new base C (2/2);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

TAGCTTCTGCTCCAGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGT  >  NC_002947/2085350‑2085644
                                                                                                                                                    |                                                                                                                                                  
tAGCTTCTGCTCCAGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATgg                                                                                                                                                   >  2:1487779/1‑150 (MQ=255)
tAGCTTCTGCTCCAGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATgg                                                                                                                                                   >  2:1244132/1‑150 (MQ=255)
   cTTCTGCTCCAGCTTCTGCTTCTGCTTCTGCTTCTGTTTCTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGTGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATg                                                                                                                                                <  2:75608/150‑1 (MQ=255)
    ttCTGCTCCAGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTCTTGCTTTGGCTTTTGCTTCTAATCGCGCTGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGc                                                                                                                                               <  1:1851461/150‑1 (MQ=255)
      cTGCTCCAGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGccc                                                                                                                                             <  1:14965/150‑1 (MQ=255)
        gCTCCAGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGt                                                                                                                                           <  1:1487779/150‑1 (MQ=255)
          tCCAGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCa                                                                                                                                         >  1:2143712/1‑150 (MQ=255)
            cAGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAg                                                                                                                                       >  2:1821460/1‑150 (MQ=255)
            cAGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGCGCATGGATGCCCGTCAAg                                                                                                                                       <  2:1347088/150‑1 (MQ=255)
                    gcttctgcttctgcttctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGc                                                                                                                               <  2:2277755/150‑1 (MQ=255)
                       tctgcttctgcttctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGcccc                                                                                                                            >  1:2088297/1‑150 (MQ=255)
                        ctgcttctgcttctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCa                                                                                                                           <  2:2088297/150‑1 (MQ=255)
                            ttctgcttctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATCCATGCCCGTc                                                                                                                                          <  1:952905/131‑1 (MQ=255)
                            ttctgcttctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATCCATGCCCGTc                                                                                                                                          >  2:952905/1‑131 (MQ=255)
                             tctgcttctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCatcctg                                                                                                                      >  1:603691/1‑145 (MQ=255)
                             tctgcttctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATg                                                                                                                      <  2:167713/150‑1 (MQ=255)
                                gcttctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATgggg                                                                                                                   >  1:1708509/1‑150 (MQ=255)
                                      gcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCa                                                                                                             <  1:774825/150‑1 (MQ=37)
                                        ttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAgc                                                                                                           >  2:493543/1‑150 (MQ=35)
                                        ttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAgc                                                                                                           >  1:2343400/1‑150 (MQ=35)
                                        ttctgctCTTGCTTTGGCTTTTGCTTCTAACCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAgc                                                                                                           >  1:1655106/1‑150 (MQ=21)
                                              tctTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGtc                                                                                                     >  2:548140/1‑150 (MQ=32)
                                               ctTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCCGCAGCGCGGtcc                                                                                                    >  2:1689490/1‑150 (MQ=17)
                                                ttGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGtcct                                                                                                   >  1:445067/1‑150 (MQ=32)
                                                ttGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGtcct                                                                                                   >  1:1403426/1‑150 (MQ=32)
                                                  gCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGtcctcc                                                                                                 <  1:493543/150‑1 (MQ=32)
                                                    tttGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCg                                                                                               >  1:188650/1‑150 (MQ=32)
                                                       ggCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCCGGATGCGGCCTGCAGCGCGGacctccaagga                                                                                            <  2:1298403/150‑6 (MQ=255)
                                                       ggCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATgggg                                                                                                                   <  1:752951/127‑1 (MQ=33)
                                                       ggCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATgggg                                                                                                                   >  2:752951/1‑127 (MQ=33)
                                                       ggCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGa                                                                                            <  1:182977/150‑1 (MQ=32)
                                                        gCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATg                                                                                                                                                    <  2:376504/93‑1 (MQ=34)
                                                        gCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATg                                                                                                                                                    >  1:376504/1‑93 (MQ=34)
                                                         cTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGc                                                                                          >  2:2354192/1‑150 (MQ=32)
                                                          ttttGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCa                                                                                         >  1:1485376/1‑150 (MQ=32)
                                                          ttttGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCa                                                                                         >  1:1816162/1‑150 (MQ=32)
                                                                ttCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCg                                                                                   >  2:2122769/1‑150 (MQ=32)
                                                                ttCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCg                                                                                   >  2:1257480/1‑150 (MQ=32)
                                                                  cTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGa                                                                                 <  1:548140/150‑1 (MQ=32)
                                                                   tAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCaa                                                                                        <  1:822174/142‑1 (MQ=32)
                                                                   tAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCaa                                                                                        >  2:822174/1‑142 (MQ=32)
                                                                    aaGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGt                                                                               <  2:2343400/150‑1 (MQ=32)
                                                                    aaGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGt                                                                               >  1:244967/1‑150 (MQ=32)
                                                                    aaGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGt                                                                               >  2:647863/1‑150 (MQ=32)
                                                                      gcgcgcGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGa                                                                             <  2:2143712/150‑1 (MQ=32)
                                                                      gcgcgcGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGa                                                                             >  2:379780/1‑150 (MQ=32)
                                                                       cgcgcgGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAg                                                                            <  1:1821460/150‑1 (MQ=32)
                                                                        gcgcgGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAgg                                                                           <  1:275130/150‑1 (MQ=32)
                                                                          gcgGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAgggg                                                                         <  2:1816162/150‑1 (MQ=32)
                                                                            ggTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGAc                                                                       >  1:436567/1‑150 (MQ=32)
                                                                                 ttCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCgg                                                                  >  2:1978878/1‑150 (MQ=32)
                                                                                  tCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGa                                                                 <  2:1708509/150‑1 (MQ=32)
                                                                                      gCGATGCAGATTGCGACTTCAGGAGGCCGCGCGCAGGGCTTGCGGAGGGAGGTGACGGGGATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGagcgt                                                             >  1:318990/1‑150 (MQ=17)
                                                                                      gCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGagcgt                                                             >  1:1430631/1‑150 (MQ=32)
                                                                                        gATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGagcgtag                                                           <  2:1503980/150‑1 (MQ=32)
                                                                                          tGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGagcgtagcg                                                         >  1:474970/1‑150 (MQ=11)
                                                                                                          aggaggCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAgg                                         <  2:436567/150‑1 (MQ=31)
                                                                                                               gCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGTCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGc                                    >  1:1489330/1‑150 (MQ=17)
                                                                                                                  gAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGa                                 <  2:1430631/150‑1 (MQ=32)
                                                                                                                     cgcAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACgg                              <  1:553440/150‑1 (MQ=32)
                                                                                                                      gcAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGc                             >  1:2020602/1‑150 (MQ=33)
                                                                                                                          ggCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGTCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTg                         <  2:1489330/150‑1 (MQ=25)
                                                                                                                          ggCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTg                         <  1:1300989/150‑1 (MQ=35)
                                                                                                                           gCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCATAGCGTCGGGGCCGGATGAGGGGAGCCGACGGCTTTgg                        <  2:1485376/150‑1 (MQ=255)
                                                                                                                            cTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGGCTGCAGCGCGGTCCTCCCGGGCGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGt                       >  2:908191/1‑150 (MQ=25)
                                                                                                                             ttGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAggg                                        >  1:896134/1‑132 (MQ=32)
                                                                                                                             ttGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAggg                                        <  2:896134/132‑1 (MQ=32)
                                                                                                                               gCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTAc                    <  1:1107650/150‑1 (MQ=35)
                                                                                                                                   agggaggTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTg                <  2:1655106/150‑1 (MQ=35)
                                                                                                                                    gggaggTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGGTGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTgg               >  2:1027324/1‑150 (MQ=21)
                                                                                                                                    gggaggTGACGGGCATGGATGCCCGTCAAGCGCTGGGCACCAGGATGGGGTCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGATTTAGTTACTTTgg               >  2:2299820/1‑150 (MQ=255)
                                                                                                                                            aCGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGAcc       >  1:1942002/1‑150 (MQ=35)
                                                                                                                                            aCGGGCATCCATGCCCGTCACGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTg                <  2:2041718/141‑1 (MQ=255)
                                                                                                                                            aCGGGCATCCATGCCCGTCACGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTg                >  1:2041718/1‑141 (MQ=255)
                                                                                                                                             cGGGCATGGATGCCCGTCAAGCTCTGGGCCCCAGGATGGGGCCGAAAACGCCGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCa      <  2:318990/150‑1 (MQ=11)
                                                                                                                                               ggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTCGCCTTGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCaaa    <  1:1978878/150‑1 (MQ=255)
                                                                                                                                                gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAg   >  1:2084646/1‑150 (MQ=35)
                                                                                                                                                 cATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGt  >  1:256039/1‑150 (MQ=35)
                                                                                                                                                    |                                                                                                                                                  
TAGCTTCTGCTCCAGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGT  >  NC_002947/2085350‑2085644

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