Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913993,9170AC45.9% 98.0 / 15.0 61V43G (GTG→GGG) ssuBaliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (29/4);  new base C (0/28);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 2.55e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATGAT  >  NC_000913/993770‑994066
                                                                                                                                                       |                                                                                                                                                     
aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCccacca                                                                                                                                                    <  2:520974/151‑1 (MQ=255)
aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCccacc                                                                                                                                                     >  2:521047/1‑150 (MQ=255)
 gCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCcaccac                                                                                                                                                   >  1:1166356/1‑151 (MQ=255)
 gCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCcaccac                                                                                                                                                   >  1:733348/1‑151 (MQ=255)
  cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCaccacc                                                                                                                                                  >  2:1083132/1‑151 (MQ=255)
   aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCg                                                                                                                                                 >  2:351222/1‑151 (MQ=255)
   aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCg                                                                                                                                                 >  1:147078/1‑151 (MQ=255)
          cATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACaa                                                                                                                                           >  2:622476/1‑150 (MQ=255)
          cATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACaa                                                                                                                                           >  1:738616/1‑150 (MQ=255)
           aTCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAAc                                                                                                                                         >  1:175374/1‑151 (MQ=255)
            tCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACt                                                                                                                                        >  1:994660/1‑151 (MQ=255)
             ccTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTg                                                                                                                                       >  2:771654/1‑151 (MQ=255)
                  aaCATCTTCGGCGTGTCTTCCTAATTTC‑CAGCCGGGGGTGGGGTCCCCGCAAACACACCGCGTGCGGTTGG‑GTTTCCAAGCCCCCCCAGAGGGCGCAGCAGGG‑TTTTTTTC‑CCCCCCCCGCGGCGGCCCCCCACCGCCACAAACTGACCTg                                                                                                                                  <  2:908230/151‑1 (MQ=255)
                   aCATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGc                                                                                                                                 >  2:1237458/1‑151 (MQ=255)
                     atcatTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCg                                                                                                                               >  1:1051501/1‑151 (MQ=255)
                      tcatTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                              >  1:259834/1‑151 (MQ=255)
                           cgcgTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAatat                                                                                                                         >  2:302454/1‑151 (MQ=255)
                           cgcgTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAatat                                                                                                                         >  2:818079/1‑151 (MQ=255)
                                    tCCTAATTTC‑CAGCCGGGGGTGTGGTCCCCGCAAACACACCGCCGGGGTTTGGCTTTCCC‑AGCCCTCCCAGCAGGCGCAGCGGGT‑TTTTTTTC‑CCCCCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                                   <  2:290863/148‑1 (MQ=255)
                                      cttaaTTTCCAGCCGGGGGTGGGGGGCCCGCTACCCCACCGCCGGGGGTGGGCGTTCCC‑AGCCCGGCCAGCGGGCGCAGCGGGG‑TTTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                                 <  2:238644/146‑1 (MQ=255)
                                       tGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                             >  1:456551/1‑151 (MQ=255)
                                         atttt‑CAGCGGGGGGTGGGGTCCCCGCAACCACTCCGCCGGGGTTGGGCTTTCCC‑ACCCCGCCCAGCGGCCGCAGCAGGT‑TTTTTTCC‑CCCCCCCCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                                  <  1:286388/142‑1 (MQ=255)
                                          aTTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                          >  2:761714/1‑151 (MQ=255)
                                                agccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                                     >  1:677316/1‑151 (MQ=255)
                                                  ccgggggTGGGGGGCCCCCAACCACACCGCCGGGGTTGGGCTTTCCC‑ACCCCTGCCAGCGGGCGCAGAGGGT‑TACTTTCG‑CCACCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGa                                                                                                    <  1:286410/145‑1 (MQ=255)
                                                  ccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                                   >  1:421602/1‑151 (MQ=255)
                                                   cggggtgggggtcCCCGCAAACCCCGCGCCTGGGTTTGGCTTTCCC‑CGCCCCCCCAGCCGGCGC‑GCAGGGTTTTTTTTC‑CCACCCCCGCGGCGCCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGtt                                                                                                         <  1:591343/135‑1 (MQ=255)
                                                       ggTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                              >  2:496071/1‑151 (MQ=255)
                                                           tggtgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTc                                                                                          >  1:511368/1‑151 (MQ=255)
                                                              tgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgc                                                                                       >  2:1188700/1‑151 (MQ=255)
                                                                cccGCAAACACCCCGCCTGCGTTGGGCTTTCCC‑ACCCCCCCCAGCAGGCGCAGCAGGT‑TATTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                                     <  2:675016/151‑1 (MQ=255)
                                                                      aaCACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                               >  2:792669/1‑151 (MQ=255)
                                                                          cctcccccggggTGGGCTTTCCC‑AGCCCTCCCAGCGGCCGCAGGGGGT‑TTTTTTCC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                            <  1:464466/140‑1 (MQ=255)
                                                                           aTCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                          >  2:146518/1‑151 (MQ=255)
                                                                              gcggggggTGGGCGTTTCC‑AGCCCGCCCAGCAGGCGCAGCGGGG‑TATTTTCC‑CCCCCCCCGCGGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                            <  1:388676/140‑1 (MQ=255)
                                                                                  ggggTTGGCTTTCCC‑AGCCCTCCCAGCGGGCGCAGCGGGG‑TTTTTTCC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                                     <  1:191613/147‑1 (MQ=255)
                                                                                  ggggTGGGCTTTCCC‑GGCCCGGCCAGCAGGCGCAGCGGGG‑TATTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCt                                                                           <  2:59390/141‑1 (MQ=255)
                                                                                      ttGGCTTTCCC‑ACCCCTGCCAGCGGGCGCAGCGGGT‑TCCTTTTC‑CCCCCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAtt                                                                  <  2:836076/148‑1 (MQ=255)
                                                                                       tgggctttcccaccccTGCCGGCGGCCCCAGCAGGT‑TCTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                          <  2:940039/126‑1 (MQ=255)
                                                                                           gTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                           >  1:775475/1‑150 (MQ=255)
                                                                                            ttccc‑AGCCCTCCCAGCAGCCGCAGCGGGT‑TTTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaaca                                                          <  2:237704/147‑1 (MQ=255)
                                                                                             ttcccaccccTCCCAGCGGGCGCGGCAGGT‑TTTTTTCG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaa                                                            <  1:260335/140‑1 (MQ=255)
                                                                                             ttccaaccccggCCGGCGGGCGCGGCGGGT‑TTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGt                                                                                <  1:216164/117‑1 (MQ=255)
                                                                                              ccc‑AGCCCGGCCAGCAGGCGCAGCAGGT‑TTTTTTCC‑CCCCCCCCGCGGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGca                                                             <  1:695277/144‑1 (MQ=255)
                                                                                                     ccgcccAGCGGGCGCGGCAGGT‑TCCTTTGC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGc                                                     <  2:928477/143‑1 (MQ=255)
                                                                                                     cccgCCAGCGGGCGCAGGGGGT‑TTTTTTTC‑CCCCCCCCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGc                                                     <  1:796105/144‑1 (MQ=255)
                                                                                                     cccgCCAACAGCCGCAGCGGGT‑TCTTTTTC‑CCCCCCCCCCTCCCCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACAATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                                    <  2:81751/145‑1 (MQ=255)
                                                                                                     ccTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGc                                                 >  2:1163372/1‑151 (MQ=255)
                                                                                                              agGCGCAGCACGG‑TACTTTAC‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATACGTAACTCCAGTTGGTTCAGGCCGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGGTc                                         >  1:205235/1‑150 (MQ=255)
                                                                                                               ggCGCAGCGGGG‑TTTTTTTC‑CCACCACCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                              <  2:1155677/144‑1 (MQ=255)
                                                                                                               ggCGCAGCAGGT‑TCTTTTCC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                         <  2:835235/149‑1 (MQ=255)
                                                                                                                   cagcagGT‑TATTTTCC‑CCCCCCCCCCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTACTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGa                                   <  1:176953/151‑1 (MQ=255)
                                                                                                                      gcggggtctttttC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAg                                  <  1:94644/141‑1 (MQ=255)
                                                                                                                        ggg‑TTCTTTTGCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAg                               <  2:798632/151‑1 (MQ=255)
                                                                                                                                 ttcccccccaccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                           <  1:877440/140‑1 (MQ=255)
                                                                                                                                  tG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc                     >  2:1116851/1‑151 (MQ=255)
                                                                                                                                     ccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt                   <  1:162706/145‑1 (MQ=255)
                                                                                                                                       ccccccgcggcccccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg                 <  1:1067912/138‑1 (MQ=255)
                                                                                                                                        cccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg                 <  1:587341/147‑1 (MQ=255)
                                                                                                                                        cccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGc                <  1:1231394/148‑1 (MQ=255)
                                                                                                                                         ctccgcggggccccACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt                   <  2:88576/136‑1 (MQ=255)
                                                                                                                                          accGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct              >  1:71513/1‑151 (MQ=255)
                                                                                                                                          accGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct              >  2:857703/1‑151 (MQ=255)
                                                                                                                                           ccccTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct              <  2:1247375/147‑1 (MQ=255)
                                                                                                                                           ccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctc             >  1:732839/1‑151 (MQ=255)
                                                                                                                                              cggcgcCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct           <  1:1154878/148‑1 (MQ=255)
                                                                                                                                                  ggCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaa       >  2:387134/1‑150 (MQ=255)
                                                                                                                                                       accaccGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTCCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATGAt  >  2:79085/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGG‑TACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATGAT  >  NC_000913/993770‑994066

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