Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,865,701 A→G 6.1% L76P (CTA→CCA)  purN ← phosphoribosylglycinamide formyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,865,7010AG 6.1% 1302.1 / 44.9 412L76P (CTA→CCA) purNphosphoribosylglycinamide formyltransferase
Reads supporting (aligned to +/- strand):  ref base A (195/191);  new base G (8/17);  total (203/209)
Fisher's exact test for biased strand distribution p-value = 9.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.13e-01

Only 100 of 418 total aligned reads displayed.

GCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTTGG  >  NC_002947/1865552‑1865850
                                                                                                                                                      |                                                                                                                                                     
    aTGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAATCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAgcgc                                                                                                                                                   <  2:7481815/150‑1 (MQ=255)
    aTGCAGGCCCTTGTACGTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAgcgc                                                                                                                                                   <  2:1533082/150‑1 (MQ=255)
     tGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGcgcg                                                                                                                                                  <  2:3232614/150‑1 (MQ=255)
        aGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCTAATGGcgcgcat                                                                                                                                               <  1:1907863/150‑4 (MQ=255)
          gCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATc                                                                                                                                             >  2:8118507/1‑150 (MQ=255)
          gCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATc                                                                                                                                             <  2:6491228/150‑1 (MQ=255)
            ccTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGa                                                                                                                                           >  2:607917/1‑150 (MQ=255)
             cTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGaa                                                                                                                                          >  1:921110/1‑150 (MQ=255)
             cTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGaa                                                                                                                                          >  2:3143785/1‑150 (MQ=255)
              ttGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCTTCGATGCCGCGCTAATGGCGCGCATCGacg                                                                                                                                         <  2:1889521/150‑3 (MQ=38)
              ttGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCTTCGATGCCGCGCTAATGGCGCGCATCGacg                                                                                                                                         <  2:5188959/150‑3 (MQ=38)
              ttGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAg                                                                                                                                         >  2:5838803/1‑150 (MQ=255)
               tGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAgg                                                                                                                                        <  2:6871468/150‑1 (MQ=255)
                gTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCTTCGATGCCGCGCTAATGGCGCGCATCGACGGc                                                                                                                                       <  2:569981/150‑1 (MQ=38)
                gTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCTTCGATGCCGCGCTAATGGCGCGCATCGACGGc                                                                                                                                       <  1:7622754/150‑1 (MQ=38)
                gTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGc                                                                                                                                       <  1:3146639/150‑1 (MQ=255)
                  aCTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCtt                                                                                                                                     <  1:4373643/150‑1 (MQ=255)
                    ttGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGTCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCa                                                                                                                                   >  1:1386697/1‑150 (MQ=255)
                    ttGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCTTCGATGCCGCGCTAATGGCGCGCATCGACGGCtttg                                                                                                                                   <  2:8522558/150‑3 (MQ=38)
                    ttGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCa                                                                                                                                   <  1:2944911/150‑1 (MQ=255)
                    ttGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCa                                                                                                                                   >  1:817517/1‑150 (MQ=255)
                    ttGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCa                                                                                                                                   >  1:3790283/1‑150 (MQ=255)
                    tacggCAGCACCGACCGGAGGAAGTTGAGCAGGAGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCAGCGAAGGCTTCa                                                                                                                                   >  2:5361516/4‑150 (MQ=255)
                     tGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCAc                                                                                                                                  <  2:7885548/150‑1 (MQ=255)
                      gggCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACg                                                                                                                                 <  1:8661799/150‑1 (MQ=255)
                      gggCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACg                                                                                                                                 <  1:3929329/150‑1 (MQ=255)
                            cagcGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTc                                                                                                                           >  1:2794325/1‑150 (MQ=255)
                             agcGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCg                                                                                                                          >  1:7002338/1‑150 (MQ=255)
                              gcGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGa                                                                                                                         >  1:6553143/1‑150 (MQ=255)
                               cGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGaa                                                                                                                        <  1:3606012/150‑1 (MQ=255)
                                   gggTGGCTGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                                    <  2:2906090/150‑1 (MQ=255)
                                   gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                                    <  2:6196913/150‑1 (MQ=255)
                                     gTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGgt                                                                                                                  >  1:840380/1‑150 (MQ=255)
                                     gTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGgt                                                                                                                  >  2:3614031/1‑150 (MQ=255)
                                      tGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGgtg                                                                                                                 <  1:4920848/150‑1 (MQ=255)
                                        gATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGtgtg                                                                                                               >  2:2861788/1‑150 (MQ=255)
                                        gATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGtgtg                                                                                                               >  2:7304573/1‑150 (MQ=255)
                                           gTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTc                                                                                                            >  1:5401788/1‑150 (MQ=255)
                                           gTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTc                                                                                                            >  1:2174615/1‑150 (MQ=255)
                                            ttGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCGAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCg                                                                                                           >  1:6972284/1‑150 (MQ=255)
                                              gAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAg                                                                                                         <  1:1797209/150‑1 (MQ=255)
                                               aGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGc                                                                                                        <  1:5215756/150‑1 (MQ=255)
                                               aGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGc                                                                                                        >  2:4792924/1‑150 (MQ=255)
                                                  aGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCacg                                                                                                     <  1:5157603/150‑1 (MQ=255)
                                                    gcggcCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCacgac                                                                                                   >  2:126900/1‑150 (MQ=255)
                                                         ccTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGt                                                                                              <  2:2260619/150‑1 (MQ=255)
                                                         ccTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGt                                                                                              <  1:4403840/150‑1 (MQ=255)
                                                          cTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTc                                                                                             >  1:5282423/1‑150 (MQ=255)
                                                             aTAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGAt                                                                                          <  1:2338508/150‑1 (MQ=255)
                                                               tgTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgc                                                                                        <  2:3514102/149‑1 (MQ=255)
                                                               aGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgc                                                                                        <  2:2720596/150‑1 (MQ=255)
                                                                gTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgcc                                                                                       <  2:2501505/150‑1 (MQ=255)
                                                                 tGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgccg                                                                                      <  2:4192971/150‑1 (MQ=255)
                                                                 tGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgccg                                                                                      >  2:5002955/1‑150 (MQ=255)
                                                                      cACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGc                                                                                 >  2:8646010/1‑150 (MQ=255)
                                                                         gAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTgc                                                                              >  2:2050065/1‑150 (MQ=255)
                                                                              cgccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGgcgc                                                                         <  2:7734614/150‑1 (MQ=255)
                                                                                   cTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGc                                                                    <  1:3263559/150‑1 (MQ=255)
                                                                                   cTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGc                                                                    >  1:8533811/1‑150 (MQ=255)
                                                                                    tGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCa                                                                   <  2:4392551/150‑1 (MQ=255)
                                                                                      aGGATGTGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAgg                                                                 >  1:6310275/1‑150 (MQ=255)
                                                                                           gcgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTa                                                            <  2:2418805/150‑1 (MQ=255)
                                                                                            cgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAg                                                           <  1:7859988/150‑1 (MQ=255)
                                                                                             gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCCTTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACGGGGTGTGGTTGACCACGACGCTGTCGATGCCGGCCGCGGCGGCGCGGTGCAGGCCGTAgg                                                          >  1:6631683/1‑150 (MQ=255)
                                                                                             gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAgg                                                          >  1:5638543/1‑150 (MQ=255)
                                                                                             gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTACAGGCCGTAgg                                                          >  1:733942/1‑150 (MQ=255)
                                                                                               aTGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCa                                                        >  1:6839108/1‑150 (MQ=255)
                                                                                                tGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAt                                                       >  1:2580009/1‑150 (MQ=255)
                                                                                                tGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAt                                                       >  1:1321600/1‑150 (MQ=255)
                                                                                                 gAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATc                                                      >  1:8143931/1‑150 (MQ=255)
                                                                                                  aaaCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCg                                                     <  2:3253335/150‑1 (MQ=255)
                                                                                                  aaaCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCg                                                     <  1:7088749/150‑1 (MQ=255)
                                                                                                  aaaCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCg                                                     <  1:3397151/150‑1 (MQ=255)
                                                                                                   aaCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCgg                                                    >  2:3759262/1‑150 (MQ=255)
                                                                                                      cGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGgcgc                                                 >  2:4496159/1‑150 (MQ=255)
                                                                                                         ccAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGt                                              <  2:8272243/150‑1 (MQ=255)
                                                                                                           aGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTg                                            <  2:340241/150‑1 (MQ=255)
                                                                                                           aGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTg                                            <  2:1091418/150‑1 (MQ=255)
                                                                                                             caccacCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGa                                          <  2:1321741/150‑1 (MQ=255)
                                                                                                              accaccAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGAGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGaa                                         <  1:5984657/150‑1 (MQ=255)
                                                                                                                caccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGTATCGGCGCGGTTGGAAAg                                       >  2:7180877/1‑150 (MQ=255)
                                                                                                                 accaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGc                                      >  2:6650667/1‑150 (MQ=255)
                                                                                                                   caGGTCCGGGGCAAAGCCGTCGATGCGCGCCTTGAGCCCGGCATCGAAGCCTTCACGGCCGACGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGACAGCAc                                    <  1:1499371/150‑1 (MQ=255)
                                                                                                                    aGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCAcc                                   <  2:6415391/150‑1 (MQ=255)
                                                                                                                     ggTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACcg                                  >  1:5543919/1‑150 (MQ=255)
                                                                                                                       tCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACcgcg                                >  2:8585642/1‑150 (MQ=255)
                                                                                                                         cGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACcgcgcg                              <  1:4507134/150‑1 (MQ=255)
                                                                                                                              cAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATAc                         <  1:6670798/150‑1 (MQ=255)
                                                                                                                               aaaGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACg                        <  2:4019305/150‑1 (MQ=255)
                                                                                                                                    cGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCg                   <  2:5113467/150‑1 (MQ=255)
                                                                                                                                     gTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCgg                  >  2:1725920/1‑150 (MQ=255)
                                                                                                                                      tCGATGCCGCGCTAATGGCGCGCATCGACGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCggg                 >  1:1013458/1‑150 (MQ=255)
                                                                                                                                      tCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTTGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCggg                 <  1:6166163/150‑1 (MQ=255)
                                                                                                                                      tCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCggg                 >  1:3539918/1‑150 (MQ=255)
                                                                                                                                      tCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCggg                 <  2:3957941/150‑1 (MQ=255)
                                                                                                                                      tCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACACTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCggg                 <  2:7846371/150‑1 (MQ=255)
                                                                                                                                               cgcCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTg         <  1:3992306/150‑1 (MQ=255)
                                                                                                                                                  cATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGt                                              <  1:5706329/110‑1 (MQ=255)
                                                                                                                                                      aGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTTgg  <  2:6530544/150‑1 (MQ=255)
                                                                                                                                                      aGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTTgg  <  2:5579901/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTTGG  >  NC_002947/1865552‑1865850

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