Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091317,0500AC45.8% 110.6 / 13.8 75intergenic (‑90/‑439)mokC/nhaAregulatory protein for HokC, overlaps CDS of hokC/sodium‑proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (29/10);  new base C (0/33);  total (29/46)
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

CTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAG  >  NC_000913/16906‑17196
                                                                                                                                                  |                                                                                                                                                  
ctctctTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACccgcc                                                                                                                                              >  2:993887/1‑151 (MQ=255)
 tctctTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACccgcc                                                                                                                                              >  2:115758/1‑150 (MQ=255)
   tcttctCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCt                                                                                                                                           >  2:396364/1‑151 (MQ=255)
         ccTTGCCCTTAGGGTCAGAAAGGGGCCCTCTAAATGTTTCCAGAATATGGGGGGCCTGGGGTTGATGATAAAATAT‑CACGCGGGGTTTTTCTCTCTCCGCCT‑TTCAAATATGGCCTGAGACAGACAGCCTCAAGCCCCCGCCGCtatta                                                                                                                                       <  1:512992/149‑1 (MQ=255)
          cTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTatat                                                                                                                                    >  2:189376/1‑151 (MQ=255)
           ttGACCTTACGGTCAGTAAGAGCCACTCTACATGTGTTCAGCATATGGGGGGCCTGGGTTTGATGTAAAAATAT‑CCCTCGGGTTTTTTCTCTACCGCCCT‑TTAAGTAAATGCCTAAGACAGCCACCCTCAAGCCCCCGCCGCtattat                                                                                                                                      <  1:332168/148‑1 (MQ=255)
                 ttACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGctctc                                                                                                                             >  2:640064/1‑151 (MQ=255)
                    cGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCttt                                                                                                                          >  1:920027/1‑151 (MQ=255)
                     ggTCAGTAAGAGGCCCTCTACATGTTTTCAGCTTATAGGAGCCCTGGGGTTGATGTAAAAATAC‑CCCTGGGGTCTTTTCTCTACCTCCCG‑TTCAGCTAATGCCTGAGACAGACACCCACAACCCCCCGCCGCTATTATATCGctct                                                                                                                              <  2:992615/146‑1 (MQ=255)
                     ggTCAGTAAGAGGCCCTCTACATGTTTTCAGCATATAGGAGGCCTCGGGTGGATGAAAAAATAT‑CCCGCGGGTTTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCCCCCGCCGCTATTATATCGCTCTCtt                                                                                                                           <  2:491421/149‑1 (MQ=255)
                     ggTAAGTAAGAGGCCCTCTACATGTTTTCAGCATATAGGAGGCCGGGGGTTGATGGAAAAATAT‑CCCTCGGGGTTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACACCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTa                                                                                                                         <  2:920077/151‑1 (MQ=255)
                        cAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAAcc                                                                                                                      >  1:844172/1‑151 (MQ=255)
                         aGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAAccc                                                                                                                     >  1:402708/1‑151 (MQ=255)
                            aaGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCAtt                                                                                                                  >  1:481294/1‑151 (MQ=255)
                            aaGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCAtt                                                                                                                  >  2:607828/1‑151 (MQ=255)
                               agGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTg                                                                                                               >  1:315195/1‑151 (MQ=255)
                               agGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTg                                                                                                               >  1:556065/1‑151 (MQ=255)
                                ggCCCTCTACTTGTTTTCAGAAAATAGGGGGCCTGGGTTTGATGAAAAAATAT‑CCCTGGGGTTTTTTCTCTATCTGCCG‑TTCAGCTAATGCCGGAGAAAGACAGCCTAAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTg                                                                                                               <  2:4221/150‑1 (MQ=255)
                                  ccctctACATTTGTTCAGCATATGGGAGCCCTGGGGTTGATGTAAAAATAT‑CCCGGGGGTTTTTTCTCTATCTGCCT‑TTCAGCTAATGCCTGAGACAGACACCCTAAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGttt                                                                                                            <  2:312530/149‑1 (MQ=255)
                                        aCATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTGTTATATCGCTCTCTTTAACCCATTTTGTTTTATCg                                                                                                       >  2:471086/1‑150 (MQ=255)
                                           tgttTTAAGCATATAGGAGGCCTGGGTTTGATGAAAAAATAC‑CCCTCGGGTTTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACACCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGAt                                                                                                     <  2:732704/149‑1 (MQ=255)
                                             tgtTCAGAATATGGGGGGCCTGGGGTTGATGAAAAAATAT‑CCCGGGGGTTTTTTCTCTTCCTCCCT‑TTCAGCTAATGCCTGAGACAGACACCCTAAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGAtt                                                                                                    <  1:195578/148‑1 (MQ=255)
                                             tgtTCACCATAGAGGGGCCCTCGGTTTGATGAAAAAATCT‑CCCTCGGGTTTTTTCTCTATCTGCCT‑TTCATCTAATCCCGGAGACAGACACCCTAAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGtt                                                                                                             <  1:557064/139‑1 (MQ=255)
                                             tgtTAAGCATATGGGAGGCCTGGGTTTGATGAAAAAATAT‑CCCTGGGGTTTTTTCTCTATCTCCCT‑TTCAGTAAATGCCTGAGACAGACAGCCTCAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCt                                                                                                  <  1:262493/150‑1 (MQ=255)
                                               ttCAGCATATAGGAGGCCTGGGGTTGATGAAAAAATAC‑CCCGGGGGTTTTTTCTCTATCTGCCT‑TTCAGCTAATGCCTGAGAAAGACACCCTCAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTa                                                                                                 <  2:556099/149‑1 (MQ=255)
                                               ttCAGAATATAGGAGGCCGGGGTTTGATGAAAAAATTT‑CACGGGGGTTTTTTCTCTTTCTCCCG‑TTCAACTAATCCCTAAAAAAGCCAGCTTCAATCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCg                                                                                                       <  1:985249/143‑1 (MQ=255)
                                                 caacATAAAGGAGGCCTGGGGTTGATGGTAAAATTC‑CCCGCGGGTTTTTTCTCTATCGGCCT‑TTCAGTTAATCCGTGAGAAAAACCGCCTCAACCGCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATc                                                                                              <  2:1054411/147‑1 (MQ=255)
                                                   atataggggaggCCGCGGGTTTTTAGAAAAATAT‑CCCGCGGGTCTTTTCTCTTCCTCCGG‑GTCATTTAGGGCCAGAACAAGAGACCCGAAATCGCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTg                                                                                           <  2:901647/144‑1 (MQ=255)
                                                       tatgggtggccttggTTGTAGGGAAAATATACACGCGGTGTTTTTCTCTATCTGCTT‑TTCATCAACTGCCGGAAACACACCTTCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTGAACCCATTTTGTTTTATCGATTCTAATCCt                                                                                            <  1:594484/134‑1 (MQ=255)
                                                            aggCCTCGGGTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCt                                                                                  >  2:342695/1‑151 (MQ=255)
                                                             ggCCTCGGGTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTc                                                                                 >  1:782111/1‑151 (MQ=255)
                                                               cccggggTTGATGAAAAAATAT‑CCCCCGGGGTTTTTTTCTTTCTCCTT‑TCCAGCTAAGGCCTGAGACAGCCACCCTAAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGAc                                                                                      <  1:490789/140‑1 (MQ=255)
                                                               ccTGGGGTTGATGTAAAAATAC‑CCCCCGGGTTTTTTCTTTATCTCCCT‑TTCAGCTAATCCCTGAAACAGACACCCTAAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCg                                                                                <  1:640023/150‑1 (MQ=255)
                                                                ctggggTTGATGAAAAAATAC‑CACTCGGGTTTTTTTCCCCTCTGCTT‑TACAATTAATCCCGGAGACAAACACCCCCACCCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCt                                                                                  <  2:280989/144‑1 (MQ=255)
                                                                  cGGGTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCAtt                                                                            >  1:815464/1‑151 (MQ=255)
                                                                   ggtttGATGGAAAAATAC‑CCCGGGGGGCTTTTCTTTATC‑GCCGTTTCAGCTACTCCCCAAGACGGCCACCCTCACCCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCAt                                                                             <  2:215593/146‑1 (MQ=255)
                                                                   gggTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCAttt                                                                           >  2:492348/1‑151 (MQ=255)
                                                                   gggTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCAttt                                                                           >  2:998308/1‑151 (MQ=255)
                                                                   gggTGGATGTAAAAATAT‑CACCCGGGTTTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACACCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCAttt                                                                           <  2:434366/151‑1 (MQ=255)
                                                                    gtttaaggaaaaaaTAC‑CACGCGGGGTTTTTCTCTTTCTCCCG‑TTCAGAAAATGCCGGAAACAGCCCCCCAAAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACACCt                                                                                  <  2:121776/133‑1 (MQ=255)
                                                                     gTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCAttttt                                                                         >  2:407416/1‑151 (MQ=255)
                                                                      tggaggaaaaaaTAT‑CCCGGGGGGTTTTTCTCTACCTGCCG‑TTCAGCTATTGCCTGAGAAAGCCACCCACAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCAtttt                                                                          <  2:445828/141‑1 (MQ=255)
                                                                       tGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTtgt                                                                       >  1:163476/1‑151 (MQ=255)
                                                                          tGGAAAAAAAT‑CCCCGGGGGCTTTTCTCTTTCTGCCT‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTtgtg                                                                      <  2:518334/149‑1 (MQ=255)
                                                                             tAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGt                                                                  >  1:585146/1‑150 (MQ=255)
                                                                              aaaaTAT‑CCCCGGGGTTTTTTCTCTTTCTGCCT‑TTCAGCTAATGCCTGAGACAGACACCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGt                                                                  <  2:809488/149‑1 (MQ=255)
                                                                              aaaaTAA‑CCCGCGGGTTTTTTCTCTATCTGCCG‑TTCAGCTATTGCCGGAGACAGACGCCCTAAGGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGt                                                                  <  2:829125/149‑1 (MQ=255)
                                                                                      ccctggggttttttCTCTTTCTGCCT‑TTCAGCTAATGCCGGAGAAAAACACCTTAAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAAttt                                                             <  1:189357/137‑1 (MQ=255)
                                                                                      cACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAAtttttt                                                          >  2:1041360/1‑150 (MQ=255)
                                                                                           gggttttttCTCTACCTCCCG‑TTCAGCTATGGCCTGAAACAGACATCTTAAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAAtttttt                                                          <  1:251449/140‑1 (MQ=255)
                                                                                           ggggTTTTTCTCTTTCTCCCT‑TTCAGCTAAGCCCTGAGACAACCAGCCACAACCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATg                                                      <  1:656163/149‑1 (MQ=255)
                                                                                                  ttttctACCTGCCG‑TTCAGCAAATGCCTGAGACAGAAACCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAatt                                              <  2:450052/147‑1 (MQ=255)
                                                                                                  ttCTCTATCTGCGT‑TTCACTAAATCCCTGAGACAAACGCCCCCAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAatt                                              <  1:378141/150‑1 (MQ=255)
                                                                                                     tctATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTAttt                                          >  2:182382/1‑151 (MQ=255)
                                                                                                       tATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTa                                         >  1:350399/1‑150 (MQ=255)
                                                                                                        aTCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAAc                                       >  1:1044569/1‑151 (MQ=255)
                                                                                                             ccG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTa                                   >  1:356403/1‑150 (MQ=255)
                                                                                                             ccG‑TTCAGCTAAGCCCGGAGACAGCCACCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAattatt                                           <  2:481314/142‑1 (MQ=255)
                                                                                                              cctttcactaaaTCCCGGAGACACACACCAAAAAGCGCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACttt                                    <  1:1031158/139‑1 (MQ=255)
                                                                                                              cctttCATCTAATCCCAGAGACACACACCCTAAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAattatt                                           <  1:128485/139‑1 (MQ=255)
                                                                                                              ccattaaGCAAATGCCTGAGAAAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACttt                                    <  2:268284/143‑1 (MQ=255)
                                                                                                              cG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAAt                                 <  2:388119/151‑1 (MQ=255)
                                                                                                                    aGCAAATGCCTGAGACAAACACCCTAAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTa                              <  2:566347/149‑1 (MQ=255)
                                                                                                                      cTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTAtctc                          >  1:1072575/1‑151 (MQ=255)
                                                                                                                            ccTGATACACACCCCCTAAACCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATc                    <  1:165626/151‑1 (MQ=255)
                                                                                                                               gagaCAGACACCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCa                 <  1:308394/151‑1 (MQ=255)
                                                                                                                                agaaagccAGCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTc                       <  2:596410/144‑1 (MQ=255)
                                                                                                                                  acaaacaccctacaccACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAatt              <  2:252658/136‑1 (MQ=255)
                                                                                                                                    agccacccccAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAatt              <  2:148047/140‑1 (MQ=255)
                                                                                                                                        cgCCCCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTga         <  1:41990/149‑1 (MQ=255)
                                                                                                                                        agCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTgac        >  2:931989/1‑151 (MQ=255)
                                                                                                                                           cTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTgacgac     >  2:147070/1‑151 (MQ=255)
                                                                                                                                              actcACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCTCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTgacga      <  1:231326/144‑1 (MQ=255)
                                                                                                                                              aaccACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAg  <  2:239858/148‑1 (MQ=255)
                                                                                                                                              aaGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAg  <  1:454553/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATAT‑CACTCGGGGCTTTTCTCTATCTGCCG‑TTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAG  >  NC_000913/16906‑17196

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: ‑