Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,491,704 G→T 5.9% intergenic (‑84/+21) yhfG ← / ← ppiA putative antitoxin for Fic/peptidyl‑prolyl cis‑trans isomerase A (rotamase A)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,491,7040GT 5.9% 237.2 / 3.2 95intergenic (‑84/+21)yhfG/ppiAputative antitoxin for Fic/peptidyl‑prolyl cis‑trans isomerase A (rotamase A)
Reads supporting (aligned to +/- strand):  ref base G (25/63);  new base T (2/3);  total (29/66)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.73e-01

CTGAAAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAA  >  NC_000913/3491555‑3491852
                                                                                                                                                      |                                                                                                                                                    
cTGAAAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAg                                                                                                                                                      <  1:1242407/150‑1 (MQ=255)
 tGAAAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGa                                                                                                                                                     >  1:110047/1‑150 (MQ=255)
  gAAAATTACGCTTACGCGGTAGCTCCGAGAGAGGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGAt                                                                                                                                                    <  1:269792/150‑1 (MQ=255)
  gAAAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCTGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGGAAGAt                                                                                                                                                    >  1:905315/1‑150 (MQ=255)
    aaaTTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTg                                                                                                                                                  <  2:216113/150‑1 (MQ=255)
    aaaTTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTg                                                                                                                                                  <  1:2899/150‑1 (MQ=255)
     aaTTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGc                                                                                                                                                 >  2:221406/1‑150 (MQ=255)
       ttacgattacgCTGAAGCTCCCAGAGACGGAACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGccc                                                                                                                                               <  2:1623213/150‑1 (MQ=255)
        tacgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCcg                                                                                                                                              >  1:314579/1‑150 (MQ=255)
          cgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCCGGGAAAAATTcccccc                                                                                                                                             >  2:1660665/1‑148 (MQ=255)
          cgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAATCTTGCCcgcg                                                                                                                                            <  1:378250/150‑1 (MQ=255)
          cgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTACACTTGATGCCTCCCTGAACGTCTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCcgcg                                                                                                                                            >  1:1300576/1‑150 (MQ=255)
          cgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGATTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTTCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAATCTTGCCcgcg                                                                                                                                            <  2:1104574/150‑1 (MQ=255)
            attacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCT‑CGGGCAATCTTGCCcgcgct                                                                                                                                          <  1:1037094/150‑2 (MQ=255)
            attacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCcgcgcg                                                                                                                                          <  2:228941/150‑1 (MQ=255)
              tacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCgag                                                                                                                                        <  2:429793/150‑1 (MQ=255)
              accgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCCCAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCgag                                                                                                                                        <  1:1776276/148‑1 (MQ=255)
               acgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGGGAGTGGCTTCACTGGATGCCTCCCTTACCGTGCCCATTTCCCCCAATTATAAGCGGCAAATCCGGGTAACGCCGGGAGAGCATAAGCG‑CGGGCAAGATTGCCCGCGCgaga                                                                                                                                       <  1:1703665/150‑1 (MQ=255)
                   tGAAGCTCCCAGACCCGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTCCCACAAGTATAATCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATc                                                                                                                                   <  1:496036/150‑1 (MQ=255)
                   ccaaGCTCCCAGAGACGGGACTTTTGCAAATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATATGCCCCAAGAATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATc                                                                                                                                   <  1:1898870/148‑1 (MQ=255)
                    gAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCa                                                                                                                                  <  2:20444/150‑1 (MQ=255)
                    gAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCa                                                                                                                                  <  1:1024200/150‑1 (MQ=255)
                    gAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATACGGGTTACGCGGGGAGATCAGAAGCG‑CGTTCAAGCTTGCCCGCGCGATAACTCa                                                                                                                                  >  2:278604/1‑150 (MQ=255)
                     atgcgcccAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGACCATTTGCCACAATTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCAt                                                                                                                                 <  1:922388/145‑1 (MQ=255)
                     aaGCTCCCAGAGCCGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTTATGCCTCCCTGAACGTGCCCATTTGCCACAAGAATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCAt                                                                                                                                 <  1:1592048/150‑1 (MQ=255)
                      aGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCAtt                                                                                                                                <  2:529839/150‑1 (MQ=255)
                      aGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGGCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTGACGCCGGGAGAGCAGAAGCG‑CGGGCAAGACTGCCCGCGCGAGAAATCAtt                                                                                                                                <  1:1540235/150‑1 (MQ=255)
                             agagACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAAATTTCCCGCGCGAGAAATAATTACGGCAg                                                                                                                         >  1:12542/1‑150 (MQ=255)
                                    ggACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCGAGATTTACCGCGGCAGCAATCATTTCGGCgggccgttt                                                                                                                  >  1:907583/1‑144 (MQ=255)
                                     gACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAg                                                                                                                 <  1:473583/150‑1 (MQ=255)
                                     gACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAg                                                                                                                 <  1:1847590/150‑1 (MQ=255)
                                        ttttGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCgg                                                                                                              <  2:1647129/150‑1 (MQ=255)
                                            gCTTATCGGTGTGTTTATACTCTTGATGCCTCCCTGAACGTTTCCATTTGCCACATCTCTAAGCGGCTAATCCGGGTTACGCAGGGAGAGCAGAAGCG‑CGGGCAAGTTTGCCCGCGCGAGAAATCATTACGGCAGTACTTTAGCGGAAAg                                                                                                          <  2:1388411/150‑1 (MQ=255)
                                            gCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAATCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAg                                                                                                          <  2:1550876/150‑1 (MQ=255)
                                            gCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAg                                                                                                          <  1:1524231/150‑1 (MQ=255)
                                             cTTATCGGTGAGGTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGACGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAgg                                                                                                         <  1:227973/150‑1 (MQ=255)
                                              ttATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGa                                                                                                        <  2:1579494/150‑1 (MQ=255)
                                              ttATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTACCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGa                                                                                                        >  2:1368098/1‑150 (MQ=255)
                                               tATCGGTGAGTTTCTTCACTTGCTGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGCTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCCCGAGAAATCATTACCGCCGGACTTTTGCGGAAAGGaa                                                                                                       >  2:991956/1‑149 (MQ=255)
                                               tATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CCGGCAAGCATGGCCGCCCGGCAAATAATTAATGGCGGCCTTTAGGGGAAAgggt                                                                                                       >  1:1867058/1‑148 (MQ=255)
                                                  gggTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAAc                                                                                                    <  2:1601288/149‑1 (MQ=255)
                                                  cGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAAc                                                                                                    <  1:605403/150‑1 (MQ=255)
                                                  cGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAAc                                                                                                    <  2:1342982/150‑1 (MQ=255)
                                                      gAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGGGAGAAATCATTACGGCAGGACGTTTGCGGAAAGGATAACGAcc                                                                                                >  2:94801/1‑150 (MQ=255)
                                                       aGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCg                                                                                               <  2:486727/150‑1 (MQ=255)
                                                        gTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATATGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGCAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCAATACGGCAGGACTTTAGCGGAAAGGATAACTGCCgg                                                                                              >  1:1222643/1‑150 (MQ=255)
                                                          ttcttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGtt                                                                                            <  1:1076601/150‑1 (MQ=255)
                                                            cttcACTTGATGCCGCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGtttt                                                                                          <  1:602319/150‑1 (MQ=255)
                                                               cACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCAATACGGCAGGACTTTAGCGGAAAGGAAAACTACCGGTTTTGAc                                                                                       >  2:14398/1‑150 (MQ=255)
                                                                   tGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGca                                                                                   <  2:314579/150‑1 (MQ=255)
                                                                   tGAAGCCTCCCTTAACGTGCCCATTTGCCACAAGTATCAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGca                                                                                   <  2:1143318/150‑1 (MQ=255)
                                                                           cccTGAACGTGTCCTTGTGCCTCAAGTCTAAGCGGCAAATCCTGGTTACGCCGGGAGAGCAGAAGCG‑CTGGCAAGATTGCCCGCGCGAGAATTCATTAAGGCAGGACTTTAGCGTAAAGGATAACTACCGGTTTTGACGGCACATTCTgg                                                                           <  1:1176467/150‑1 (MQ=255)
                                                                           cccTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CTGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTgg                                                                           <  1:1175320/150‑1 (MQ=255)
                                                                           cccTGAACGTGGCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTgg                                                                           <  2:877065/150‑1 (MQ=255)
                                                                           cccTGAACGGGGCCATTTGCCACCAGTATAAGCGGCAAATCCGGGCTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGACAGGATAACTACCGGTTTTGACGGCACATTCTgg                                                                           <  1:519440/150‑1 (MQ=255)
                                                                               caaCGTGTCCAATTGCCACAAGTAAAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACg                                                                       <  2:1670556/149‑1 (MQ=255)
                                                                                 aCGTGTCCATTTGCCACAAGTATAAGCGGCAACTCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGa                                                                     <  2:86706/150‑1 (MQ=255)
                                                                                    tgtCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTTGACGCGCGACAAAAAAATTATGCAGGTCATTATGGGAGACGATAACAACCCGGTTTGACGGCCCCTTGAGGTACgggacc                                                                  >  2:924721/1‑146 (MQ=255)
                                                                                    tgtCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCa                                                                  >  1:1359855/1‑150 (MQ=255)
                                                                                           ttGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCAt                                                           <  2:845286/150‑1 (MQ=255)
                                                                                             gCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCCAGATTTCCCGCGCGAGAAAACATTACGGCAGGTCATTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGa                                                         >  1:1258614/1‑150 (MQ=255)
                                                                                              ccACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAAGTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAg                                                        <  2:104893/150‑1 (MQ=255)
                                                                                               cacaAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGGTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGt                                                       <  2:389697/150‑1 (MQ=255)
                                                                                                  caGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTGTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCgg                                                    <  1:1406631/149‑1 (MQ=255)
                                                                                                   aGTATAAGCGTCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGTGCAAGATTGCCCGCGCGAGAAATCATTACGTCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGc                                                   <  2:1338989/150‑1 (MQ=255)
                                                                                                        aaGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTg                                              <  2:220315/150‑1 (MQ=255)
                                                                                                               aaaTCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGATAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATc                                       <  2:1005198/150‑1 (MQ=255)
                                                                                                                aaTCCGGGTTACGCCGGGAGCGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCt                                      <  1:377545/150‑1 (MQ=255)
                                                                                                                aaTCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCt                                      <  2:1106262/150‑1 (MQ=255)
                                                                                                                    cGGGTTACGCCTGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATc                                  <  2:980923/150‑1 (MQ=255)
                                                                                                                    cGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATc                                  <  1:1102692/150‑1 (MQ=255)
                                                                                                                     gggTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATTAGTCGGCACCTGGGAAATCTTATCg                                 <  2:1300576/150‑1 (MQ=255)
                                                                                                                     gggTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCg                                 <  1:221406/150‑1 (MQ=255)
                                                                                                                           cGCCGGGAGAGCAGAAGCG‑CGGGCCAGATTGCCCGCGCGAGAAATCATTACGGCAGGGCTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGGCATGAGTCGGCCCCGGGGAAATCTTTTCGGCAACg                           >  1:1262305/1‑150 (MQ=255)
                                                                                                                                   gagCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACATTAGCGGAAAGGATAACTACCGGTATTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATAGGCAACGTCCATGcc                   >  1:68182/1‑150 (MQ=255)
                                                                                                                                    atttctcgcg‑cGGGCAATCTTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCt                  >  2:1064663/8‑150 (MQ=255)
                                                                                                                                      cAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCttt                >  2:1259365/1‑150 (MQ=255)
                                                                                                                                       aGAATCG‑CGGGCAAGATTTCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTc               <  2:905315/150‑1 (MQ=255)
                                                                                                                                       aGAAGCG‑CGGGCAAGAGTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTc               <  1:1344244/150‑1 (MQ=255)
                                                                                                                                        gAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTca              <  2:399166/150‑1 (MQ=255)
                                                                                                                                           gcg‑cGGGCAATCTTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGATTTGCCGGCACATTCTGGTACGGCCCAACGGCAGGAGTCGGCACCTGGGAAATGTTATCGGCAACTTCAATGCCTTTcccca           >  2:1340839/1‑150 (MQ=255)
                                                                                                                                            cg‑cgGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTcaccac          <  2:1772453/150‑1 (MQ=255)
                                                                                                                                             g‑cgGGCATGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACt         <  2:1716286/150‑1 (MQ=255)
                                                                                                                                             g‑cgGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACt         <  2:524017/150‑1 (MQ=255)
                                                                                                                                               cgGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACtt        >  1:996277/1‑150 (MQ=255)
                                                                                                                                               cgGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACATGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACtt        >  1:462120/1‑150 (MQ=255)
                                                                                                                                                gggCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACttt       >  1:1121146/1‑150 (MQ=255)
                                                                                                                                                gggCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACttt       >  1:1396995/1‑150 (MQ=255)
                                                                                                                                                gggCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACttt       >  1:159984/1‑150 (MQ=255)
                                                                                                                                                gggCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGAGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACttt       >  1:1724848/1‑150 (MQ=255)
                                                                                                                                                gggCAAGATTGCACGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGCATAACTACCGGTTTTGACGGCACATTCTGGCACTGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACttt       <  1:1164349/150‑1 (MQ=255)
                                                                                                                                                   cAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAcc    <  2:1567636/150‑1 (MQ=255)
                                                                                                                                                   cAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAcc    <  2:907583/150‑1 (MQ=255)
                                                                                                                                                    aaGATTGCCCGCGCGATAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTAATTGTTTTGACAGTACATTCTGGTACGGACCAAAGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCAAGCCTTTCCCCACTAAACCa   <  2:290734/150‑1 (MQ=255)
                                                                                                                                                    aaGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCa   <  1:219508/150‑1 (MQ=255)
                                                                                                                                                     aGATTGCCCGCGCGAGAAATCATTACTGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTTGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCaa  <  2:110047/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CTGAAAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAA  >  NC_000913/3491555‑3491852

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