Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,483,2540GA29.5% 87.1 / 10.5 44S29L (TCA→TTA) emrKtripartite efflux pump membrane fusion protein EmrK
Reads supporting (aligned to +/- strand):  ref base G (23/8);  new base A (0/13);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 3.92e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTC‑GAACGTAGTTCGTATCTTTATGATTAACGACAGTGAC‑ACTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTT  >  NC_000913/2483006‑2483503
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         
ataGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTC‑GAACGTAGTTCGTATCTTTATGATTAACGACAGTGAC‑ACTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGa                                                                                                                                                                                                                                                          >  1:227279/1‑250 (MQ=255)
     ttATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTC‑GAACGTAGTTCGTATCTTTATGATTAACGACAGTGAC‑ACTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGc                                                                                                                                                                                                                                                    >  1:265114/1‑251 (MQ=255)
          cGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTC‑GAACGTAGTTCGTATCTTTATGATTAACGACAGTGAC‑ACTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGGGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCaata                                                                                                                                                                                                                                                >  2:743974/1‑250 (MQ=255)
           gTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTC‑GAACGTAGTTCGTATCTTTATGATTAACGACAGTGAC‑ACTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCaataa                                                                                                                                                                                                                                               >  1:243899/1‑250 (MQ=255)
             ttGCCGAACAATATGTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTC‑GAACGTAGTTCGTATCTTTATGATTAACGACAGTGAC‑ACTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCaataaat                                                                                                                                                                                                                                             <  2:557567/250‑1 (MQ=255)
                      aaTATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTC‑GAACGTAGTTCGTATCTTTATGATTAACGACAGTGAC‑ACTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCg                                                                                                                                                                                                                                   <  1:622775/251‑1 (MQ=255)
                                aGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTC‑GAACGTAGTTCGTATCTTTATGATTAACGACAGTGAC‑ACTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATaa                                                                                                                                                                                                                            >  1:794047/1‑248 (MQ=255)
                                     attTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTC‑GAACGTAGTTCGTATCTTTATGATTAACGACAGTGAC‑ACTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGt                                                                                                                                                                                                                    >  2:335557/1‑251 (MQ=255)
                                        tttAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTC‑GAACGTAGTTCGTATCTTTATGATTAACGACAGTGAC‑ACTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTAtt                                                                                                                                                                                                                 >  1:235049/1‑251 (MQ=255)
                                           agttttATTGAGGGCGATGGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCCCCTTGTC‑GAACGTAGTTCGTATTTTTATGATTAACGACAGTGAC‑ACTACCTGAGACTT‑GGGCAGAAATTGGATCTGCATTCCCCGTGACATGGGCGTCATCTGTACTAATCATGTCTTCTAAT‑CCCATTGACCAATAGGCAT‑AGCCCCCTAAAAACGAAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTttc                                                                                                                                                                                                              <  1:433329/248‑1 (MQ=255)
                                                      gTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTC‑GAACGTAGTTCGTATCTTTATGATTAACGACAGTGAC‑ACTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGa                                                                                                                                                                                                    >  1:491112/1‑250 (MQ=255)
                                                          gatggtgtcatcgtttttATCCATTAAAACTAAAAGCTC‑CCCTGTCCGAACGTATTTCTTTTTTTTATGATTATCGCCAGTAAC‑ACTCCCGGAGATTT‑GTGCAAAAATTTAATCTGCATCCCCCTTAACATGGCCGTCATCTTTACTATTCATGTCTTC‑AATACCCATTGCCCAATAACAAA‑AGACCCTTAAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGtt                                                                                                                                                                                               <  2:265114/237‑1 (MQ=255)
                                                             atggactgctgTTTTACCCTGTGAAACAAAAATGTCCCTTTGTC‑GAACGTAGTTCTTTTCTTTGTGATTAACTACATGGCC‑CCTACATGAGATTT‑GTGCAGAAATTGTATCTTCATTCCCCTTGAAATAGGCTCCTTCTGTACTAATCATGCCTTCTAAT‑TCTATTAACAAATAGGAAT‑AGGCACCTAAAACCGAAAAAAAAAATCCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGttttt                                                                                                                                                                                            <  2:233682/241‑1 (MQ=255)
                                                                            tcccaaggAAACTAAAATGTCACCTTGTC‑GAACTTAGTTCTTCTCTTTTTGATTACCGACAGTGAC‑ACCCCCGGACACTT‑GTGCATAAATTGGATCGCCATTCCCAGTGACACAGTCGTCCTCTGTCCTAAGCTTGTCTTCTAAT‑TCCATTGCCCAATAG‑CATTAGCCCCCTGAAAACGCAAAAAAAAAAAGTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAGTTAATc                                                                                                                                                                              <  1:793212/243‑1 (MQ=255)
                                                                             cccAGTGAAATTAAATTGTCCCCTTGTC‑GAACGTATTTCTTTTCTTTTTAATTAACGACAGTGAC‑CCTTCCTGA‑ACTTTGTGCAAAAATTGGATCTGCATTCCCCGTGACATAGGCGTCTTCTGTACTAATCATGTCTTCTAAT‑CCCATTGACCAATGGGCAA‑AGGCCCCGAAAAACCAAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTg                                                                                                                                                                            <  1:564825/250‑1 (MQ=255)
                                                                                 gtaaaaCTAAAATGTCACCTTGTC‑AAACGTAGTTGGTATTTTTATGATTAACGACAGGGCC‑CCTCCCGGAGATTT‑TTGAAGAATTTGGTTCTCCATTCCCCGTGACAGAGGCTTCATCTTTACTAATCATGTCTTCAAAT‑TCAATTGACCAATAGGCAT‑AGCCACCTAAAAACTCAATAAATAAACTTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTcc                                                                                                                                                                        <  2:568027/248‑1 (MQ=255)
                                                                                        tAAAATGTCACTTTGTC‑GTACGTTGTTCGTTTCTTTAGGATTAACGACAGTGAC‑TCTACCTAAGACTT‑GTGAAAATATTTGATCTGCATTCCCCGTTACATAGGCTTATTCTGTACTAATCATGTCTGCTAAT‑TCAATTGACCAAAAGGCAA‑AGCCACCTGAAAACGAAATAAATAATACTCCCGCCAATAAATAAAAGTATTTTCGTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACtatt                                                                                                                                                                  <  1:260401/250‑1 (MQ=255)
                                                                                          gaccctcttccgatgtaaacttAGTTCTTATTTTTATGATTAACGCCAG‑GACCACTCCCTGGGATTT‑TGGCAAAAATTGGATTTCCTTTCCCGGTAACATAGGCTTCATTTGTACAAATCATGTTTTCAAAT‑TCCATTGACAAATAGCAAT‑AGGCACCTAAAACCGCAAAAAATAAATCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACtattat                                                                                                                                                                <  2:804358/230‑1 (MQ=255)
                                                                                                cACCTTGTC‑GAACGTAGTTCGTATCTTTATGATTAACGACAGTGAC‑ACTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCAt                                                                                                                                                         >  1:560931/1‑251 (MQ=255)
                                                                                                   cTTGTC‑GAACGTAGTTCGTATCTTTATGATTAACGACAGTGAC‑ACTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTtc                                                                                                                                                      >  1:748852/1‑251 (MQ=255)
                                                                                                          gAACGTAGTTCGTATCTTTATGATTAACGACAGTGAC‑ACTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACtattat                                                                                                                                                                >  1:804358/1‑235 (MQ=255)
                                                                                                                                aTTAACGACAATGAC‑ACACCCTGAGATTT‑GTGGAAAAGTTTGATCTGCATTCCCCGGAACGTAGGCGTCATCTGTACTAATCGTTTCTTCTAAT‑CCCATAGCCCAATAGCCGA‑AGGCCCCTAAAACCGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATtgaatga                                                                                                                          <  1:298418/251‑1 (MQ=255)
                                                                                                                                        cAGTGAC‑ACTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGc                                                                                                                   >  1:791747/1‑250 (MQ=255)
                                                                                                                                         aGTGAC‑ACTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGaa                                                                                                                                                                                    >  1:272310/1‑184 (MQ=255)
                                                                                                                                         aGTGAC‑ACTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGcat                                                                                                                 >  1:65138/1‑251 (MQ=255)
                                                                                                                                              acaatACCTGAGACTT‑GTGAAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACAAATCATGTCTTCAAAT‑TCCATTGACCAAAAGGCAT‑AGGCACCTAAAAACGAAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGcattcat                                                                                                             <  2:245184/247‑1 (MQ=255)
                                                                                                                                                acTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCattcattc                                                                                                           >  1:207739/1‑251 (MQ=255)
                                                                                                                                                   cccTGAGACTG‑GTGAAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACAAATCATGTCTTAAAAT‑CCCATTGACAAATAGGCAT‑AGGCCCCTGAAAACGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCAttctt                                                                                                         <  1:289354/249‑1 (MQ=255)
                                                                                                                                                     cagggtcgtgttgCAGAATTTGGACCTGCATTCCCCGGTACATAGGCGTCATTTGTACGAATCTTGTCTTCTAAT‑GCCATTGACCAATAGACAT‑GGGCCCCTAAAAACACAATAAATAAATTTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCa                                                                                                       <  2:587897/240‑1 (MQ=255)
                                                                                                                                                                   aGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTg                                                                                         >  1:721106/1‑251 (MQ=255)
                                                                                                                                                                                gCATCCCCCGTGCCATAGGCGTAATCTTTCTTAATAAGGTCTTCTATT‑CCCATTGAACAATAGGCAG‑AGCCACCTGAAAACACAAAAAATAAATCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATGCATTCTTCATATCAAACACAGTTCATGTGTTGcac                                                                             <  2:298846/250‑1 (MQ=255)
                                                                                                                                                                                      cccatgaaaTAGGCGTCATCTTTACAAATAATGTCTTAAAAT‑CCCATTGACCAATAGGCAA‑AGCCACCTAAAACCGAAAAAAAAAAAACTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATcc                                                                      <  2:235049/243‑1 (MQ=255)
                                                                                                                                                                                       ccGTGACATAGGCGTCATCTGTACTAATCATTTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTAAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATccc                                                                     <  1:437624/251‑1 (MQ=255)
                                                                                                                                                                                         gTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATccc                                                                     >  1:506771/1‑249 (MQ=255)
                                                                                                                                                                                                           gTACTAATCATGTCTCCTAAT‑TCCATTGCCCAATAGGCAT‑AGGCACCTGAAAACGCAAAAAAAAAAACTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTa                                                   <  1:85700/249‑1 (MQ=255)
                                                                                                                                                                                                              cTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAAttt                                               >  1:763997/1‑250 (MQ=255)
                                                                                                                                                                                                                attCATGTTTTCAAAT‑CCCATTGCCCAATAGGCTT‑AGCCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCAGTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTaa                                             <  1:743974/248‑1 (MQ=255)
                                                                                                                                                                                                                aaTCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTg                                                                                                                                                                                                         >  1:453097/1‑94 (MQ=255)
                                                                                                                                                                                                                   cATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATg                                          >  1:224770/1‑250 (MQ=255)
                                                                                                                                                                                                                              aT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCaaaaa                               >  1:839219/1‑250 (MQ=255)
                                                                                                                                                                                                                                  ccATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAaac                           >  2:842711/1‑251 (MQ=255)
                                                                                                                                                                                                                                   acatgtaccaataggattgagcaCCTAAATACGAAATAAATAAAATTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAAGACAGTGCATGTGTTGCACATTATACCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAACAAaac                           <  2:565230/231‑1 (MQ=255)
                                                                                                                                                                                                                                           aaaaggcat‑agccaACTAAAAACGAAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGa                  <  2:603404/248‑1 (MQ=255)
                                                                                                                                                                                                                                                  at‑aggcaCCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAAcc           >  1:289964/1‑251 (MQ=255)
                                                                                                                                                                                                                                                             gAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTAATATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCACAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGCCTAAACCGCGGCtttt  >  2:478356/1‑250 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         
ATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTC‑GAACGTAGTTCGTATCTTTATGATTAACGACAGTGAC‑ACTACCTGAGACTT‑GTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAAT‑TCCATTGACCAATAGGCAT‑AGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTT  >  NC_000913/2483006‑2483503

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