Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913501,5640AC39.7% 114.7 / 16.6 68*559E (TAA→GAA) ybaLinner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (33/8);  new base C (0/27);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 3.62e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATC  >  NC_000913/501416‑501709
                                                                                                                                                    |                                                                                                                                                 
aCTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGATCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCTTGTTTACTTATTTTTCATCTTTTACATCTGGGAACCACCGATTa                                                                                                                                                   >  2:708553/1‑149 (MQ=37)
   gggATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCg                                                                                                                                              <  2:630890/151‑1 (MQ=255)
     gATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTc                                                                                                                                            <  2:832172/151‑1 (MQ=255)
         gTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcacca                                                                                                                                        <  1:1148899/151‑1 (MQ=255)
          ttAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcaccac                                                                                                                                       >  1:492990/1‑151 (MQ=255)
           tAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcaccac                                                                                                                                       <  2:1239546/150‑1 (MQ=255)
           tAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCaccacc                                                                                                                                      >  1:181650/1‑151 (MQ=255)
             aGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTc                                                                                                                                    >  1:587595/1‑151 (MQ=255)
                 aTCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCAccc                                                                                                                                >  2:721555/1‑151 (MQ=255)
                 aTCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCAccc                                                                                                                                <  2:1116023/151‑1 (MQ=255)
                   cgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccg                                                                                                                               <  1:1135573/150‑1 (MQ=255)
                       ggTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGt                                                                                                                          >  2:714654/1‑151 (MQ=255)
                             tAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGtt                                                                                                                    >  1:278214/1‑151 (MQ=255)
                                       aaGGCTTCCCCGCCTCATCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCGGGCGCTTTTTATCATCCCTACTTTTTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGt                                                                                                           <  1:405144/150‑1 (MQ=255)
                                        gggttttCCCGCCGTATCTGGCGTTTGGCTCTCGTTGCTTGAGCCGGCGCGGGCGCTTTTTTTCATCCCTACATTTTTTTAATTTTTTCCACCCGGAACCCACCGTTTCCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGt                                                                                                           <  1:443596/144‑1 (MQ=255)
                                             ttCCCGCCGCATCTGGCATTTGCCTCTCGATGCTTGATGCGCCGCTGGCGCGTTTTATCATCCCTACATTTTTTTCATATTTTACACCCGCCACCCACCGTTTACCCCGCCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCa                                                                                                       <  1:137874/148‑1 (MQ=255)
                                              tCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                                                                                   >  2:1172975/1‑151 (MQ=255)
                                               cccgccgcATCTGGCTTTTGGCTCTCGATGCCTGAGGCGCCGCTGGCGCGTCTTATCATGCCAACATATTTTTCATTTTTTACCCCCGGAAACCCCCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                                                                                    <  2:456362/147‑1 (MQ=255)
                                               cccgccgcATCTGGCTTTTGGCCCTCGATGCCTGAGGCGACGCGGGCGCGTCTTTTCATGCCTACATATTTTTCATATTTTACATCCGGCACCCCCCTTTTCCCCCTCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                                                                                    <  1:300586/147‑1 (MQ=255)
                                               cccgccccTTCTGGCTTTTGCCTCCCGATCCCTGATGCGCCGCGGGCGCGTCTTACAATCCCTACATTTTTTTAATTTTTTACACCCGGCACCCACCTTTTCCCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                                                                                   <  1:831244/148‑1 (MQ=255)
                                               ccccccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGCCGCGGGCGCGTTTTATCATCCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCa                                                                                                       <  2:1124558/142‑1 (MQ=255)
                                               cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATg                                                                                                  >  2:39372/1‑151 (MQ=255)
                                               cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATg                                                                                                  >  2:217156/1‑151 (MQ=255)
                                                 cgccgcCTCTGCCATTGGGCTCCGGATGCGTGATGCGACGCGGGCGCGTCTTATCATCCCTACATTTTTTTTTTATTTCCCCCCCGGAAACCCCCTTTTCCCCCGTCCCCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCa                                                                                                       <  2:220812/144‑1 (MQ=255)
                                                    cgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGGCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCCTGGTACg                                                                                             >  1:956050/1‑151 (MQ=255)
                                                        tCTGGCTTTTGGCTCTCGATGCCTGATGCGCCGCGGGCGCGTCTTATCATGCCTACATATTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGc                                                                                          <  2:412478/150‑1 (MQ=255)
                                                          tGGCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCTTCTTATCATGCCTACATATTTTTCATTTTTTACATCCGGCACCCCCCGTTTCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACgg                                                                                            <  2:182451/146‑1 (MQ=38)
                                                           ggCTTTTGGCTCCCGATGCCGGATGCGACGCTGGCGGTTCTTATCATGCCTACATATTTTTCATATTTTACATCCGACACCCCCCGTTTCCCCCTTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGc                                                                                          <  2:566599/147‑1 (MQ=255)
                                                              ttttGGCTCTCGATGCCTGACGCTGCCCGGGCGCTTCTAATCATGCTTACTTTTTTTTCAAATTTTAGACCCGGCACCCACCTTTGCCCCGGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAttt                                                                                     <  1:1101080/148‑1 (MQ=255)
                                                                    ctctCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCg                                                                             >  1:687694/1‑151 (MQ=255)
                                                                          aTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGTTACGGGCGATTTCACGCTCGCCCAt                                                                        >  1:663159/1‑150 (MQ=255)
                                                                            gCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                                     >  2:986349/1‑151 (MQ=255)
                                                                             cttgagcggccgcgGGCGTGTTTTATCTTCCCTACTTTTTTTCAATATTTCCCATCCGTCACCCCCCGTTTCCCCCGCCCCCCCCCACCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                                     <  1:943824/140‑1 (MQ=255)
                                                                             ccTGATGCGACGCTGGCGCGACTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCTACGGAGGTGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCCCGCTCGCCCATCACt                                                                    >  1:568360/1‑151 (MQ=255)
                                                                               tggtGCGCCGCTGGCGCTTTTTACCATCCCTACTTTTTTTTAATATTTTCCACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                                    <  1:644663/146‑1 (MQ=255)
                                                                                        gCTGGCGCCTTTTATCATGCCTACATTTTTTTCATTTTTTACCTCCGGCACCCCCCTTTTCCCCCTTCCCCCCCCCGCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCt                                                                <  1:972685/144‑1 (MQ=255)
                                                                                         cTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGc                                                         >  2:1072247/1‑150 (MQ=255)
                                                                                           ggCGCGTCTTATCTTCCCTATATTTTTTCCATTTTTTCCACCCGCCACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAt                                                             <  1:259801/144‑1 (MQ=255)
                                                                                           ggCCCGTTTTACCATGCCTCTATTTTTTTAATTTTTTCCACCCGGCACCCCCCTTTCCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAt                                                             <  1:1044512/144‑1 (MQ=255)
                                                                                                gTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                 >  2:626310/1‑151 (MQ=255)
                                                                                                  cTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                               >  1:688409/1‑151 (MQ=255)
                                                                                                   ttaccttcccTACATTTTTTTCATTTTTCACACCCGGAACCCACCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGca                                                        <  1:932723/133‑1 (MQ=255)
                                                                                                   ttATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg                                              >  2:181144/1‑151 (MQ=255)
                                                                                                     atctgcccaactttttttcCATTTTTAACACCCGGCACCCCCCTTTTCCCCCCCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg                                              <  1:777420/135‑1 (MQ=255)
                                                                                                        aTGCCTACATTTTTTTCATTTTTTACTCCCGGCACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                             <  2:701056/147‑1 (MQ=255)
                                                                                                         tGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                        >  2:7884/1‑151 (MQ=255)
                                                                                                           ccTACATTTTTTTCATTTTTTCCACCCGGCACCCCCCTTTTCCCCCGTCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                        <  2:225262/149‑1 (MQ=255)
                                                                                                             tACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc                                    >  1:1012948/1‑151 (MQ=255)
                                                                                                             tACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGAGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc                                    >  1:1053106/1‑151 (MQ=255)
                                                                                                               cATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGcc                                   >  1:37567/1‑150 (MQ=255)
                                                                                                                 tatTTTTCATATTTTACACCCGCAAACCCCCGTTTCCCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                      <  2:597938/145‑1 (MQ=255)
                                                                                                                  atTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGACGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACt                                <  2:33672/150‑1 (MQ=255)
                                                                                                                      ttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATc                           >  1:1139500/1‑151 (MQ=255)
                                                                                                                       tCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCg                          >  2:97459/1‑151 (MQ=255)
                                                                                                                           tttttCCACCCGCCACCCCCCTTTTCCCCCCCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGt                         <  1:832352/147‑1 (MQ=255)
                                                                                                                           atTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCa                       >  1:519419/1‑150 (MQ=255)
                                                                                                                           atTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                      >  1:623379/1‑151 (MQ=255)
                                                                                                                            ttttACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGt                         <  2:10691/147‑1 (MQ=255)
                                                                                                                            ttttACACCCGGCACCCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTc                        <  2:587624/148‑1 (MQ=255)
                                                                                                                              ttACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt                   >  1:862851/1‑151 (MQ=255)
                                                                                                                              ttACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt                   >  1:994807/1‑151 (MQ=255)
                                                                                                                               tACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATg                  >  2:578832/1‑151 (MQ=255)
                                                                                                                                 cgcccGCCACCCCCCTTTTCCCCCCCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATggg                <  1:768148/148‑1 (MQ=255)
                                                                                                                                 cccccGGCACCCCCCTTTTCCCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATggg                <  1:641292/148‑1 (MQ=255)
                                                                                                                                 cATCCGGCAACCACCGGTTACCCCGTCAGCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATggg                >  1:1176786/1‑151 (MQ=255)
                                                                                                                                  aTCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgc               >  2:1206459/1‑151 (MQ=255)
                                                                                                                                       gCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCggg          >  2:1170729/1‑151 (MQ=255)
                                                                                                                                                gTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtc  >  1:148559/1‑150 (MQ=255)
                                                                                                                                                 tttCCCCCGCCCCTACCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtc  <  2:740214/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATC  >  NC_000913/501416‑501709

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