Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913993,9130GC45.8% 103.6 / 26.1 72G44G (GGC→GGG) ssuBaliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (29/10);  new base C (0/33);  total (29/43)
Fisher's exact test for biased strand distribution p-value = 5.63e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCA  >  NC_000913/993768‑994060
                                                                                                                                                   |                                                                                                                                                   
gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCccac                                                                                                                                                <  2:848230/151‑1 (MQ=255)
gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCccac                                                                                                                                                >  1:590718/1‑151 (MQ=255)
  aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCccacc                                                                                                                                               <  1:4818/150‑1 (MQ=255)
    cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCaccacc                                                                                                                                            <  2:907553/151‑1 (MQ=255)
    cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCaccacc                                                                                                                                            <  1:832305/151‑1 (MQ=255)
                  gAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGAcc                                                                                                                              >  2:208384/1‑151 (MQ=255)
                    accatcatCCGCTTTTCTTCCTAATTTCCACCCAGCGGTGGGGTCCCCGTAAACACACCCCCTGGGTTTGGCTTTT‑CCAGCCCGGCCAGCAGGCGCAGCAGGGTATTTTC‑CCCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTg                                                                                                                            <  1:743092/149‑1 (MQ=255)
                    aaCATCTTCCGCGTTTCTTCCGAAATTTCAGCCGGGGGGGGGGTCCCCGCAACCACACCG‑CGGGGGTTGGCCTTTCCCAGCCCGGCCAGCAGGCGCAGAGGGGTTCTTTT‑GCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTg                                                                                                                            <  2:797880/151‑1 (MQ=255)
                      catcatCCGCGTGTCTTCCTGAATTTCAGCCGGCGGTGTGGGGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGCCCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGc                                                                                                                           <  1:716670/150‑1 (MQ=255)
                       atcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCg                                                                                                                         >  1:779366/1‑151 (MQ=255)
                        tcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGTTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                        >  2:417066/1‑151 (MQ=255)
                          atTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGa                                                                                                                       >  1:280744/1‑150 (MQ=255)
                                 tgtCTTCCTAAATTTCAGCCGGCGGTGTGGTGCCCGCAAACACACCGCCTGGGTTGGGCTTTC‑CCAGCCCTGCCGGCGGGCGCGGCGGGTTTCTTTT‑GCCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAata                                                                                                                    <  1:395192/146‑1 (MQ=255)
                                  gtCTTCCTGAATTTCAGCCAGCGGCGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTaa                                                                                                               >  1:206464/1‑150 (MQ=255)
                                    cTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                            >  1:307987/1‑151 (MQ=255)
                                     ttCCTAATTTCCAGCCGGGGGGGGGGGCCCCGCAACCCCATCCCCTGCGGTTGGCTTTC‑CCAGCCCTGCCGGCAGGCGCAGCAGGGTATTTTT‑GCCCCCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTa                                                                                                                <  2:290146/146‑1 (MQ=255)
                                      tCCTGATTTCCACCCGGGGGTGTGGTCCCCGCTAACACCCCGCCGGGGGTTGGCTTTT‑CCAGCCCTCCCAGCGGGCGCAGCAGGGTATTTTT‑CCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                             <  1:214025/148‑1 (MQ=255)
                                        cTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt                                                                                                        >  2:587280/1‑151 (MQ=255)
                                          gAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                      >  2:507469/1‑151 (MQ=255)
                                             ttccacccgggggTGGGGTCCCCGCTACCCCACCGCCGGCGTTTGGCTTTC‑CCAGCCCTCCCAGCGGGCGCGGCGGGTTTTTTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                      <  2:947571/137‑1 (MQ=255)
                                                 agccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                               >  2:485608/1‑151 (MQ=255)
                                                   ccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                             >  1:635591/1‑151 (MQ=255)
                                                        ggTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                        >  1:408123/1‑151 (MQ=255)
                                                        ggTGTGGTCCCCGCTAACACACCGCCGGGGGTTGGCGTTT‑CCAGACCGGCCAGCGGGCGCAGCGGGTTCTTTTT‑GCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGa                                                                                           <  1:541486/148‑1 (MQ=255)
                                                          tgggggccccGCTAACACACCGCCTGCGGTGGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCGGGTTCCTTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                       <  1:38503/143‑1 (MQ=255)
                                                           gtggtcCCCGCTAACACCTCGCCTGGGTTTGGCTTTC‑CCACCCCTCCCAGCAGGCGCAGCGGGGTCCTTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                       <  1:712174/149‑1 (MQ=255)
                                                             ggtcCCCGCAAACCCACCGCCGGGGGTTGGCTTTC‑CCAGCCCTGCCAGCAGGCGCAGCGGGGTACTTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTc                                                                                    <  1:226317/150‑1 (MQ=255)
                                                             ggcgCCCGCAAACACACCGCCGGCGTTTGGCGTTC‑CCAGCCCCCCCAGCAGGCGCAGAAGGTTCTTTTT‑CCCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt                                                                                     <  2:1003275/146‑1 (MQ=255)
                                                                 cccGCAAACACTCCGCTGGGGTTTGGCTTTC‑CCAGCCCGGCCAGCAGGCGCAGCAGGGTTTTTTT‑GCCCCCACCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgcg                                                                                <  1:353435/150‑1 (MQ=255)
                                                                  ccGCTAACCCATCGCCTGCGGTTGGCTTTC‑CCAGCCCTGCCAGCAGGCGCAGCGGGTTATTTTT‑CCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTa                                                                              <  1:802810/151‑1 (MQ=255)
                                                                   cGCAAACCCACCCCCGGCGGTTGGCTTTC‑CCAGCCCGGCCAGCAGGCGCAGCGGGTTCTTTTT‑GCCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTaa                                                                             <  1:6673/151‑1 (MQ=255)
                                                                       aaCACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                         >  2:958899/1‑151 (MQ=255)
                                                                       aaCACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                         >  1:145895/1‑151 (MQ=255)
                                                                            atcgcggggggTGGGCGTTC‑CCAGCCCTCCCGGCGGGCGCAGCAGGGTACTTTT‑GCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                      <  1:179179/140‑1 (MQ=255)
                                                                            aTCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                    >  2:688194/1‑151 (MQ=255)
                                                                                  ggctgTTGGCTTTC‑CCCGCCCCCCCCGCAGCCGCAGCAGTGTCTTTTC‑GCCGCCCCCGCGGCCCCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                               <  2:743892/146‑1 (MQ=255)
                                                                                   gCGGTTGGCTTTC‑CCAGCCCTGCCAGCGGGCGCAGCAGGGTACTTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                                <  1:748129/148‑1 (MQ=255)
                                                                                     ggTTGGCTTTC‑CCACCCCGCCCAGCAGGCGCAGCGGGGTATTTTT‑CCCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                             <  1:617210/149‑1 (MQ=255)
                                                                                     ggTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                           >  1:557397/1‑151 (MQ=255)
                                                                                     ggTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                           >  2:286871/1‑151 (MQ=255)
                                                                                     ggTTGCCTTTC‑CCAGCCCCGCCAGGGGGCGCGGGAGGGTTCTTTT‑CCCCCCCCCGGGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                             <  2:786603/149‑1 (MQ=255)
                                                                                      gTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                           >  1:415020/1‑150 (MQ=255)
                                                                                             ttt‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaaca                                                    >  2:635972/1‑150 (MQ=255)
                                                                                                  cAGCCCTCCCAGCAGGCGCAGCGGGGTATTTTT‑GCCCCCCCCGGTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                  <  2:824437/148‑1 (MQ=255)
                                                                                                     cccTGCCGGCAGGCGCGGCGGGGTATTTTT‑GCCCCCACCGCGGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                <  1:844160/146‑1 (MQ=255)
                                                                                                      ccTGCCAGCGGGCGCAGCAGGGTTTTTTC‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                            <  2:785814/150‑1 (MQ=255)
                                                                                                                ggCGCAGC‑GGGTTCTTTTCCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAg                                 <  1:962468/151‑1 (MQ=255)
                                                                                                                  cgcAGCGGGTTTTTTTT‑CCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                        <  1:320866/142‑1 (MQ=255)
                                                                                                                  cgcAGCGGGGTTTTTTT‑CCCCCCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                               <  2:280709/151‑1 (MQ=255)
                                                                                                                    cggcgggttccTTTT‑CCCACCCCCTCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGa                             <  1:988887/141‑1 (MQ=255)
                                                                                                                     agcagGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                              >  2:786678/1‑149 (MQ=255)
                                                                                                                       cagGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                           >  2:688756/1‑150 (MQ=255)
                                                                                                                         ggtttttttt‑CCCCCCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGa                             <  1:316772/140‑1 (MQ=255)
                                                                                                                           gTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                      >  1:69799/1‑151 (MQ=255)
                                                                                                                           gTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                      >  1:988932/1‑151 (MQ=255)
                                                                                                                               tttt‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                  >  2:366788/1‑151 (MQ=255)
                                                                                                                               tttt‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                    <  1:507515/149‑1 (MQ=255)
                                                                                                                                tttgcccaccgCCGTTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                    <  1:748210/144‑1 (MQ=255)
                                                                                                                                ttt‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                  <  1:258285/150‑1 (MQ=255)
                                                                                                                                     cccccccggCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATa                 <  2:206449/139‑1 (MQ=255)
                                                                                                                                     ccccccccGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc               <  2:831542/143‑1 (MQ=255)
                                                                                                                                     ccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt             <  2:630210/145‑1 (MQ=255)
                                                                                                                                      cccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg            <  1:644483/146‑1 (MQ=255)
                                                                                                                                      caccacCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg            >  1:669415/1‑151 (MQ=255)
                                                                                                                                       accaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg           >  2:14911/1‑151 (MQ=255)
                                                                                                                                           ccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct        <  1:139531/150‑1 (MQ=255)
                                                                                                                                           ccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctc       <  2:444007/151‑1 (MQ=255)
                                                                                                                                            cgcgGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcc      <  2:539109/151‑1 (MQ=255)
                                                                                                                                              cggCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct     <  1:636023/148‑1 (MQ=255)
                                                                                                                                              cTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt    >  2:794140/1‑151 (MQ=255)
                                                                                                                                               tGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTc   >  1:402293/1‑151 (MQ=255)
                                                                                                                                                 cGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCa  >  2:274010/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCA  >  NC_000913/993768‑994060

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