Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,043,722 T→A 7.2% intergenic (+18/+16) gmk → / ← PP_5297 deoxyguanylate kinase/guanylate kinase/amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,043,7220TA 7.2% 243.7 / 8.6 83intergenic (+18/+16)gmk/PP_5297deoxyguanylate kinase/guanylate kinase/amino acid transporter
Reads supporting (aligned to +/- strand):  ref base T (39/38);  new base A (2/4);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

ACTACGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAACACCAG  >  NC_002947/6043574‑6043869
                                                                                                                                                    |                                                                                                                                                   
aCTACGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAtt                                                                                                                                                    >  1:581021/1‑150 (MQ=255)
  tACGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGc                                                                                                                                                  >  2:1717738/1‑150 (MQ=255)
    cGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAg                                                                                                                                                <  1:72994/150‑1 (MQ=255)
     gTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGc                                                                                                                                               <  1:26031/150‑1 (MQ=255)
     gTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGc                                                                                                                                               <  1:734361/150‑1 (MQ=255)
     gTCATCATCAATGATGATTTCGACACGACTCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGc                                                                                                                                               >  2:401219/1‑150 (MQ=255)
      tcatcatcaATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAGTCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGcc                                                                                                                                              <  1:1224231/150‑1 (MQ=255)
      tcatcatcaATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGcc                                                                                                                                              >  2:114229/1‑150 (MQ=255)
       catcatcaATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAATGCAGccc                                                                                                                                             <  2:1728078/150‑1 (MQ=255)
           atcaATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAAt                                                                                                                                         >  1:1761153/1‑150 (MQ=255)
              aaTGATGATTTCGACACGGCGCTGGCAGACTTGCGAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGt                                                                                                                                      >  1:1590637/1‑150 (MQ=255)
              aaTGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAATGCAGCCCCatcgat                                                                                                                                      <  1:246311/150‑6 (MQ=255)
               atgatgATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                                     <  1:1541465/150‑1 (MQ=255)
                    gatTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTc                                                                                                                                <  1:1717738/150‑1 (MQ=255)
                     atTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTcc                                                                                                                               >  2:1565456/1‑150 (MQ=255)
                     atTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTcc                                                                                                                               <  2:1615615/150‑1 (MQ=255)
                           acacGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGc                                                                                                                         >  1:695714/1‑150 (MQ=255)
                                gcgcTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCaa                                                                                                                    >  1:1369155/1‑150 (MQ=255)
                                       cAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGcc                                                                                                             >  1:1682353/1‑150 (MQ=255)
                                          aCTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAATGCAGCCCCat                                                                                                                                          <  1:1042642/118‑2 (MQ=255)
                                          aCTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAATGCAGCCCCat                                                                                                                                          >  2:1042642/1‑117 (MQ=255)
                                            ttGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAACAGAACGGCGGTAGGCCTGCAGCGGCAAGGTAGCctgctg                                                                                                        >  2:628281/1‑150 (MQ=255)
                                             tGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAAATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGa                                                                                                       <  2:1211532/150‑1 (MQ=255)
                                              gAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGAc                                                                                                      >  1:1734967/1‑150 (MQ=255)
                                                      tgtTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTa                                                                                              <  1:1749850/150‑1 (MQ=255)
                                                      tgtTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTa                                                                                              >  2:519998/1‑150 (MQ=255)
                                                            gCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCg                                                                                        >  2:911367/1‑150 (MQ=255)
                                                             cTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCTTCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGa                                                                                       <  1:628281/150‑1 (MQ=255)
                                                              tCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAg                                                                                      <  1:911367/150‑1 (MQ=255)
                                                              tCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAg                                                                                      <  1:156947/150‑1 (MQ=255)
                                                                gAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGa                                                                                    >  1:1695619/1‑150 (MQ=255)
                                                                 aaTCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCCCAGTCCGTTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGaa                                                                                   <  2:1590637/150‑1 (MQ=255)
                                                                  aTCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAg                                                                                                  >  1:1906007/1‑134 (MQ=255)
                                                                  aTCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAg                                                                                                  <  2:1906007/134‑1 (MQ=255)
                                                                     ggCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTAc                                                                               <  2:1734967/150‑1 (MQ=255)
                                                                     ggCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTAc                                                                               <  2:1813934/150‑1 (MQ=255)
                                                                      gctgctgctGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGTAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACg                                                                              >  2:1278781/1‑150 (MQ=255)
                                                                         gctgctgAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCg                                                                           <  2:1369155/150‑1 (MQ=255)
                                                                             ctgAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAg                                                                                             >  2:1115465/1‑128 (MQ=255)
                                                                             ctgAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAg                                                                                             <  1:1115465/128‑1 (MQ=255)
                                                                               gaagaaGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGcgcac                                                                     <  1:274270/150‑1 (MQ=255)
                                                                                aagaagCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                                     <  1:1819458/85‑1 (MQ=255)
                                                                                aagaagCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                                     >  2:1819458/1‑85 (MQ=255)
                                                                                aagaagCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGGTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGggcacg                                                                    >  2:125118/1‑150 (MQ=255)
                                                                                  gaagCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGcgcacgca                                                                  >  2:178322/1‑150 (MQ=255)
                                                                                   aagCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCAt                                                                 <  1:178322/150‑1 (MQ=255)
                                                                                   aagCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCAt                                                                 <  2:1695619/150‑1 (MQ=255)
                                                                                    agcagcagcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCa                                                                                                   <  1:1181944/115‑1 (MQ=255)
                                                                                    agcagcagcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCa                                                                                                   >  2:1181944/1‑115 (MQ=255)
                                                                                     gcagcagcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATgc                                                               >  1:1321996/1‑150 (MQ=255)
                                                                                     gcagcagcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATgc                                                               <  2:1397468/150‑1 (MQ=255)
                                                                                       agcagcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATgcgc                                                             <  1:519998/150‑1 (MQ=255)
                                                                                          agcagcGTCAGGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGCACGCCGGCGCACGCATGCCCGAc                                                          >  2:845413/1‑150 (MQ=255)
                                                                                           gcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGTTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTTCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACa                                                         <  1:667465/150‑1 (MQ=255)
                                                                                           gcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACa                                                         <  2:581021/150‑1 (MQ=255)
                                                                                           gcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACa                                                         <  1:401219/150‑1 (MQ=255)
                                                                                             agcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGc                                                       >  2:293859/1‑150 (MQ=255)
                                                                                              gcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGca                                                      <  1:788743/150‑1 (MQ=255)
                                                                                                 tCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTa                                                                                                                                   <  1:1107697/70‑1 (MQ=255)
                                                                                                 tCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTa                                                                                                                                   >  2:1107697/1‑70 (MQ=255)
                                                                                                   aTGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGCAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGcataagc                                                 >  2:932232/1‑146 (MQ=255)
                                                                                                    tGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCAtcatca                                                <  1:1477937/150‑1 (MQ=255)
                                                                                                        gCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCCGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGcc                                            >  2:1343361/1‑150 (MQ=255)
                                                                                                          aTTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGCCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGCAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGa                                          <  1:1036378/150‑1 (MQ=255)
                                                                                                              cTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGcc                                      >  2:1029279/1‑150 (MQ=255)
                                                                                                               tCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGcc                                                                                                                                              >  2:1184563/1‑45 (MQ=255)
                                                                                                               tCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGcc                                                                                                                                              <  1:1184563/45‑1 (MQ=255)
                                                                                                               tCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCa                                     >  1:544961/1‑150 (MQ=255)
                                                                                                                aaaGGAAATGATCGGCTGAAGATTGATGGGGCTGCATTGCCGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGCGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGCCCGAGGCCaa                                    <  2:1761153/149‑1 (MQ=255)
                                                                                                                 aaGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAg                                   >  2:1470996/1‑150 (MQ=255)
                                                                                                                 aaGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGAAGAGTAGAACCAAATAAGTACTCCCGGGCGCTCATGCGCGACAGCATCATCAGGCCGATGCCAAg                                   <  1:1470996/150‑1 (MQ=255)
                                                                                                                           tCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCga                             >  2:358261/1‑146 (MQ=255)
                                                                                                                           tCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCga                             <  1:358261/146‑1 (MQ=255)
                                                                                                                             ggCTGAAGAACGATGGGGCCGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAgatgat                       >  2:174642/1‑150 (MQ=255)
                                                                                                                              gCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATc                      >  2:1570872/1‑150 (MQ=255)
                                                                                                                                  aaGACCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCgg                  <  1:114229/150‑1 (MQ=255)
                                                                                                                                     aTCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCaa               >  2:1347379/1‑150 (MQ=255)
                                                                                                                                         ttGGGGCTGCAATGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGa                                          >  1:829334/2‑119 (MQ=255)
                                                                                                                                         ttGGGGCTGCAATGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGa                                          <  2:829334/118‑1 (MQ=255)
                                                                                                                                           ggggCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATcgcg                   >  2:208664/1‑140 (MQ=255)
                                                                                                                                           ggggCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATcgcg                   <  1:208664/140‑1 (MQ=255)
                                                                                                                                               cTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAacac     >  2:600507/1‑150 (MQ=255)
                                                                                                                                                tGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAACAcc    >  1:903781/1‑150 (MQ=255)
                                                                                                                                                  cATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAACAccag  >  1:167165/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACTACGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAACACCAG  >  NC_002947/6043574‑6043869

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