Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,350,508 T→A 8.4% K11I (AAA→ATA) ‡ PP_1176 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,350,5080TA 8.4% 259.7 / 12.6 85K11I (AAA→ATA) ‡PP_1176hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (41/36);  new base A (3/4);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

AGCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTAGTCCCC  >  NC_002947/1350363‑1350657
                                                                                                                                                 |                                                                                                                                                     
tcccGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGc                                                                                                                                                   >  2:183357/3‑150 (MQ=255)
    gTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGcc                                                                                                                                               >  1:739278/1‑150 (MQ=255)
         tGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAAATGGCAGCCttttt                                                                                                                                          <  2:1360217/150‑1 (MQ=255)
          gACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCttgttg                                                                                                                                         >  1:1604953/1‑150 (MQ=255)
           aCGAGAATGGGATTTCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCTATTTGGCAGCTTTGTTGa                                                                                                                                        >  2:1295752/1‑150 (MQ=255)
            cGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGAt                                                                                                                                       <  2:1604953/150‑1 (MQ=255)
            cGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGAt                                                                                                                                       >  2:554951/1‑150 (MQ=255)
            cGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGAt                                                                                                                                       <  2:127186/150‑1 (MQ=255)
             gagaACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATg                                                                                                                                      <  2:1409803/150‑1 (MQ=255)
                 aCGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCa                                                                                                                                  <  2:1003913/150‑1 (MQ=255)
                        tCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGt                                                                                                                           <  2:2281512/150‑1 (MQ=255)
                          cGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGtca                                                                                                                         >  2:46776/1‑150 (MQ=255)
                            cTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGTAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATCTGGCAGCCTTGTTGATGCTCAGCCCGGtcaat                                                                                                                       <  1:1295752/150‑1 (MQ=255)
                            cTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGtcaat                                                                                                                       <  1:832258/150‑1 (MQ=255)
                                gtgAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTcaatcaat                                                                                                                   <  2:1547121/150‑1 (MQ=255)
                                  gAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATcc                                                                                                                 <  1:774947/150‑1 (MQ=255)
                                  gAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGAACAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATcc                                                                                                                 <  1:2344867/150‑1 (MQ=255)
                                   aGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATccc                                                                                                                >  2:912487/1‑150 (MQ=255)
                                            aaGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACg                                                                                                       <  2:2001575/150‑1 (MQ=255)
                                                gcCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGcccc                                                                                                   >  1:1311921/1‑150 (MQ=255)
                                                      aaaCGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAAATGGCAGCCt                                                                                                                                              >  2:192078/1‑101 (MQ=255)
                                                      aaaCGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAAATGGCAGCCt                                                                                                                                              <  1:192078/101‑1 (MQ=255)
                                                        aCGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATTCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCt                                                                                           <  1:554951/150‑1 (MQ=255)
                                                         cGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTg                                                                                          >  1:1898675/1‑150 (MQ=255)
                                                           ttATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGCAACGTCCTGTTCTTCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGAt                                                                                        >  2:1594866/1‑150 (MQ=255)
                                                            tATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATa                                                                                       >  1:194721/1‑150 (MQ=255)
                                                            tATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATa                                                                                       <  1:183357/150‑1 (MQ=255)
                                                               cgcgATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGtcaat                                                                                                                       <  1:1054697/115‑1 (MQ=255)
                                                               cgcgATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGtcaat                                                                                                                       >  2:1054697/1‑115 (MQ=255)
                                                                  gATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGAt                                                                                 >  2:13572/1‑150 (MQ=255)
                                                                     gTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCttgt                                                                                                                                           >  2:954400/1‑89 (MQ=255)
                                                                     gTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCttgt                                                                                                                                           <  1:954400/89‑1 (MQ=255)
                                                                        cAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACAcg                                                                           >  2:1688524/1‑150 (MQ=255)
                                                                        cAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACAcg                                                                           >  1:569171/1‑150 (MQ=255)
                                                                        cAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACAcg                                                                           >  2:1014362/1‑150 (MQ=255)
                                                                         aaaaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCGGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACAcgc                                                                          >  2:1662682/1‑150 (MQ=255)
                                                                           aaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAAATGGCAGCCt                                                                                                                                              <  1:2259574/80‑1 (MQ=255)
                                                                           aaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAAATGGCAGCCt                                                                                                                                              >  2:2259574/1‑80 (MQ=255)
                                                                                        cGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGttc                                                           <  1:1662682/150‑1 (MQ=255)
                                                                                          gCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCTACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGttctt                                                         <  1:46776/150‑1 (MQ=255)
                                                                                          gCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGttctt                                                         >  2:1448815/1‑150 (MQ=255)
                                                                                             gCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTtctc                                                      <  1:1448815/150‑1 (MQ=255)
                                                                                                gTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAaagaag                                                              <  2:2413878/139‑1 (MQ=255)
                                                                                                gTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAaagaag                                                              >  1:2413878/1‑139 (MQ=255)
                                                                                                 tAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCaaa                                                                  <  2:1655611/134‑1 (MQ=255)
                                                                                                 tAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCaaa                                                                  >  1:1655611/1‑134 (MQ=255)
                                                                                                   ggACACGTCCTGTTCTGCTTCTGCCACGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCcttctt                                                <  1:342407/150‑1 (MQ=255)
                                                                                                     aaaCGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATcaca                                                                             >  2:1963474/1‑119 (MQ=255)
                                                                                                     aaaCGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATcaca                                                                             <  1:1963474/119‑1 (MQ=255)
                                                                                                      aaCGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCaa                                             <  1:912487/150‑1 (MQ=255)
                                                                                                          tCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTcc                                         >  2:1734898/1‑150 (MQ=255)
                                                                                                           ccTGTTCTGCTTCTGCCAAGAGATAAAATTGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCa                                        <  2:553366/150‑1 (MQ=255)
                                                                                                             tGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATcaca                                                                             <  2:1529909/111‑1 (MQ=255)
                                                                                                             tGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATcaca                                                                             >  1:1529909/1‑111 (MQ=255)
                                                                                                             tGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATg                                      >  2:1684939/1‑150 (MQ=255)
                                                                                                                tctgcttctgcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGa                                                                                  <  2:518343/103‑1 (MQ=255)
                                                                                                                tctgcttctgcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGa                                                                                  >  1:518343/1‑103 (MQ=255)
                                                                                                                tctgcttctgcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACCGTGGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGAtc                                   >  1:1205873/1‑150 (MQ=255)
                                                                                                                  tgcttctgcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGAtctc                                 >  1:1034041/1‑150 (MQ=255)
                                                                                                                   gcttctgcCAGGCAATCAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCg                                <  2:1723325/150‑1 (MQ=255)
                                                                                                                    attcaaccAGGCAAGAAACCGGCTGCAATTTGGTAGCATTGTTGATGATCAGCTCGGTCAATCAATCTCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGc                               >  2:1027042/7‑150 (MQ=255)
                                                                                                                     ttctgcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCt                              >  2:2008917/1‑150 (MQ=255)
                                                                                                                     ttctgcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCt                              >  1:277790/1‑150 (MQ=255)
                                                                                                                       ctgcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAACCAATCCCAGCTCAACGCCCCACCCGTc                                                                                            >  2:1771373/1‑86 (MQ=255)
                                                                                                                       ctgcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAACCAATCCCAGCTCAACGCCCCACCCGTc                                                                                            <  1:1771373/86‑1 (MQ=255)
                                                                                                                         gcCNGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATc                          >  2:1505551/1‑150 (MQ=255)
                                                                                                                           cAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCa                        >  2:1299980/1‑150 (MQ=255)
                                                                                                                           cAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCa                        >  1:1576139/1‑150 (MQ=255)
                                                                                                                             ggCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTc                                                                                    >  2:150081/1‑88 (MQ=255)
                                                                                                                             ggCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTc                                                                                    <  1:150081/88‑1 (MQ=255)
                                                                                                                              gCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATgg                     <  1:1684939/150‑1 (MQ=255)
                                                                                                                              gCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATgg                     <  2:739278/150‑1 (MQ=255)
                                                                                                                                aaTAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGa                   >  2:616706/1‑150 (MQ=255)
                                                                                                                                  tAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAAc                 <  1:1734898/150‑1 (MQ=255)
                                                                                                                                     aaaGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACggg              <  2:194721/150‑1 (MQ=255)
                                                                                                                                      aaGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGt             <  1:1014362/150‑1 (MQ=255)
                                                                                                                                      aaGGCTGCCAAATGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGt             <  1:939366/150‑1 (MQ=255)
                                                                                                                                      aaGGCTGCCAAATGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGt             >  2:939366/1‑150 (MQ=255)
                                                                                                                                        ggCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTg           <  1:1594866/150‑1 (MQ=255)
                                                                                                                                          cTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGt         >  2:658290/1‑150 (MQ=255)
                                                                                                                                           tGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTa        >  1:1688369/1‑150 (MQ=255)
                                                                                                                                              cATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGttc                                                           <  2:1403900/96‑1 (MQ=255)
                                                                                                                                              cATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGttc                                                           >  1:1403900/1‑96 (MQ=255)
                                                                                                                                               aTTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTtctc                                                      <  2:105494/100‑1 (MQ=255)
                                                                                                                                               aTTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTtctc                                                      >  1:105494/1‑100 (MQ=255)
                                                                                                                                                 ttGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTAGTcccc  <  1:1505551/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
AGCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTAGTCCCC  >  NC_002947/1350363‑1350657

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