Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,860,752 A→G 5.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,860,7520AG 5.1% 948.9 / 23.3 294L76P (CTA→CCA) purNphosphoribosylglycinamide formyltransferase
Reads supporting (aligned to +/- strand):  ref base A (158/120);  new base G (12/3);  total (170/124)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

Only 100 of 300 total aligned reads displayed.

GCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTTGG  >  NC_002947/1860603‑1860901
                                                                                                                                                      |                                                                                                                                                     
gCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTa                                                                                                                                                       <  2:5049301/150‑1 (MQ=255)
     tGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGcgcg                                                                                                                                                  <  2:2492332/150‑1 (MQ=255)
       cAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGc                                                                                                                                                >  2:6557439/1‑150 (MQ=255)
       cAGGCCCTCGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGc                                                                                                                                                >  1:3219141/1‑150 (MQ=255)
        aGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCa                                                                                                                                               >  2:8938813/1‑150 (MQ=255)
          gCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATc                                                                                                                                             >  1:7008259/1‑150 (MQ=255)
             cTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGaa                                                                                                                                          >  1:6332592/1‑150 (MQ=255)
             cTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGACGATGCGCGCCATTAGCGCGGCATCGaa                                                                                                                                          >  2:8137179/1‑150 (MQ=255)
                    ttGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCa                                                                                                                                   >  1:1988097/1‑150 (MQ=255)
                        gcagcaGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGc                                                                                                                               >  2:3400303/1‑150 (MQ=255)
                         cagcagCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACGGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGcc                                                                                                                              >  1:1806641/1‑150 (MQ=255)
                             agcGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCg                                                                                                                          <  2:7244849/150‑1 (MQ=255)
                                 aCGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACt                                                                                                                      >  1:382698/1‑150 (MQ=255)
                                   gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                                    <  2:4955951/150‑1 (MQ=255)
                                   gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                                    <  1:6299854/150‑1 (MQ=255)
                                   gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                                    <  1:5725312/150‑1 (MQ=255)
                                   gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                                    <  2:1271829/150‑1 (MQ=255)
                                   gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                                    <  1:2726754/150‑1 (MQ=255)
                                   gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                                    <  2:1357546/150‑1 (MQ=255)
                                     gTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGACCTGGgt                                                                                                                  >  2:7854849/1‑150 (MQ=255)
                                     gTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGgt                                                                                                                  <  1:2834771/150‑1 (MQ=255)
                                        gATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGtgtg                                                                                                               <  2:5646147/150‑1 (MQ=255)
                                         aTGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTgg                                                                                                              >  2:7515224/1‑150 (MQ=255)
                                         aTGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTgg                                                                                                              <  1:3499726/150‑1 (MQ=255)
                                           gTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTc                                                                                                            >  2:3193775/1‑150 (MQ=255)
                                           gTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTc                                                                                                            <  1:7786476/150‑1 (MQ=255)
                                           gTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTc                                                                                                            <  2:4079529/150‑1 (MQ=255)
                                           gTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTc                                                                                                            <  1:2186960/150‑1 (MQ=255)
                                            ttGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCg                                                                                                           >  2:7568618/1‑150 (MQ=255)
                                             tGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGa                                                                                                          <  2:5785683/150‑1 (MQ=255)
                                              gAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAg                                                                                                         >  2:3116785/1‑150 (MQ=255)
                                              gAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAg                                                                                                         >  1:3001993/1‑150 (MQ=255)
                                               aGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGc                                                                                                        <  1:5085279/150‑1 (MQ=255)
                                                  aGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCacg                                                                                                     >  1:5688483/1‑150 (MQ=255)
                                                    gcggcCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCacgac                                                                                                   >  2:2595415/1‑150 (MQ=255)
                                                    gcggcCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCacgac                                                                                                   >  1:1539974/1‑150 (MQ=255)
                                                         ccTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGt                                                                                              >  2:1966200/1‑150 (MQ=255)
                                                         ccTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGt                                                                                              <  2:9074098/150‑1 (MQ=255)
                                                          cTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTc                                                                                             >  1:3953659/1‑150 (MQ=255)
                                                           tGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACAACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCg                                                                                            >  1:488395/1‑150 (MQ=255)
                                                               aGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTTGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgc                                                                                        >  2:8363002/1‑150 (MQ=255)
                                                                   gcgcACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgccggc                                                                                    >  1:8210449/1‑150 (MQ=255)
                                                                       aCGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCt                                                                                >  1:4013894/1‑150 (MQ=255)
                                                                        ttcgtgcgccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTg                                                                               >  1:8309067/7‑150 (MQ=255)
                                                                         gAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTgc                                                                              >  2:3003433/1‑150 (MQ=255)
                                                                         gAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTgc                                                                              <  2:2836362/150‑1 (MQ=255)
                                                                            gccgccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTgcggc                                                                           >  2:601528/1‑150 (MQ=255)
                                                                             ccgccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTgcggcg                                                                          >  2:5360333/1‑150 (MQ=255)
                                                                               gccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGcgcg                                                                        >  1:4530225/1‑150 (MQ=255)
                                                                                ccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGt                                                                       <  2:7766631/150‑1 (MQ=255)
                                                                                 cACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGtt                                                                      >  2:825238/1‑150 (MQ=255)
                                                                                 cACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGtt                                                                      >  1:4308539/1‑150 (MQ=255)
                                                                                      aGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAgg                                                                 <  1:1054468/150‑1 (MQ=255)
                                                                                        gATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGcc                                                               <  2:1654524/150‑1 (MQ=255)
                                                                                            cgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAg                                                           <  2:3895114/150‑1 (MQ=255)
                                                                                             gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAgg                                                          >  2:6369335/1‑150 (MQ=255)
                                                                                             gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAgg                                                          >  1:7509398/1‑150 (MQ=255)
                                                                                              cATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGc                                                         <  2:6459336/150‑1 (MQ=255)
                                                                                               aTGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCa                                                        >  1:6911918/1‑150 (MQ=255)
                                                                                                tGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAt                                                       <  2:1028321/150‑1 (MQ=255)
                                                                                                 gAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATc                                                      >  1:2512190/1‑150 (MQ=255)
                                                                                                   aaCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCgg                                                    >  1:872805/1‑150 (MQ=255)
                                                                                                   aaCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCgg                                                    >  2:5192262/1‑150 (MQ=255)
                                                                                                    aCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGgc                                                   >  2:6275644/1‑150 (MQ=255)
                                                                                                      cGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGgcgc                                                 >  1:8665687/1‑150 (MQ=255)
                                                                                                         ccAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGt                                              <  2:6098436/150‑1 (MQ=255)
                                                                                                          cAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGtt                                             <  1:6568708/150‑1 (MQ=255)
                                                                                                          cAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGtt                                             <  1:7096173/150‑1 (MQ=255)
                                                                                                           aGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTg                                            <  2:1229842/150‑1 (MQ=255)
                                                                                                              accaccAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGaa                                         <  2:4346637/150‑1 (MQ=255)
                                                                                                                caccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAATGCTTCACGGCCGTCGAACTGGGTGTGGTCGATCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGaaac                                       >  1:5404389/1‑149 (MQ=255)
                                                                                                                caccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAgg                                                                 >  1:7901053/1‑124 (MQ=255)
                                                                                                                caccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAgg                                                                 <  2:7901053/124‑1 (MQ=255)
                                                                                                                caccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAg                                       >  1:4624665/1‑150 (MQ=255)
                                                                                                                caccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAg                                       <  2:7637685/150‑1 (MQ=255)
                                                                                                                caccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCAGTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCGGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAg                                       >  2:4526673/1‑150 (MQ=255)
                                                                                                                 accaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGTGTGTGGTCGAGCATGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGc                                      <  2:4435809/150‑1 (MQ=255)
                                                                                                                  ccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCa                                     <  1:8273963/150‑1 (MQ=255)
                                                                                                                  ccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCa                                     >  2:815739/1‑150 (MQ=255)
                                                                                                                    aGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCAcc                                   >  2:882696/1‑150 (MQ=255)
                                                                                                                    aGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCAcc                                   <  1:5501349/150‑1 (MQ=255)
                                                                                                                     ggTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACcg                                  >  2:6978233/1‑150 (MQ=255)
                                                                                                                      gTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACcgc                                 <  2:3332908/150‑1 (MQ=255)
                                                                                                                       tCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACcgcg                                >  2:1511952/1‑150 (MQ=255)
                                                                                                                         cGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACcgcgcg                              >  1:5872327/1‑150 (MQ=255)
                                                                                                                          ggggCAAAGCCGTCGATGCGCGCCAGTAGTGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCgg                             >  2:2035929/1‑150 (MQ=255)
                                                                                                                           ccgtgaaGCCTTCGATGCCGCGCTAATGGCGCGCATCGACGGCTTTACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGa                            >  1:1592794/6‑150 (MQ=37)
                                                                                                                           ccgtgaaGCCTTCGATGCCGCGCTAATGGCGCGCATCGACGGCTTTACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGa                            >  1:6186443/6‑150 (MQ=37)
                                                                                                                             gCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATa                          <  2:4303292/150‑1 (MQ=255)
                                                                                                                             gCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATa                          <  2:3441700/150‑1 (MQ=255)
                                                                                                                              cAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATAc                         <  1:1771193/150‑1 (MQ=255)
                                                                                                                                     ttCGATGCCGCGCTAATGGCGCGGCATCGAAGGCTTCACGGCCGTCGACCTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCGGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCg                   >  2:6912001/2‑150 (MQ=255)
                                                                                                                                        gATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCt               >  2:3634220/1‑150 (MQ=255)
                                                                                                                                          tgccgcgcTAATGGCGCGCATCGACGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGt             >  1:6833732/4‑150 (MQ=255)
                                                                                                                                           gcgcgcCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTc            <  1:2117615/150‑1 (MQ=255)
                                                                                                                                           gcgcgcCATGAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTc            <  1:413594/150‑1 (MQ=255)
                                                                                                                                              gcgcCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCt          >  2:4661115/1‑150 (MQ=255)
                                                                                                                                                 ccATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGcc       >  1:8948948/1‑150 (MQ=255)
                                                                                                                                                      aGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTTgg  <  2:7414808/150‑1 (MQ=255)
                                                                                                                                                      aGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTTgg  <  2:6751656/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTTGG  >  NC_002947/1860603‑1860901

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: