Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,622,8160GA44.3% 32.3 / 61.6 62intergenic (+75/+144)yneM/dgcZinner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base G (34/0);  new base A (0/27);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 6.33e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAAC  >  NC_000913/1622668‑1622954
                                                                                                                                                       |                                                                                                                                          
ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGa                                                                                                                                           >  1:1017752/1‑150 (MQ=255)
     gCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGaaaaaa                                                                                                                                      >  1:1154726/1‑150 (MQ=255)
       cGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAAc                                                                                                                                   >  2:931167/1‑151 (MQ=255)
       cGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAAc                                                                                                                                   >  2:853574/1‑151 (MQ=255)
             gggAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGcc                                                                                                                               >  1:455940/1‑149 (MQ=255)
                      tttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcggc                                                                                                                    >  1:989775/1‑151 (MQ=255)
                         attttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcggcgg                                                                                                                  >  1:1109413/1‑150 (MQ=255)
                          ttttttgtttttctgcCCGCTTATTCCCCCAACAAAGGGGATGACAAATAAACGGAGATACCCCCTCACCTACCCGCCCCCTTTTAAC‑GGTGGTGA‑ACAGGGGGCT‑TTTTTTTTATGACGACAAAAAAAAACCGCCAGTAAACcggcgg                                                                                                                  <  1:925437/149‑1 (MQ=255)
                               tggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGc                                                                                                           >  1:85667/1‑151 (MQ=255)
                               tgggtttctggCCGCGCATTTCGGCTACAACTGGGCTGACTAATGAACGGAGATAATCCCTCATATAACCGGCACCTTGTTAC‑TGCTGTGT‑AGAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGACTg                                                                                                            >  1:974403/1‑150 (MQ=255)
                                ggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGCGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCCGCGgagactttt                                                                                                          >  1:615889/1‑146 (MQ=255)
                                   tttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                       >  1:302945/1‑151 (MQ=255)
                                   tttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                       >  2:610450/1‑151 (MQ=255)
                                    ttctggCCGCTTATTTCAGCAAAAAAGGGGATAATAAAGGAAGGGAAATAACCCCCCCCCTACCCGCCCCTTTTTTAC‑ATTTGTGT‑AAAGGGGGCT‑TATTTTTTAGGGGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                       <  2:351407/150‑1 (MQ=255)
                                        ggccgtttttttCACCCAAAAGGGGGATAACAAATGAGGGGAAAAAACCCCCCCCCAACCCGCCCCTTTGTTAC‑AGTTGTGT‑AAAGGGGGCC‑TGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCa                                                                                                      <  2:139201/138‑1 (MQ=255)
                                        gccggctttttcccgCCAAAAAGGGGATGATAAAGAAAGGGAAATAACCCCCCCCCAACCCGGCCCTTTTTTAC‑GGTTGTGA‑AAAGGGGGCT‑TTTTTTTTAGGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCAt                                                                                                     <  1:154331/134‑1 (MQ=255)
                                          ccGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGa                                                                                                  >  2:450685/1‑149 (MQ=255)
                                           cgcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TCATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                 >  1:276204/1‑149 (MQ=255)
                                            gcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                 >  1:161732/1‑148 (MQ=255)
                                              tttttccacccAAAAATGGTATGATAAATAAAGGAAAATACCCCCCCCCCAACCCGCCCCTTTTTTAC‑AGTTGTGA‑AAAGGGGGCT‑TGTTTTTTAGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCa                                                                                                      <  1:176294/132‑1 (MQ=255)
                                               tttttCAGCCAAAAGGGGGATGACTAATGAGGGGAGATACCCCCCCCCCTACCCGCCCCTTTTTTAC‑AGTTGTGT‑ACAGGGGGCT‑TTTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt                                                                                            <  1:68712/148‑1 (MQ=255)
                                                 tttCACCCAAAAAGGGGATGACAAATGAAGGAAGATACCCCCCCCCCAACCCGCCCCTTTTTTAC‑ATTTGTG‑‑AAAAGGGGCCTTTTTTTTTAGGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAg                                                                                               <  2:631995/145‑1 (MQ=255)
                                                  ttCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTtgt                                                                                         >  2:569836/1‑150 (MQ=255)
                                                     aGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGttt                                                                                     >  1:879543/1‑151 (MQ=255)
                                                     aGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAATTGTGttt                                                                                     >  1:436597/1‑151 (MQ=255)
                                                      gCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGtttt                                                                                    >  1:912871/1‑151 (MQ=255)
                                                         aaaaaGGGGATGAAAAAGGAAGGGAAATAACCCCCCCCCTCCCCGCCCCTTTTTTAC‑ATTTGTGT‑AAAGGGGGCT‑TTTTTTTTATGACGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGtttt                                                                                    <  2:557724/146‑1 (MQ=255)
                                                           aaaTGGGATGACTAATGAACGGAGATAATACCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGC‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGAGGAGAATGATTGCATGGATAGATTTGTGATTTGCtt                                                                                >  1:1083116/1‑150 (MQ=255)
                                                           aaaGGGGATGATAAATGAAGGGAAAAACCCCCTCCCCAACCCGCCCCTTTTTTAC‑AGTGGT‑T‑AAAAGGGGCCTTTTTTTTTATGCGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtt                                                                                <  1:3543/150‑1 (MQ=255)
                                                             aTGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGGTGATGCTTGCATGGGTAGATTTTTGTTTTTCtttt                                                                              >  1:780176/1‑150 (MQ=255)
                                                                         aTGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAGAGATTTGTGTTTTGCTTTTACGCCAACAGGc                                                                  >  1:530438/1‑150 (MQ=255)
                                                                           gAACGGAGATAA‑CCCCCCCCAAACCCGCCCCTTTTTCA‑GGGTGTATCGCAGGGGTGC‑TGTTTTTAATGACGCAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAgg                                                                   <  2:234801/147‑1 (MQ=255)
                                                                            aggggAAATAACCCCTCACC‑AACCCGCCCCTTTTTTCGAGTTGTGT‑AAAGGGGGCT‑TGTTTTTTAGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAt                                                                <  1:938565/146‑1 (MQ=255)
                                                                             aCGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtttt                                                             >  2:1011332/1‑151 (MQ=255)
                                                                                 aaaaacccccTCCCCTACCCGCCCCTTTTTTAC‑GGTTGTGT‑AAAGGGGGCT‑TTTTTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTc                                                            <  1:704510/141‑1 (MQ=255)
                                                                                     acccccTCACAAACCCGCCCCTTTGTTAC‑GGTTGTGA‑ACAGGGGGCT‑TTTTTTTTAGGCGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCAc                                                      <  2:638596/147‑1 (MQ=255)
                                                                                        cccccaccaccccGCCCCTTTGTAAC‑TGTTGTAT‑AAAGGGGTCT‑TTTTTTTTATGGGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTg                                                    <  2:1037986/138‑1 (MQ=255)
                                                                                          cTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATa                                                 >  2:483876/1‑150 (MQ=255)
                                                                                            cACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGTCGATAAAAAACCGCCAGTAAACCGGCGGCTAATGCCTGCATGGATAGATTTGTGTTTTGCTTTTATGCTAACAGGCATTTTCCTGCACTGATaa                                                >  1:924160/1‑149 (MQ=255)
                                                                                              ccTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa                                             >  1:1010631/1‑150 (MQ=255)
                                                                                               cTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAAt                                           >  2:533450/1‑151 (MQ=255)
                                                                                                 acccGCCCCTTTTTTAC‑ATTTGTGT‑AAAGGGGGTT‑TTTTTTTTATGGCGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACg                                              <  1:857023/144‑1 (MQ=255)
                                                                                                 acccGCCCCTTTTTTAC‑AGTTGTGT‑A‑AAGGGGCCTTGTTTTTTATGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa                                             <  2:122527/145‑1 (MQ=255)
                                                                                                 aaCCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCg                                         >  2:340550/1‑151 (MQ=255)
                                                                                                    cgcccctTTTTTAC‑AGTTGT‑A‑AAAAGGGGCCTTGTTTTTTATGACGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATaa                                                <  1:470072/138‑1 (MQ=255)
                                                                                                       ccctTTTTTAC‑GGTTGTGA‑AAAGGGGGCC‑TTTTTTTTATGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCg                                         <  1:397734/145‑1 (MQ=255)
                                                                                                        cccTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGacac                                  >  2:203307/1‑151 (MQ=255)
                                                                                                         ccTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGacac                                  >  1:103667/1‑150 (MQ=255)
                                                                                                              ttAC‑GGTTGTGT‑ACAGGGGGCT‑TTTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAg                             <  2:989796/150‑1 (MQ=255)
                                                                                                              taaC‑ATTTGTGA‑ACAGGGGGGT‑TTTTTTTTTGGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTa                              <  1:551009/147‑1 (MQ=255)
                                                                                                                     ttGTGT‑ACAGGGGCCC‑TATTTTTTAGGGGGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCa                        <  1:690565/149‑1 (MQ=255)
                                                                                                                     ttGTGT‑AAAGGGGGCT‑TGTTTTTTATGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAg                       <  2:359776/150‑1 (MQ=255)
                                                                                                                      tgtgT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtt                     >  1:48535/1‑151 (MQ=255)
                                                                                                                      tggtaaaCAGGGGGCC‑TGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCa                        <  1:170442/143‑1 (MQ=255)
                                                                                                                            aaagggggCT‑TGTTTTTTAGGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct                 <  2:977240/146‑1 (MQ=255)
                                                                                                                             cAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtctc                >  2:181182/1‑150 (MQ=255)
                                                                                                                                tggacc‑TGTTTTTTATAGGGGCTAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACATTAGCATCAGtttt                   <  1:796616/140‑1 (MQ=255)
                                                                                                                                 gggCT‑TGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatga           <  2:1044244/151‑1 (MQ=255)
                                                                                                                                   gCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaat         >  2:765196/1‑151 (MQ=255)
                                                                                                                                       tGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt      >  2:1136171/1‑151 (MQ=255)
                                                                                                                                          tttttAGGACGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTaa    <  1:1163109/150‑1 (MQ=255)
                                                                                                                                            tttATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAAc  <  2:630437/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                          
TTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTAC‑AGTTGTGT‑ACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAAC  >  NC_000913/1622668‑1622954

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