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

TTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGAT‑GACTAATGAACGGAGATAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACT  >  NC_000913/1622667‑1622966
                                                                                                                                                           |                                                                                                                                                       
tttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGAT‑GACTAATGAACGGAGATAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCg                                                                                                                                                         >  2:550860/1‑150 (MQ=255)
 ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGAT‑GACTAATGAACGGAGATAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGa                                                                                                                                                        >  1:528386/1‑150 (MQ=255)
    tGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGAT‑GACTAATGAACGGAGATAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGaaaa                                                                                                                                                     <  1:744004/150‑1 (MQ=255)
     ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGAT‑GACTAATGAACGGAGATAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGaaaaa                                                                                                                                                    >  2:896864/1‑150 (MQ=255)
      gCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGAT‑GACTAATGAACGGAGATAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGaaaaaaa                                                                                                                                                  <  2:638940/151‑1 (MQ=255)
          tACTGGGAATAATTTTATTTTCTGGTTTTCTGGCTGCGTATTTCAGCCACAAATGGGAT‑GACTAATGAACGGAGATAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCg                                                                                                                                             >  1:352427/1‑151 (MQ=255)
              gggAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGAT‑GACTAATGAACGGAGATAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGcc                                                                                                                                           <  2:998802/149‑1 (MQ=255)
              gggAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGAT‑GACTAATGAACGGAGATAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCa                                                                                                                                          <  1:820836/150‑1 (MQ=255)
                 aataatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGAT‑GACTAATGAACGGAGATAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTa                                                                                                                                       >  1:129386/1‑150 (MQ=255)
                    aatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGAT‑GACTAATGAACGGAGATAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑CAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAAcc                                                                                                                                   >  1:804672/1‑151 (MQ=255)
                    aatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGAT‑GACTAATGAACGGAGATAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAAcc                                                                                                                                   >  1:106292/1‑151 (MQ=255)
                      ttttattttCGGTTTTTCTGGCTGCTTTTTGCACCCA‑AAATTGGATGGATAAATGAAGGGA‑AAAATCCCCCCCCCAACCCGCCCCCTTGTAAAAT‑TTGTGTCA‑AGGGGGCT‑TTT‑TTTTTATGCCGGAAAAAAAAA‑ACCGCCAGTAAAcc                                                                                                                                   <  1:896616/149‑1 (MQ=255)
                           ttttctggttttctggCCGCGTATTTCAGCCACAAATGGGAT‑GACTAATGAACGGAGATAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTg                                                                                                                            >  1:519475/1‑151 (MQ=255)
                            tttctggttttctggCCGCGTATTTCAGCCACAAATGGGAT‑GACTAATGAACGGAGATAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTg                                                                                                                            >  1:494224/1‑150 (MQ=255)
                                 ggttttctggCCGCGTATTTCAGCCACAAATGGGAT‑GACTAATGAACGGAGATAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCt                                                                                                                      >  2:861297/1‑151 (MQ=255)
                                       cggcccGCGTTTTTCAGCCAAAAAGGGGAT‑GACAAATAAACGGGAATAACCCCCCCCCAAC‑CCGCCCCTTTTTTACAG‑TTGTGTAC‑AGGGGGCT‑TGA‑TTTTTAAGAGGGAAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCAt                                                                                                                 <  2:228093/146‑1 (MQ=255)
                                         gccggcttATTTCACCCAAAAAGGGGAT‑AACAAATAAGGGAAAATAACCCCTC‑CCCAA‑CCCGCCCCTTTTTAAAG‑TTTTTGAA‑AAGGGGCCTTGT‑TTTTTAGGCCGGAAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCAt                                                                                                                 <  2:152925/141‑1 (MQ=255)
                                              cgTATTTCAGCCACAAATGGGAT‑GACTAATGAACGGAGATAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                             >  2:937234/1‑147 (MQ=255)
                                                tttttCACCCAAAATTGGGAG‑GCTTAATG‑ACGGGGGAAATCCCCCCACCAA‑CCCCCCCTTTTTTTCAGTTTGTGTAG‑AGGGGGCT‑TGT‑TTTTGATGAGGGGAAAAAAAC‑CCCGCCAGTAAACCGGCGGTGAATGCTTGCa                                                                                                                  <  1:356841/138‑1 (MQ=255)
                                                        ccACAAATGGGAT‑GACTAAAGAGGGGAGATATCCCCTCCCCTAC‑CCGGCCCTTTTTTACGT‑TTGTGTAC‑AGGGGGCT‑TGA‑TTTTTAGGAGGGAAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTg                                                                                               <  2:594783/151‑1 (MQ=255)
                                                          aaaaaGGGAAT‑GATAAATGAAGGGAAAAAACCCCTCCCCTAC‑CCGCCCCTTTTTTCCAG‑TTGTGTAC‑AGGGGCCT‑TGT‑TTTTTAGGACGGCAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGc                                                                                              <  2:528397/148‑1 (MQ=255)
                                                            aaaTGGGAT‑GACTAATGAACGGAGATAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCCTGCATGGATAGATTTGTGTTTTGCttt                                                                                           >  1:61376/1‑151 (MQ=255)
                                                             aaaggaggtgataaaTAAACGGAAATACCCCCTCACCAAC‑CCGCCCCTTTGTTA‑GG‑TT‑TGTACAAAGGGGGC‑TTT‑TTTTTAAGACGGGAAAAAAAACCCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCttt                                                                                           <  1:77860/138‑1 (MQ=255)
                                                                gggAT‑GAAAAAGGAACGGAAATAACCCCTCCCCTAC‑CCGCCCCTTTTTTACGT‑TTGTGTCC‑AGGGGCCT‑TTT‑TTTTTATGAGGGAAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACg                                                                                       <  1:79001/151‑1 (MQ=255)
                                                                       aaaaaTGAAGGGAAATAACCCCCCCCCAAC‑CCGCCCCTTTGTTAAAG‑TGGTTAAA‑AGGGGGCT‑TTT‑TTTTTAGGACGGAAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACa                                                                                 <  1:990785/148‑1 (MQ=255)
                                                                          aaTGAAGGGAGATACCCCCTCACCTAC‑CCGGCCCTTTTTTACGG‑TGGTGTAG‑AGGGGCCT‑TTT‑TTTTTATGACGGAAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAgg                                                                               <  2:790438/150‑1 (MQ=255)
                                                                                     aTAATCCCCCCCCTAC‑CCGGCCCCTTTTTACGG‑TTGTGAAC‑AAGGGGCC‑TGT‑TTTTTATGACGGAAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCt                                                                      <  1:721707/148‑1 (MQ=255)
                                                                                      tAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCAc                                                                  >  1:864308/1‑151 (MQ=255)
                                                                                          cccTCCCCTAC‑CCGCCCCTTTTTTACAG‑TTGTGAAC‑AGGGGGCT‑TTT‑TTTTTAGGACGGAAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTg                                                                <  2:575798/149‑1 (MQ=255)
                                                                                           ccTCCCCTAC‑CCGCCCCTTTTTTACGT‑TTGTGAAA‑AGGGGGCT‑TGT‑TTTTTATGACGGAAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATa                                                             <  2:61961/151‑1 (MQ=255)
                                                                                               cccccacccGCCCCCTTTTTATAG‑TTGTGAAC‑AGGGGGCC‑TTT‑TTTTTATGAGGGAAAAAAAAC‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATa                                                             <  2:201329/141‑1 (MQ=255)
                                                                                                 ccaacccGCCCCTTTTTTACTG‑TTGTTTAA‑AGGGGGCT‑TTT‑TTTTTATGGCGGAAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAAc                                                           <  1:549745/143‑1 (MQ=255)
                                                                                                 cTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGaa                                                        >  1:470085/1‑150 (MQ=255)
                                                                                                 cTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAAt                                                       >  1:139548/1‑151 (MQ=255)
                                                                                                 cTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAAt                                                       >  2:667369/1‑151 (MQ=255)
                                                                                                  tAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATc                                                      >  2:137639/1‑151 (MQ=255)
                                                                                                       cGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGtt                                                   >  1:1110119/1‑150 (MQ=255)
                                                                                                          ccctTTTTTACAT‑TTGTG‑AA‑AAGGGGCCTTGT‑TTTTTATGAGGGAAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg                                                  <  1:263280/148‑1 (MQ=255)
                                                                                                                 ttACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAg                                         >  1:261817/1‑150 (MQ=255)
                                                                                                                 ttACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGc                                        >  1:1115384/1‑151 (MQ=255)
                                                                                                                  tACAG‑TTGTGTAC‑AGGGGTCT‑TTT‑TTTTTAGGACGGAAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTa                                          <  1:82925/148‑1 (MQ=255)
                                                                                                                  tACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGc                                        >  2:805900/1‑150 (MQ=255)
                                                                                                                    cAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc                                     >  2:328662/1‑151 (MQ=255)
                                                                                                                             tgaaaaGGGGCCTTGT‑TTTTTATGGCGGCAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGttt                                <  2:1036984/146‑1 (MQ=255)
                                                                                                                             tataaaGGGGCC‑TTTTTTTTTAGGAAGGAAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGt                                  <  1:133172/144‑1 (MQ=255)
                                                                                                                             tAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct                             >  1:46318/1‑151 (MQ=255)
                                                                                                                             tAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct                             >  2:6210/1‑151 (MQ=255)
                                                                                                                             tAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct                             >  2:17617/1‑151 (MQ=255)
                                                                                                                                 aaGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa                           >  2:602803/1‑150 (MQ=255)
                                                                                                                                 aaGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa                           >  2:1060682/1‑150 (MQ=255)
                                                                                                                                  aGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaat                         >  1:967018/1‑151 (MQ=255)
                                                                                                                                   ggggCT‑TGT‑TTTTTAGGACGGAAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaa                          <  2:90016/149‑1 (MQ=255)
                                                                                                                                   agggggttttttttttATGAGGGAAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtttt                               <  1:668206/135‑1 (MQ=255)
                                                                                                                                    ggggtgtttttttttACGACGGAAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtctc                            <  1:223393/138‑1 (MQ=255)
                                                                                                                                    gggggctttttttttATGACGGAAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa                           <  1:550849/139‑1 (MQ=255)
                                                                                                                                     ggggctttttttttAGGGCGGAAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa                           <  1:785137/139‑1 (MQ=255)
                                                                                                                                      gCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaat                     >  2:142853/1‑151 (MQ=255)
                                                                                                                                        gctttttttttATGCGGGAAAAAAAAA‑CCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaa                      <  1:114914/144‑1 (MQ=255)
                                                                                                                                          tGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt                  >  1:459863/1‑151 (MQ=255)
                                                                                                                                              tttttAGGACGGCAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt                  <  2:152983/148‑1 (MQ=255)
                                                                                                                                                        ggCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTa      >  2:1039042/1‑150 (MQ=255)
                                                                                                                                                        ggCAAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTa      <  1:248787/150‑1 (MQ=255)
                                                                                                                                                           gAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACt  >  2:756447/1‑151 (MQ=255)
                                                                                                                                                           |                                                                                                                                                       
TTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGAT‑GACTAATGAACGGAGATAATCCCTCACCTAA‑CCGGCCCCTTGTTACAG‑TTGTGTAC‑AAGGGGCC‑TGA‑TTTTTATGACGGCGAAAAAAA‑ACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACT  >  NC_000913/1622667‑1622966

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