Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913825,2350TA30.5% 76.8 / 21.5 61M292L (ATG→TTG) ybhRputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (27/14);  new base A (3/15);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 6.31e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATGAAGATCTGCC  >  NC_000913/825003‑825462
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                    
tACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGGATCCTAAAGAc                                                                                                                                                                                                                                >  1:75408/1‑239 (MQ=255)
 aCATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGGGAATGGCGGGGACCAAAAAAACa                                                                                                                                                                                                                               >  1:593985/1‑239 (MQ=255)
       ccTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCCAAAAGACAAaaac                                                                                                                                                                                                                          >  1:1589613/1‑238 (MQ=255)
                 aaaCATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCAAAAAGACAAACACCCC‑GCTAAAc                                                                                                                                                                                                               >  1:630218/1‑238 (MQ=255)
                       cgcgTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGCCTCGCATCCTTCAAATAAAT‑CTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAAAACACC‑GGTAAAAGCct                                                                                                                                                                                                           >  1:111865/1‑236 (MQ=255)
                                     cctgtGTCGTGGCCGTAACCACCAGTGGCGACCACAAACTATACCCAAAATTATCCGTACCGGCATCATT‑AAATAAATGCTGCGTGGTAATGTCCGAAAAGGGACAAATAGGGTTAAACCAGGTCGGTTTTGGCAGCCGTACCGGCATTTTTCCCACCGGAAAACCGTAACAGAAAAGGAGAAGGGGGGGCATAATAAAGACAACCACGCC‑GATAAACGCCTGCTGTTGTGTTGAACa                                                                                                                                                                                           <  2:111865/234‑1 (MQ=255)
                                               ggCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCATGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCCTAAAGACAAACACCCC‑GAAAAAACCCTTCTTTTGGGTTGAAAAGAGggagaaaa                                                                                                                                                                                >  1:2089221/1‑234 (MQ=255)
                                               ggCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCACCATAAAGACAAACAAGCC‑CCAAAACGCCTGATGTTGTGTTGTACAAAGGGGTGaaa                                                                                                                                                                                >  1:1444758/1‑239 (MQ=255)
                                                gCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAAAACGCC‑GATAAAAGCCCGCTGTTGTGTTGAACAGAGGGATGaaa                                                                                                                                                                                >  1:499899/1‑238 (MQ=255)
                                                    ttATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAAAACGCC‑GATAAACGCCTGCTGTTGTGTTGAACACAGTGATGAAaacaac                                                                                                                                                                           >  1:1576129/1‑239 (MQ=255)
                                                              tAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCC‑GATAAACGCCCGCTGTTGTGTTGAACAAAGGGATGAAATCAACAAACCGAAt                                                                                                                                                                  >  1:329016/1‑238 (MQ=255)
                                                               agcgtgaacaAACGACTCCAAACAATACCCGTACTCGTATCCTTAAAATAGAT‑CTGCTGGGTAATGACCTTAAAGTGGCGAAGAGGGTTATTCCACTGCATGTTTAGCATCCGTCGCGGCATGTTTCCCACCGGCGAAACGGAACTGAAAAGAAGATGGGCGGTAATGAGAAAAAAAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGCCCGaa                                                                                                                                                                   <  2:1505664/229‑1 (MQ=255)
                                                                              ttCCACAAAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGTTTGGAAATTTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTACCCGGAAAGGAGAATGGCGGGCATCAAAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGaa                                                                                                                                                                                 <  1:842079/207‑1 (MQ=255)
                                                                                       ttatcagatcTCGCATCTTGAAAATATAT‑CGGCTTGGTAATGTCCGAAAGGGGGCGAATAGGGTAAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTACCCCCGGAGAAACGTAACGGAAAAGGAGAATGGCGGGCAAAAAAAAGAAAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCAc                                                                                                                                        <  2:630218/230‑1 (MQ=255)
                                                                                       ataccaaGACTCTCATCCTTCAAATAAAT‑CTGCTTGTTAAGCGCCGCAAAGTGTCGAATATGGTTACTCCACGGCAGGTGGTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCAc                                                                                                                                        <  2:1139591/233‑1 (MQ=255)
                                                                                          tCCAGCCTCGCATCCTTAAAATAAAT‑CTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCCCCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGAAAAACACGCC‑GATAAACGCctgctg                                                                                                                                                                                                       <  1:749013/173‑1 (MQ=255)
                                                                                                    cATCCTTCAAATAAAT‑CTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGGATCATAAAGGCAAAAACCCCGGATAACCCCCGGC‑GTTGTGGTGaa                                                                                                                                                                                             >  2:749013/1‑173 (MQ=255)
                                                                                                             aaaaaat‑CTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCACCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGAAAAGGAGAAGGGCGGGCATCATAAAGACAAACACGCC‑GATAAACGCCTGCtgttgtg                                                                                                                                                                                                  <  1:1640131/156‑1 (MQ=255)
                                                                                                                        cTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGGATCATAAAGACAAAAACGCC‑GAGAAAACGCTGCTGTTGTGTTGAACagag                                                                                                                                                                                        >  2:1640131/1‑159 (MQ=255)
                                                                                                                            gTAATGTCCGTAAGGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATTTTTTCCACCGGAAAAACGTAACCGGAAAGGAGAATGGCGGGCTTCAAAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGaa                                                                                                    <  1:314979/239‑1 (MQ=255)
                                                                                                                                  tCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGAAAAAAACGCC‑GATAAACCCCTGCtgttgtgtt                                                                                                                                                                                                >  1:176404/1‑141 (MQ=255)
                                                                                                                                  tCCGTAAAGGGGCAAATAGGGTTAACCCACGTCAGTTTTTGCAGCCATCCCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAAGGGGGGGCATCATAAAGACAAACACGCC‑GATAAAc                                                                                                                                                                                                               <  2:148304/126‑1 (MQ=255)
                                                                                                                                      tAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAAAAGCGCCAGCGATCCGGCGTAGGGTATTTGa                                                                                          >  1:1498018/1‑239 (MQ=255)
                                                                                                                                            ggCGAATGGGTTTAATCCACGTCAGTTTTTGCAGCCATCCCGGCATGTTTTCCCCCGGAAAAACGTAACCGGAAAGGAGAATGGCGGGCATCAAAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATa                                                                                        <  1:1228598/235‑1 (MQ=255)
                                                                                                                                            ggCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCC‑GATAAACGCCTGCtgtt                                                                                                                                                                                                     >  1:148304/1‑126 (MQ=255)
                                                                                                                                             gggAATAGGGTGAAACCCGTTGAGGTTTTGGAGCAATACGGGCATGTTTACCACCGGAGAAACGTAACCGAAAAGGGGAAGGGGGGGAATAAAAAAGAAAAACACGCC‑GATAAACGCctg                                                                                                                                                                                                          <  2:1245331/118‑1 (MQ=255)
                                                                                                                                              cGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGGATCATAAAAACAAAAACCCA‑GCAAAACACCCGTCGTGTGGTTTAACACAGGGATGAAAAAAAACAACAGAAACCCACCACCCATAAACCAAAAAACAACAACCTAAACTAGAAAACAGCCAGCCATCCGGCGAAGGGGACTTTATACgaaca                                                                                  >  2:977816/1‑235 (MQ=255)
                                                                                                                                                      gtaaacacACGTAATGTTGTCCAGCCATACCGGCATGTTTTCCACCGGAGAACCGTAACCGGAAAGGAGAATGGGGGACATCAAAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATGCGCCCCGATACCa                                                                           <  2:1589613/231‑1 (MQ=255)
                                                                                                                                                       ttAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGAc                                                                                                                                                                       >  1:1428024/1‑145 (MQ=255)
                                                                                                                                                       ttAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCC‑CATAAACGCCCTCtgttgggtt                                                                                                                                                                                                >  1:1245331/1‑120 (MQ=255)
                                                                                                                                                          aTCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCg                                                                       >  1:436529/1‑238 (MQ=255)
                                                                                                                                                               cGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGcatcat                                                                                                                                                                                                                                      <  2:1717168/75‑1 (MQ=255)
                                                                                                                                                               cGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGcatcat                                                                                                                                                                                                                                      >  1:1717168/1‑75 (MQ=255)
                                                                                                                                                               cGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATa                                                                                                                                                                                                                                     <  1:236980/76‑1 (MQ=255)
                                                                                                                                                               cGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATa                                                                                                                                                                                                                                     >  2:236980/1‑76 (MQ=255)
                                                                                                                                                               cGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGaa                                                                                                                                                                                             >  2:895367/1‑115 (MQ=255)
                                                                                                                                                               cGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGaa                                                                                                                                                                                             <  1:895367/115‑1 (MQ=255)
                                                                                                                                                               cGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGGGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGc                                                                  >  1:1271939/1‑238 (MQ=255)
                                                                                                                                                               cGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGGGATGAAATCAACAGAACGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATAACAATCGCCAGca                                                                 >  1:1967350/1‑239 (MQ=255)
                                                                                                                                                               cGTCAGGTTTGGCAGCCTACCCGGCATGTTTCCCACGGGAGAAACGTAACGGGAAAGGAGAATGGCGGGAATCAAAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGca                                                                 <  1:284135/239‑1 (MQ=255)
                                                                                                                                                                      ttttGCAGCAATCCCGGATTGTTTCCCACCGGAAAAACGTAACCGAAAAGGAAAATGGCGGGAATCAAAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAAtg                                                            <  1:834980/237‑1 (MQ=255)
                                                                                                                                                                      ttttGCAACCATACCGGAGTGTTTTCCACCGGAAAAGCGTAGCCGGAAAGGAGAATGGCGGGCATCATAAAAAAAAACACGCC‑GATAAACGCCTGCGGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAg                                                                                                                      <  2:1011943/179‑1 (MQ=255)
                                                                                                                                                                              ccATACCGGCATTTTTTCCACCGGAGAAACGTACCCGGAAAGGAGAAGGGGGGGCATCAAAAAGACAAACACGCC‑GATAAACGCCTGCtgttgt                                                                                                                                                                                                   <  1:1423391/94‑1 (MQ=255)
                                                                                                                                                                                aTACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAAACAAACCCCCC‑GATAAACGCCTGCTGTTTTGTTGAACAGAGGGATGAAATCAACAGACCGAATCCCACCACCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGAt                                                                                                            >  1:1011943/1‑179 (MQ=255)
                                                                                                                                                                                        aTGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAAACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACaga                                                                                                                                                                                         >  2:1423391/1‑94 (MQ=255)
                                                                                                                                                                                           ttttCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCACCATAAAGAAAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGcc                                                                                    <  1:475241/192‑1 (MQ=255)
                                                                                                                                                                                           ttttCCACCGGAAAACCTTACCCGAAAAGGAGAATGGCGGGCATCAAAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATa                                                                                        <  1:2064171/188‑1 (MQ=255)
                                                                                                                                                                                                   cgaggaAACGTACCCGGAAAGTAGAATGGCGGGCATCAAAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGAt                                                                                              <  1:1284515/170‑1 (MQ=255)
                                                                                                                                                                                                   cGGAGAAACGAACCGGGAAGGGAAAAGGGGGGGAAAAAAAAAAACAACCACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGc                               <  1:752148/237‑1 (MQ=255)
                                                                                                                                                                                                     gagaAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAAAACGCC‑GATAAACGCCTGCTGTTGTGTTTAACAGAGTGATGAAATCAACAGACCGAATCCCCCCAGCGATAAAACATAAAATACCCACGGAAAGGAGAACAACGCCCGGGGTCCGGGGGAGGGGATTTGATACGCCCAGATa                                                                              >  2:2064171/1‑188 (MQ=255)
                                                                                                                                                                                                         aaCGAAACGGAAAGGGAGAAGGGCGGTCATCAAAAAGAAAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAAc                                                                                                                                                  <  2:1349778/116‑1 (MQ=255)
                                                                                                                                                                                                           cGTAACCGGAAAGGAGAATGGCGGGCATCAAAAAGACAAACACGCC‑GGTAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCa                                                                           <  1:1327254/185‑1 (MQ=255)
                                                                                                                                                                                                                   gAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAAAAAAGGGATGAAATCAACAGAACGAATCCCACCAGCGATAAACCATAAATCAc                                                                                                                                        >  1:1349778/1‑116 (MQ=255)
                                                                                                                                                                                                                   aaaaGGAAAATGGAGAGAACCAGAAAAACAACCACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCGCCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGc                 <  1:1159746/234‑1 (MQ=255)
                                                                                                                                                                                                                    aaaGGAGAATGGGGGGAATCATAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGAt             <  1:388449/238‑1 (MQ=255)
                                                                                                                                                                                                                    aaaGAAGAATGGCGGGCATCATAAAGACAAACACGCC‑TATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCCACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATg            <  1:1669230/239‑1 (MQ=255)
                                                                                                                                                                                                                     aaGGAGAATGGCGGGCATCATAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCa                                                                                                                                                <  2:705794/106‑1 (MQ=255)
                                                                                                                                                                                                                     aaGGAGAATGGCGGGCATCATAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCa                                                                                                                                                >  1:705794/1‑106 (MQ=255)
                                                                                                                                                                                                                     aaGGAGAATGGCGGGCATCATAAAAACAAACACGCC‑GGTAAACGCCTGCCGTTGTGTTGAACAGAGTGATGAAAACAACAGACCGAAACCCACCAGCGCTAAACCATAAATCACCCACGTAAAACAGAAAAGCGCCAGCGATCCGGCGACGGGGACTTGATACCCCCAGAAACCAAACGCCAAca                                                                 >  2:1327254/1‑185 (MQ=255)
                                                                                                                                                                                                                      aGGAGAATGGCGGGAATAAAAAAAAAAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATGaa          <  2:436529/239‑1 (MQ=255)
                                                                                                                                                                                                                        gagaATGGCGGGCATCATAAAGACAACCACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAAtg                                                            <  1:209408/187‑1 (MQ=255)
                                                                                                                                                                                                                           aaTGGCGGGCATCATAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAAAAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCaa                                                                          >  1:333740/1‑170 (MQ=255)
                                                                                                                                                                                                                               gggggAATCATAAAAACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATGAAGATCTGcc  <  1:1268994/236‑1 (MQ=255)
                                                                                                                                                                                                                               gcggaaatcatAAACACAACCACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATGAAGATCTGc   <  1:650231/231‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                    
TACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAAT‑CTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCC‑GATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATGAAGATCTGCC  >  NC_000913/825003‑825462

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑