Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,098,3910AC22.9% 99.2 / 11.8 48Y219S (TAC→TCC) rhaStranscriptional activator of rhaBAD and rhaT
Reads supporting (aligned to +/- strand):  ref base A (26/11);  new base C (0/11);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 3.12e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGT  >  NC_000913/4098160‑4098624
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                         
cGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCGGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGaa                                                                                                                                                                                                                                     >  1:81114/1‑238 (MQ=255)
                            gCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGGACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGAGACCTGAACCGCCTGGGACCGAGTAAAACCCGGCATc                                                                                                                                                                                                        >  2:112396/1‑239 (MQ=255)
                                ttACTGCTGTGGCGGCTAAAAAGCATTTGGCAGGAAACCCTGAAAAACAGCGCTTCACGTCTAAACTTGCTTCGGCCCTGGCTGGAGGACCATTTTGCCGAGGAGGTAAATTGGGATGCCGGGGCGAATCAATTTTCTCTTTCACGGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACGGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCt                                                                                                                                                                                                    <  2:31283/239‑1 (MQ=255)
                                ttACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCt                                                                                                                                                                                                    >  1:31283/1‑239 (MQ=255)
                                     gCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCt                                                                                                                                                                                                       >  1:397767/1‑231 (MQ=255)
                                                gTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGGCTGATGAAAACCCCACATTCGCTACGCCCCAGCCAGGc                                                                                                                                                                                     >  2:666534/1‑238 (MQ=255)
                                                         gTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCCCAGCGAGGCCAACGtt                                                                                                                                                                              >  1:309571/1‑236 (MQ=255)
                                                           ttGCAGGAGAACCTGGAAACCAGCGCATCCCGTCTCAACTTGCTTCTGGCCTGGCTGGAGGCCCATTTTGCCGATGAGGTAAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGa                                                                                                                                                                         <  1:112396/239‑1 (MQ=255)
                                                                   taaCTTGGAATACAGCTCATCGCGTCTCAACTTGGTTCTGGCCTGGTGGAAGGCCCATTTTGCGAATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTCCTCTTTCACTGCGTACGCTACAC‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATCCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGcc                                                                                                                                                                  <  2:464403/237‑1 (MQ=255)
                                                                    aaCCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCt                                                                                                                                                                 >  1:155523/1‑238 (MQ=255)
                                                                                aGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCCTTGGAtt                                                                                                                                                    >  1:147014/1‑239 (MQ=255)
                                                                                      tcgcGTCGCAGCTTGCTTGTGGCGTGGCTGGAGGACCATTTTGCCGATGAGGTAAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGa                                                                                                                                              <  2:81114/236‑1 (MQ=255)
                                                                                                    gCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTg                                                                                                                                                                                                                              >  2:538108/1‑145 (MQ=255)
                                                                                                    gCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTg                                                                                                                                                                                                                              <  1:538108/145‑1 (MQ=255)
                                                                                                     cTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATCCCTGAACCGCCTGCGACTGATGAAAGc                                                                                                                                                                                                                <  1:506602/158‑1 (MQ=255)
                                                                                                     cTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGc                                                                                                                                                                                                                >  2:506602/1‑158 (MQ=255)
                                                                                                                         ccATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTTACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTa                                                                                                                                                                                                   >  2:500610/1‑151 (MQ=255)
                                                                                                                         ccATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTa                                                                                                                                                                                                   <  1:500610/151‑1 (MQ=255)
                                                                                                                          cATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGAACTCGCCTATCGCTGTGGATTCAGCGCAAGTAACCACTTTTCGAAGCTTTTTCGCCGAGAGtt                                                                                                          >  2:691572/1‑239 (MQ=255)
                                                                                                                                              ttGGGATCCCGTGGCGAATCAATTCTCTGTGTCACGGCGTAC‑CTTCATCCGGCAGCTTAGGCAGCAACC‑GGGACTGACCCCTACGCGATCCCTAAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGata                                                                                      <  2:147014/239‑1 (MQ=255)
                                                                                                                                                 ggATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGa                                                                                                                                                                         >  2:33139/1‑153 (MQ=255)
                                                                                                                                                 ggATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGa                                                                                                                                                                         <  1:33139/153‑1 (MQ=255)
                                                                                                                                                 ggATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATAtt                                                                                    >  1:373109/1‑238 (MQ=255)
                                                                                                                                                      ccGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATCCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCtg                                                                                                                                                          <  2:735687/163‑1 (MQ=255)
                                                                                                                                                      ccGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCtg                                                                                                                                                          >  1:735687/1‑163 (MQ=255)
                                                                                                                                                       ctgcgCGAATCAATTTACTTTTCCCCTGCGCACGCTACCT‑CGGCCGCTTAAGCAACAAAC‑GGGACGGACGCCGCAGCGATCCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAg                                                                              <  2:155523/234‑1 (MQ=255)
                                                                                                                                                                tCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATg                                                                    >  1:478554/1‑239 (MQ=255)
                                                                                                                                                                tCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGGGGATTCAGCGACAGTAACCACTTTTCGACGCCTTTTTGCCGAGAGGTTAACTGGTTACCGCGGGGTATTTGCCAGGGACGGGATg                                                                    >  2:61247/1‑239 (MQ=255)
                                                                                                                                                                                     tACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGTACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAt                                               >  2:676136/1‑239 (MQ=255)
                                                                                                                                                                                     tACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGAATTTCGCCAGGGACGGGATGGCTTTCTGCATTAACGCGAAt                                               >  2:299246/1‑239 (MQ=255)
                                                                                                                                                                                        gCTAAAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtctt                                            <  1:691572/239‑1 (MQ=255)
                                                                                                                                                                                         cTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATCCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGc                                                                                 <  1:622865/201‑1 (MQ=255)
                                                                                                                                                                                         cTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGc                                                                                 >  2:622865/1‑201 (MQ=255)
                                                                                                                                                                                          tACAT‑CGGCAGCTTAGGCCGCAAGC‑GGGACTGACGCCTCAGCGATCCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtcttct                                          <  1:666534/239‑1 (MQ=255)
                                                                                                                                                                                                 ggCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCa                                        <  1:483993/235‑1 (MQ=255)
                                                                                                                                                                                                     gttaaaGCAGAAAGG‑GGAACTGCCGCTTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTAttt                                <  2:618320/235‑1 (MQ=255)
                                                                                                                                                                                                     atttAAACCGCAAACGGGGCCTGCCGCCTGGGGGATCCCGGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTAt                                  <  2:373109/236‑1 (MQ=255)
                                                                                                                                                                                                                aaaC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCga                                                                                                               <  1:151964/149‑1 (MQ=255)
                                                                                                                                                                                                                aaaC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCga                                                                                                               >  2:151964/1‑149 (MQ=255)
                                                                                                                                                                                                                        aCTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGc                                                                                                                                                                                 <  2:591433/76‑1 (MQ=255)
                                                                                                                                                                                                                        aCTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGc                                                                                                                                                                                 >  1:591433/1‑76 (MQ=255)
                                                                                                                                                                                                                        aCTGACGCCTAAGCGATCCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGc                  <  1:299246/235‑1 (MQ=255)
                                                                                                                                                                                                                            aCGCCTCAGCGATCCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGa                                                                                        <  2:718502/161‑1 (MQ=255)
                                                                                                                                                                                                                            aCGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGa                                                                                        >  1:718502/1‑161 (MQ=255)
                                                                                                                                                                                                                            aCGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGTGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATc          >  2:439946/1‑239 (MQ=255)
                                                                                                                                                                                                                             cGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCa         >  2:16327/1‑239 (MQ=255)
                                                                                                                                                                                                                             cGCCTAAGCAATCCCTGAGCCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGa                <  1:61247/232‑1 (MQ=255)
                                                                                                                                                                                                                                        tACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGt  >  1:414978/1‑235 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                         
CGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCTACAT‑CGGCAGCTTAAGCAGCAAAC‑GGGACTGACGCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGT  >  NC_000913/4098160‑4098624

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