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.

TACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATGAAGATCTGCC  >  NC_000913/825003‑825462
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                    
tACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGGATCCTAAAGAc                                                                                                                                                                                                                                >  1:75408/1‑239 (MQ=255)
 aCATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGGGAATGGCGGGGACCAAAAAAACa                                                                                                                                                                                                                               >  1:593985/1‑239 (MQ=255)
       ccTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCCAAAAGACAAaaac                                                                                                                                                                                                                          >  1:1589613/1‑238 (MQ=255)
                 aaaCATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCAAAAAGACAAACACCCCGCTAAAc                                                                                                                                                                                                               >  1:630218/1‑238 (MQ=255)
                       cgcgTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGCCTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAAAACACCGGTAAAAGCct                                                                                                                                                                                                           >  1:111865/1‑236 (MQ=255)
                                     cctgtGTCGTGGCCGTAACCACCAGTGGCGACCACAAACTATACCCAAAATTATCCGTACCGGCATCATTAAATAAATGCTGCGTGGTAATGTCCGAAAAGGGACAAATAGGGTTAAACCAGGTCGGTTTTGGCAGCCGTACCGGCATTTTTCCCACCGGAAAACCGTAACAGAAAAGGAGAAGGGGGGGCATAATAAAGACAACCACGCCGATAAACGCCTGCTGTTGTGTTGAACa                                                                                                                                                                                           <  2:111865/234‑1 (MQ=255)
                                               ggCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCATGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCCTAAAGACAAACACCCCGAAAAAACCCTTCTTTTGGGTTGAAAAGAGggagaaaa                                                                                                                                                                                >  1:2089221/1‑234 (MQ=255)
                                               ggCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCACCATAAAGACAAACAAGCCCCAAAACGCCTGATGTTGTGTTGTACAAAGGGGTGaaa                                                                                                                                                                                >  1:1444758/1‑239 (MQ=255)
                                                gCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAAAACGCCGATAAAAGCCCGCTGTTGTGTTGAACAGAGGGATGaaa                                                                                                                                                                                >  1:499899/1‑238 (MQ=255)
                                                    ttATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAAAACGCCGATAAACGCCTGCTGTTGTGTTGAACACAGTGATGAAaacaac                                                                                                                                                                           >  1:1576129/1‑239 (MQ=255)
                                                              tAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCCGCTGTTGTGTTGAACAAAGGGATGAAATCAACAAACCGAAt                                                                                                                                                                  >  1:329016/1‑238 (MQ=255)
                                                               agcgtgaacaAACGACTCCAAACAATACCCGTACTCGTATCCTTAAAATAGATCTGCTGGGTAATGACCTTAAAGTGGCGAAGAGGGTTATTCCACTGCATGTTTAGCATCCGTCGCGGCATGTTTCCCACCGGCGAAACGGAACTGAAAAGAAGATGGGCGGTAATGAGAAAAAAAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGCCCGaa                                                                                                                                                                   <  2:1505664/229‑1 (MQ=255)
                                                                              ttCCACAAAATATCCAGACTCGCATCCTTCAAATAAATCTGTTTGGAAATTTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTACCCGGAAAGGAGAATGGCGGGCATCAAAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGaa                                                                                                                                                                                 <  1:842079/207‑1 (MQ=255)
                                                                                       ttatcagatcTCGCATCTTGAAAATATATCGGCTTGGTAATGTCCGAAAGGGGGCGAATAGGGTAAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTACCCCCGGAGAAACGTAACGGAAAAGGAGAATGGCGGGCAAAAAAAAGAAAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCAc                                                                                                                                        <  2:630218/230‑1 (MQ=255)
                                                                                       ataccaaGACTCTCATCCTTCAAATAAATCTGCTTGTTAAGCGCCGCAAAGTGTCGAATATGGTTACTCCACGGCAGGTGGTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCAc                                                                                                                                        <  2:1139591/233‑1 (MQ=255)
                                                                                          tCCAGCCTCGCATCCTTAAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCCCCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGAAAAACACGCCGATAAACGCctgctg                                                                                                                                                                                                       <  1:749013/173‑1 (MQ=255)
                                                                                                    cATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGGATCATAAAGGCAAAAACCCCGGATAACCCCCGGCGTTGTGGTGaa                                                                                                                                                                                             >  2:749013/1‑173 (MQ=255)
                                                                                                             aaaaaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCACCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGAAAAGGAGAAGGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCtgttgtg                                                                                                                                                                                                  <  1:1640131/156‑1 (MQ=255)
                                                                                                                        cTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGGATCATAAAGACAAAAACGCCGAGAAAACGCTGCTGTTGTGTTGAACagag                                                                                                                                                                                        >  2:1640131/1‑159 (MQ=255)
                                                                                                                            gTAATGTCCGTAAGGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATTTTTTCCACCGGAAAAACGTAACCGGAAAGGAGAATGGCGGGCTTCAAAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGaa                                                                                                    <  1:314979/239‑1 (MQ=255)
                                                                                                                                  tCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGAAAAAAACGCCGATAAACCCCTGCtgttgtgtt                                                                                                                                                                                                >  1:176404/1‑141 (MQ=255)
                                                                                                                                  tCCGTAAAGGGGCAAATAGGGTTAACCCACGTCAGTTTTTGCAGCCATCCCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAAGGGGGGGCATCATAAAGACAAACACGCCGATAAAc                                                                                                                                                                                                               <  2:148304/126‑1 (MQ=255)
                                                                                                                                      tAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAAAAGCGCCAGCGATCCGGCGTAGGGTATTTGa                                                                                          >  1:1498018/1‑239 (MQ=255)
                                                                                                                                            ggCGAATGGGTTTAATCCACGTCAGTTTTTGCAGCCATCCCGGCATGTTTTCCCCCGGAAAAACGTAACCGGAAAGGAGAATGGCGGGCATCAAAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATa                                                                                        <  1:1228598/235‑1 (MQ=255)
                                                                                                                                            ggCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCtgtt                                                                                                                                                                                                     >  1:148304/1‑126 (MQ=255)
                                                                                                                                             gggAATAGGGTGAAACCCGTTGAGGTTTTGGAGCAATACGGGCATGTTTACCACCGGAGAAACGTAACCGAAAAGGGGAAGGGGGGGAATAAAAAAGAAAAACACGCCGATAAACGCctg                                                                                                                                                                                                          <  2:1245331/118‑1 (MQ=255)
                                                                                                                                              cGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGGATCATAAAAACAAAAACCCAGCAAAACACCCGTCGTGTGGTTTAACACAGGGATGAAAAAAAACAACAGAAACCCACCACCCATAAACCAAAAAACAACAACCTAAACTAGAAAACAGCCAGCCATCCGGCGAAGGGGACTTTATACgaaca                                                                                  >  2:977816/1‑235 (MQ=255)
                                                                                                                                                      gtaaacacACGTAATGTTGTCCAGCCATACCGGCATGTTTTCCACCGGAGAACCGTAACCGGAAAGGAGAATGGGGGACATCAAAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATGCGCCCCGATACCa                                                                           <  2:1589613/231‑1 (MQ=255)
                                                                                                                                                       ttAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGAc                                                                                                                                                                       >  1:1428024/1‑145 (MQ=255)
                                                                                                                                                       ttAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCCATAAACGCCCTCtgttgggtt                                                                                                                                                                                                >  1:1245331/1‑120 (MQ=255)
                                                                                                                                                          aTCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCg                                                                       >  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)
                                                                                                                                                               cGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGaa                                                                                                                                                                                             >  2:895367/1‑115 (MQ=255)
                                                                                                                                                               cGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGaa                                                                                                                                                                                             <  1:895367/115‑1 (MQ=255)
                                                                                                                                                               cGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGGGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGc                                                                  >  1:1271939/1‑238 (MQ=255)
                                                                                                                                                               cGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGGGATGAAATCAACAGAACGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATAACAATCGCCAGca                                                                 >  1:1967350/1‑239 (MQ=255)
                                                                                                                                                               cGTCAGGTTTGGCAGCCTACCCGGCATGTTTCCCACGGGAGAAACGTAACGGGAAAGGAGAATGGCGGGAATCAAAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGca                                                                 <  1:284135/239‑1 (MQ=255)
                                                                                                                                                                      ttttGCAGCAATCCCGGATTGTTTCCCACCGGAAAAACGTAACCGAAAAGGAAAATGGCGGGAATCAAAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAAtg                                                            <  1:834980/237‑1 (MQ=255)
                                                                                                                                                                      ttttGCAACCATACCGGAGTGTTTTCCACCGGAAAAGCGTAGCCGGAAAGGAGAATGGCGGGCATCATAAAAAAAAACACGCCGATAAACGCCTGCGGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAg                                                                                                                      <  2:1011943/179‑1 (MQ=255)
                                                                                                                                                                              ccATACCGGCATTTTTTCCACCGGAGAAACGTACCCGGAAAGGAGAAGGGGGGGCATCAAAAAGACAAACACGCCGATAAACGCCTGCtgttgt                                                                                                                                                                                                   <  1:1423391/94‑1 (MQ=255)
                                                                                                                                                                                aTACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAAACAAACCCCCCGATAAACGCCTGCTGTTTTGTTGAACAGAGGGATGAAATCAACAGACCGAATCCCACCACCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGAt                                                                                                            >  1:1011943/1‑179 (MQ=255)
                                                                                                                                                                                        aTGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAAACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACaga                                                                                                                                                                                         >  2:1423391/1‑94 (MQ=255)
                                                                                                                                                                                           ttttCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCACCATAAAGAAAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGcc                                                                                    <  1:475241/192‑1 (MQ=255)
                                                                                                                                                                                           ttttCCACCGGAAAACCTTACCCGAAAAGGAGAATGGCGGGCATCAAAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATa                                                                                        <  1:2064171/188‑1 (MQ=255)
                                                                                                                                                                                                   cgaggaAACGTACCCGGAAAGTAGAATGGCGGGCATCAAAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGAt                                                                                              <  1:1284515/170‑1 (MQ=255)
                                                                                                                                                                                                   cGGAGAAACGAACCGGGAAGGGAAAAGGGGGGGAAAAAAAAAAACAACCACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGc                               <  1:752148/237‑1 (MQ=255)
                                                                                                                                                                                                     gagaAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAAAACGCCGATAAACGCCTGCTGTTGTGTTTAACAGAGTGATGAAATCAACAGACCGAATCCCCCCAGCGATAAAACATAAAATACCCACGGAAAGGAGAACAACGCCCGGGGTCCGGGGGAGGGGATTTGATACGCCCAGATa                                                                              >  2:2064171/1‑188 (MQ=255)
                                                                                                                                                                                                         aaCGAAACGGAAAGGGAGAAGGGCGGTCATCAAAAAGAAAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAAc                                                                                                                                                  <  2:1349778/116‑1 (MQ=255)
                                                                                                                                                                                                           cGTAACCGGAAAGGAGAATGGCGGGCATCAAAAAGACAAACACGCCGGTAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCa                                                                           <  1:1327254/185‑1 (MQ=255)
                                                                                                                                                                                                                   gAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAAAAAAGGGATGAAATCAACAGAACGAATCCCACCAGCGATAAACCATAAATCAc                                                                                                                                        >  1:1349778/1‑116 (MQ=255)
                                                                                                                                                                                                                   aaaaGGAAAATGGAGAGAACCAGAAAAACAACCACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCGCCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGc                 <  1:1159746/234‑1 (MQ=255)
                                                                                                                                                                                                                    aaaGGAGAATGGGGGGAATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGAt             <  1:388449/238‑1 (MQ=255)
                                                                                                                                                                                                                    aaaGAAGAATGGCGGGCATCATAAAGACAAACACGCCTATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCCACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATg            <  1:1669230/239‑1 (MQ=255)
                                                                                                                                                                                                                     aaGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCa                                                                                                                                                <  2:705794/106‑1 (MQ=255)
                                                                                                                                                                                                                     aaGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCa                                                                                                                                                >  1:705794/1‑106 (MQ=255)
                                                                                                                                                                                                                     aaGGAGAATGGCGGGCATCATAAAAACAAACACGCCGGTAAACGCCTGCCGTTGTGTTGAACAGAGTGATGAAAACAACAGACCGAAACCCACCAGCGCTAAACCATAAATCACCCACGTAAAACAGAAAAGCGCCAGCGATCCGGCGACGGGGACTTGATACCCCCAGAAACCAAACGCCAAca                                                                 >  2:1327254/1‑185 (MQ=255)
                                                                                                                                                                                                                      aGGAGAATGGCGGGAATAAAAAAAAAAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATGaa          <  2:436529/239‑1 (MQ=255)
                                                                                                                                                                                                                        gagaATGGCGGGCATCATAAAGACAACCACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAAtg                                                            <  1:209408/187‑1 (MQ=255)
                                                                                                                                                                                                                           aaTGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAAAAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCaa                                                                          >  1:333740/1‑170 (MQ=255)
                                                                                                                                                                                                                               gggggAATCATAAAAACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATGAAGATCTGcc  <  1:1268994/236‑1 (MQ=255)
                                                                                                                                                                                                                               gcggaaatcatAAACACAACCACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATGAAGATCTGc   <  1:650231/231‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                    
TACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCCGATGAAGATCTGCC  >  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: