Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,193,996 T→G D70A (GAT→GCT)  thiE ← thiamine phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,193,9960TG95.5% 204.1 / NA 66D70A (GAT→GCT) thiEthiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base G (28/35);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCCTC  >  NC_000913/4193751‑4194241
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
gCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTg                                                                                                                                                                                                                                                  <  1:150693/251‑1 (MQ=255)
 ccGTGCCATCTGGTCCAGCCCCTGCGGTGCAGAAGGCGTCTGTTTGTTTTGCGTCGGGACCCCATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCTCAGCCGCAGTCTTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGCCGCCATACGCCTGATGCTTGAGCGCCAGCCGCCAGTAAGCGTTGa                                                                                                                                                                                                                                                 <  1:255832/251‑1 (MQ=255)
        aggggtcccAGCCCCGGCGGTGAGGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATAGAAGAGGGGCGGGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACaa                                                                                                                                                                                                                                          <  2:276994/244‑1 (MQ=255)
        aGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGCTTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCATGATCGCCCGCCGCCAGTACGCGTTGATAAACa                                                                                                                                                                                                                                           >  2:92719/1‑250 (MQ=255)
            tttCCAGCCCCTTCGGTGCATAAGGAGTCGTTTTGGTTTGCGGCGGGAACCCATGTCCCAGCGAGATACAAGCGGGGCGGGCTGCCGCCGCACCTTCGTTTTCCAGATCGTCATGTGGCGAAACGCCCAGCAGCAGGCCTGCCGCGCGGATGGCATTGAGCGCGGTGGCTTGCAAATCTTCCTGCCCAAAAGGGACGCCATACGCCTGATGCTGGATCGCCAGCCGCCAGTAATCGTTGATAAACAATcg                                                                                                                                                                                                                                       <  2:295531/249‑1 (MQ=255)
               ccAGCCCCTGCGGTGCAGAAGGCTTCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATcgcgcg                                                                                                                                                                                                                                   <  2:85733/251‑1 (MQ=255)
                 agccccctcGGTGCAGAAGGCTTCGGTTTGGGTTGCGTCGGGAACAAATGCCCAAGGCCGATAGAAGAGGGGCGTGCTGCCGGCGCGGCGTCGATTTCCATATCGTCATGGTTCGAACCGCCCAGCCGCAGGCCTGCCGCGCGGATGTCATTGAGATCGGTGGCTGGAAGATCTTCCTGCCCCAAATGGACGTCATACGCCGGATGCTTGATCGCCGGCGGCCAGTAATCGTTGATAAACAATCGCGCGtt                                                                                                                                                                                                                                 <  2:47160/243‑1 (MQ=255)
                         cGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATAACGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAgcggc                                                                                                                                                                                                                         >  1:132067/1‑251 (MQ=255)
                           gTGCAGAAGGCATTTGGGTGGTTTGCGGCGGGAAAACAGGTCACAGCGCGGTATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATGGGTGGCTTGCAAATCTTCCGGCCCCAAATGGACGCCATACGCCTGATGCTTGATGGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGcggcgg                                                                                                                                                                                                                       <  1:12704/251‑1 (MQ=255)
                                 aaGGCATCTTTTTGGTTGGCGTCGGGAACACATGTCCCAGCGCTATATAAGAGGGTCGGTCCGCCAGCGCGCCGTTGATTCCCATAGCGTCATGTGTCGAACCGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATGCGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTGGATAAACAATCGCGCGTTATAgcggcgttccagc                                                                                                                                                                                                                 <  2:130207/251‑8 (MQ=255)
                                      aTCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAAt                                                                                                                                                                                                            >  1:66818/1‑251 (MQ=255)
                                         tGTTTGGTTTGCGGCGGGAACACATGTCCCAGCGCGATAAAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGAGGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTg                                                                                                                                                                                                          <  2:22463/250‑1 (MQ=255)
                                          gtttggtttgCGTCGGGACCACATGCCCCAGCGCTATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATGTCTTCATGTGTCTAAACGCCCAGCCGCAGGCCTGCCGCGGGGATGGCATTTAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCTTACGCCTGATGCTTGATCGCCAGCCGCCATTAATCGTTGATAAACAACCGCGCGTTATAGCGGCGGCCCAGCGCAATTgc                                                                                                                                                                                                         <  1:295365/250‑1 (MQ=255)
                                           tttggtttgCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCCATTgccc                                                                                                                                                                                                       >  2:248634/1‑250 (MQ=255)
                                                   gcgccGGGAACCCATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGACCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACa                                                                                                                                                                                               <  1:241282/247‑1 (MQ=255)
                                                     gTCGGGCACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCACGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATc                                                                                                                                                                                             >  1:233762/1‑251 (MQ=255)
                                                       cGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCgg                                                                                                                                                                                           <  1:9027/251‑1 (MQ=255)
                                                          gAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCCCGACATCGGCtt                                                                                                                                                                                        >  1:219497/1‑251 (MQ=255)
                                                                    cccAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAAGCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTc                                                                                                                                                                              <  1:113597/251‑1 (MQ=255)
                                                                      cAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCAt                                                                                                                                                                            >  1:310800/1‑251 (MQ=255)
                                                                        gcgcGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATcg                                                                                                                                                                          <  1:16777/251‑1 (MQ=255)
                                                                           cgATATAAGAGGGGCGTGCTGCCAGCGCGCCGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCgcgc                                                                                                                                                                       >  2:182108/1‑251 (MQ=255)
                                                                                taAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCACCCCGCCGATc                                                                                                                                                                  >  1:257771/1‑251 (MQ=255)
                                                                                    gggggCGTGCTGCCATCGCGACGTCGATTCCCTAACCTCCAGGTGCCAAAACGCCCGGCCGCAGGCCGGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACTCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTg                                                                                                                                                              <  2:220784/250‑1 (MQ=255)
                                                                                        gCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGc                                                                                                                                                          >  1:264652/1‑251 (MQ=255)
                                                                                          gTGCTTCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGACCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGt                                                                                                                                                        <  2:163015/251‑1 (MQ=255)
                                                                                          gTGCTGCCAGCGCGACGTCGATTTCCCTATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGcc                                                                                                                                                                                                                     >  2:39296/1‑190 (MQ=255)
                                                                                          gTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCGTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGcc                                                                                                                                                                                                                     <  1:39296/190‑1 (MQ=255)
                                                                                             ctgcCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCCGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGc                                                                                                                                                     >  1:132649/1‑251 (MQ=255)
                                                                                               gcCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTg                                                                                                                                                   >  1:8959/1‑251 (MQ=255)
                                                                                                       gACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAgtac                                                                                                                                           >  1:189146/1‑251 (MQ=255)
                                                                                                         cGTCGATTTCCATCTCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGCGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCCGCGCAATTGCCGCCCCGACCTCGGCTTCCCCCTCTTCATCGCGCCCCTCTTTGATGCGTAGCTGGAGAgtacgt                                                                                                                                         >  2:179117/1‑251 (MQ=255)
                                                                                                             caTTTCCATAACGTCAGGTGTCAAAACGCCCAGCCGCAGGCCTGCCGCTCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAACGGACGCCATACGCCTGATGCTTGATCCCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCTCCACGACATCGTCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGc                                                                                                                                     <  1:302728/250‑1 (MQ=255)
                                                                                                                  ccATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGc                                                                                                                                 >  1:7575/1‑250 (MQ=255)
                                                                                                                    atatCGGCATGTGTCGAAACGCCCAGCCGCAGGCCAGCCGCGCGGATGGCATTGAGATCGGTGGATTGCAAATCTTCCTGCCCCAAAAAGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATc                                                                                                                              >  1:24701/1‑251 (MQ=255)
                                                                                                                    atatCGACGTGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGCTCGGTGGCTTGCAAATCTTCCTGCCCCAAATTGACGCCATACGCTTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCTCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCGCGACATCGGCGTCCGCCTCGGCATCGCGTCGATCTTTGATGCGTAGCTGGAGAGGACGTACGCCTGCATc                                                                                                                              <  1:132913/251‑1 (MQ=255)
                                                                                                                             cccgatcaaaaCGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTTCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGAc                                                                                                                      <  1:34346/243‑1 (MQ=255)
                                                                                                                                ggtggaaccGCCCAGCCGCAGGCCGGCCGCGCGGATGGCATTGAGATCGGTGGCTTACAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTGGACCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCGTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGtt                                                                                                                   <  2:146471/243‑1 (MQ=255)
                                                                                                                                 gtCGAAACGCCCGGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGAACGCCTGCATCCAACAGACGTTc                                                                                                                  <  1:110691/251‑1 (MQ=255)
                                                                                                                                   cgaacagCCCAGCCGGATGCCTGCCGCGCGGATGGCATTGAGATGGGTGGCTCGCAAGTCTTCCTGCCCCAAAGGGACGCCATACGCCTGAGGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTAGAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTACCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCGGGAGAGTACGTACGCCTGCATCCAACAGAGGTTCg                                                                                                                 <  2:266068/244‑1 (MQ=255)
                                                                                                                                     aaaCGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGGGAGTACGTACGCCTGCCTCCAACAGAc                                                                                                                      >  2:34346/1‑243 (MQ=255)
                                                                                                                                             gccgtcAGGCCCGCCGCGCGGTTGGTTTTGATATCGGGGGCTTGCAATTCTTCCTGCCCAAACTGGACGCCATACTCCTTGTGCTGGAGCGCCAGCCGCCAGTTAGCGTTGATAACCAATCGCGCGTTATATCGGCGGCCCCGCGCAATTGCCGACACGACTTCGGCTTCCACCTCTTCATCGCGCCGATCTTTGTTGCGTAGCTGGAGAGTACGTACGCCTGCATCCATCAGTCGTTCGATCCACTGTAc                                                                                                      <  2:244414/246‑1 (MQ=255)
                                                                                                                                                 gCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTg                                                                                                  <  1:182108/251‑1 (MQ=255)
                                                                                                                                                     gCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACg                                                                                                     >  2:154281/1‑244 (MQ=255)
                                                                                                                                                     gCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACg                                                                                                     <  1:154281/244‑1 (MQ=255)
                                                                                                                                                          ccGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCaccacc                                                                                         <  1:248634/251‑1 (MQ=255)
                                                                                                                                                             cgcgGATGTCAGTGAGATCGGGGGCTTGCAAATCTCCCCCCCCCAATTGGAACCCATACCCCGGATGCTTGATCGCCAGCCGCCAGTAAGCCGTGAGAAACAACCGCGCGTTATAGCGGCGGCCCAGCGCAATGGCCGCCACGACATCGGCTTCCACCGCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGTCGTTCGATCCACTGTACGCTGTCCACCTCCggg                                                                                      <  2:132067/251‑1 (MQ=255)
                                                                                                                                                              gcgGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGCGAGTACGTACGCCTGCCTCCAACAGACGTTCGCTCCACTGTACGCTGTCCACCACCgggg                                                                                     >  2:57977/1‑250 (MQ=255)
                                                                                                                                                               cgGATGGCATTGAGATCAGAGGCTTTCAAACCTTCATGCCCCAAATGGAAGCCATACGCCTGATGCTTGCTCGCCAGACGCCAGTAAGCGTTGATAAACAATCGCGCGTTCTAGCGGCGGCCAAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTTATGCGTCGATGGTGAGCACGTACGCCTGAATCCCACACACGTTCGATACCCTTAACGCTGTCCACCACCgggc                                                                                     >  2:45972/1‑249 (MQ=255)
                                                                                                                                                                       aTTGAGCTCGGAGGATTGCCACTCTTTCTGCCCCAAATGGACGCCATACGCCTGATGATTGAACGCCAGCCGCCAGTAAGCGTTGCTAAACAATAGCGAGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATTCGTCGCTGGCGAGTACGTAAGCCTGCATCACACAGACGTTCGATCCACTGTACGCTGTCCAACACAGGGTACCGTCCTg                                                                            >  2:47997/1‑251 (MQ=255)
                                                                                                                                                                          gaaaTCGGTGGCTTGCAAACCTTCCTCCCCAAAATGGACGCCATACGCCGGATGTTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAAc                                                                         <  2:189146/248‑1 (MQ=255)
                                                                                                                                                                           gtatgggggctgtgaaaagattCCTGCCCCAAAGGAACCCCTTACGCCGGATGCTTGAGCGCCAGCCGCCAGTAAGCGTTGATAACCAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGGTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACg                                                                        <  2:52086/231‑1 (MQ=255)
                                                                                                                                                                             aTCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGa                                                                       >  2:296143/1‑250 (MQ=255)
                                                                                                                                                                                  tGGCTTGCAAATCTTCCTGACCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCCGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGAAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAATGTACGCTGTCCACCCCCggggagcatcctgcaagaaaaggg                                                                 >  2:105159/1‑230 (MQ=255)
                                                                                                                                                                                        gCAAAGCTTCCTGCCCCAAAGGGACGCAATACGCCTGTTGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCGTCCACCTCTTCACCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGGACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACagga                                                           <  2:257771/251‑1 (MQ=255)
                                                                                                                                                                                            aTCTTCCTGCCCCAAATGGACGCAATAAGCCTCCTCCTTGATCGCCAGCCGCCAGTAAGCGTTAATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCCTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCCCCGGGGCCAATCCTGAACGAAACGGGACCggagga                                                        >  2:176326/1‑250 (MQ=255)
                                                                                                                                                                                                tccggCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGGACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGaaa                                                      <  1:176326/244‑1 (MQ=255)
                                                                                                                                                                                                                cGCCATACGCCGTATGCTTGACCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACg                                   <  2:132649/251‑1 (MQ=255)
                                                                                                                                                                                                                    aTACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTtc                               <  1:248073/251‑1 (MQ=255)
                                                                                                                                                                                                                      aCGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTtcct                             >  1:134672/1‑251 (MQ=255)
                                                                                                                                                                                                                      aCGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTtcc                              <  1:57977/250‑1 (MQ=255)
                                                                                                                                                                                                                          cTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCtt                         >  1:245591/1‑251 (MQ=255)
                                                                                                                                                                                                                                      cGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGAt             <  1:42100/251‑1 (MQ=255)
                                                                                                                                                                                                                                                cAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAggagga                                                        <  1:169591/198‑1 (MQ=255)
                                                                                                                                                                                                                                                cAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGAACAggagga                                                        >  2:169591/1‑198 (MQ=255)
                                                                                                                                                                                                                                                 aGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCctc  <  1:269299/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
GCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCCTC  >  NC_000913/4193751‑4194241

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