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→GGGssuBaliphatic 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

GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCA  >  NC_000913/993768‑994060
                                                                                                                                                   |                                                                                                                                                   
gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCccac                                                                                                                                                <  2:848230/151‑1 (MQ=255)
gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCccac                                                                                                                                                >  1:590718/1‑151 (MQ=255)
  aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCccacc                                                                                                                                               <  1:4818/150‑1 (MQ=255)
    cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCaccacc                                                                                                                                            <  2:907553/151‑1 (MQ=255)
    cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCaccacc                                                                                                                                            <  1:832305/151‑1 (MQ=255)
                  gAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGAcc                                                                                                                              >  2:208384/1‑151 (MQ=255)
                    accatcatCCGCTTTTCTTCCTAATTTCCACCCAGCGGTGGGGTCCCCGTAAACACACCCCCTGGGTTTGGCTTTTCCAGCCCGGCCAGCAGGCGCAGCAGGGTATTTTCCCCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTg                                                                                                                            <  1:743092/149‑1 (MQ=255)
                    aaCATCTTCCGCGTTTCTTCCGAAATTTCAGCCGGGGGGGGGGTCCCCGCAACCACACCGCGGGGGTTGGCCTTTCCCAGCCCGGCCAGCAGGCGCAGAGGGGTTCTTTTGCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTg                                                                                                                            <  2:797880/151‑1 (MQ=255)
                      catcatCCGCGTGTCTTCCTGAATTTCAGCCGGCGGTGTGGGGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGCCCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGc                                                                                                                           <  1:716670/150‑1 (MQ=255)
                       atcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCg                                                                                                                         >  1:779366/1‑151 (MQ=255)
                        tcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGTTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                        >  2:417066/1‑151 (MQ=255)
                          atTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGa                                                                                                                       >  1:280744/1‑150 (MQ=255)
                                 tgtCTTCCTAAATTTCAGCCGGCGGTGTGGTGCCCGCAAACACACCGCCTGGGTTGGGCTTTCCCAGCCCTGCCGGCGGGCGCGGCGGGTTTCTTTTGCCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAata                                                                                                                    <  1:395192/146‑1 (MQ=255)
                                  gtCTTCCTGAATTTCAGCCAGCGGCGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTaa                                                                                                               >  1:206464/1‑150 (MQ=255)
                                    cTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                            >  1:307987/1‑151 (MQ=255)
                                     ttCCTAATTTCCAGCCGGGGGGGGGGGCCCCGCAACCCCATCCCCTGCGGTTGGCTTTCCCAGCCCTGCCGGCAGGCGCAGCAGGGTATTTTTGCCCCCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTa                                                                                                                <  2:290146/146‑1 (MQ=255)
                                      tCCTGATTTCCACCCGGGGGTGTGGTCCCCGCTAACACCCCGCCGGGGGTTGGCTTTTCCAGCCCTCCCAGCGGGCGCAGCAGGGTATTTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                             <  1:214025/148‑1 (MQ=255)
                                        cTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt                                                                                                        >  2:587280/1‑151 (MQ=255)
                                          gAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                      >  2:507469/1‑151 (MQ=255)
                                             ttccacccgggggTGGGGTCCCCGCTACCCCACCGCCGGCGTTTGGCTTTCCCAGCCCTCCCAGCGGGCGCGGCGGGTTTTTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                      <  2:947571/137‑1 (MQ=255)
                                                 agccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                               >  2:485608/1‑151 (MQ=255)
                                                   ccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                             >  1:635591/1‑151 (MQ=255)
                                                        ggTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                        >  1:408123/1‑151 (MQ=255)
                                                        ggTGTGGTCCCCGCTAACACACCGCCGGGGGTTGGCGTTTCCAGACCGGCCAGCGGGCGCAGCGGGTTCTTTTTGCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGa                                                                                           <  1:541486/148‑1 (MQ=255)
                                                          tgggggccccGCTAACACACCGCCTGCGGTGGGCGTTTCCAGACCTGCCAGCAGGCGCAGCGGGTTCCTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                       <  1:38503/143‑1 (MQ=255)
                                                           gtggtcCCCGCTAACACCTCGCCTGGGTTTGGCTTTCCCACCCCTCCCAGCAGGCGCAGCGGGGTCCTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                       <  1:712174/149‑1 (MQ=255)
                                                             ggtcCCCGCAAACCCACCGCCGGGGGTTGGCTTTCCCAGCCCTGCCAGCAGGCGCAGCGGGGTACTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTc                                                                                    <  1:226317/150‑1 (MQ=255)
                                                             ggcgCCCGCAAACACACCGCCGGCGTTTGGCGTTCCCAGCCCCCCCAGCAGGCGCAGAAGGTTCTTTTTCCCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt                                                                                     <  2:1003275/146‑1 (MQ=255)
                                                                 cccGCAAACACTCCGCTGGGGTTTGGCTTTCCCAGCCCGGCCAGCAGGCGCAGCAGGGTTTTTTTGCCCCCACCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgcg                                                                                <  1:353435/150‑1 (MQ=255)
                                                                  ccGCTAACCCATCGCCTGCGGTTGGCTTTCCCAGCCCTGCCAGCAGGCGCAGCGGGTTATTTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTa                                                                              <  1:802810/151‑1 (MQ=255)
                                                                   cGCAAACCCACCCCCGGCGGTTGGCTTTCCCAGCCCGGCCAGCAGGCGCAGCGGGTTCTTTTTGCCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTaa                                                                             <  1:6673/151‑1 (MQ=255)
                                                                       aaCACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                         >  2:958899/1‑151 (MQ=255)
                                                                       aaCACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                         >  1:145895/1‑151 (MQ=255)
                                                                            atcgcggggggTGGGCGTTCCCAGCCCTCCCGGCGGGCGCAGCAGGGTACTTTTGCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                      <  1:179179/140‑1 (MQ=255)
                                                                            aTCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                    >  2:688194/1‑151 (MQ=255)
                                                                                  ggctgTTGGCTTTCCCCGCCCCCCCCGCAGCCGCAGCAGTGTCTTTTCGCCGCCCCCGCGGCCCCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                               <  2:743892/146‑1 (MQ=255)
                                                                                   gCGGTTGGCTTTCCCAGCCCTGCCAGCGGGCGCAGCAGGGTACTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                                <  1:748129/148‑1 (MQ=255)
                                                                                     ggTTGGCTTTCCCACCCCGCCCAGCAGGCGCAGCGGGGTATTTTTCCCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                             <  1:617210/149‑1 (MQ=255)
                                                                                     ggTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                           >  1:557397/1‑151 (MQ=255)
                                                                                     ggTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                           >  2:286871/1‑151 (MQ=255)
                                                                                     ggTTGCCTTTCCCAGCCCCGCCAGGGGGCGCGGGAGGGTTCTTTTCCCCCCCCCGGGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                             <  2:786603/149‑1 (MQ=255)
                                                                                      gTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                           >  1:415020/1‑150 (MQ=255)
                                                                                             tttCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaaca                                                    >  2:635972/1‑150 (MQ=255)
                                                                                                  cAGCCCTCCCAGCAGGCGCAGCGGGGTATTTTTGCCCCCCCCGGTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                  <  2:824437/148‑1 (MQ=255)
                                                                                                     cccTGCCGGCAGGCGCGGCGGGGTATTTTTGCCCCCACCGCGGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                <  1:844160/146‑1 (MQ=255)
                                                                                                      ccTGCCAGCGGGCGCAGCAGGGTTTTTTCCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                            <  2:785814/150‑1 (MQ=255)
                                                                                                                ggCGCAGCGGGTTCTTTTCCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAg                                 <  1:962468/151‑1 (MQ=255)
                                                                                                                  cgcAGCGGGTTTTTTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                        <  1:320866/142‑1 (MQ=255)
                                                                                                                  cgcAGCGGGGTTTTTTTCCCCCCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                               <  2:280709/151‑1 (MQ=255)
                                                                                                                    cggcgggttccTTTTCCCACCCCCTCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGa                             <  1:988887/141‑1 (MQ=255)
                                                                                                                     agcagGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                              >  2:786678/1‑149 (MQ=255)
                                                                                                                       cagGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                           >  2:688756/1‑150 (MQ=255)
                                                                                                                         ggttttttttCCCCCCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGa                             <  1:316772/140‑1 (MQ=255)
                                                                                                                           gTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                      >  1:69799/1‑151 (MQ=255)
                                                                                                                           gTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                      >  1:988932/1‑151 (MQ=255)
                                                                                                                               ttttGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                  >  2:366788/1‑151 (MQ=255)
                                                                                                                               ttttCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                    <  1:507515/149‑1 (MQ=255)
                                                                                                                                tttgcccaccgCCGTTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                    <  1:748210/144‑1 (MQ=255)
                                                                                                                                tttCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                  <  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)
                                                                                                                                                   |                                                                                                                                                   
GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCA  >  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: