Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,622,816	0	G	A	43.1%	38.4 / 57.7	58	intergenic (+75/+144)	yneM/dgcZ	inner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base G (30/3);  new base A (0/25);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.88e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACT  >  NC_000913/1622666‑1622966                                                                                                                                                            |                                                                                                                                                        ttttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCg                                                                                                                                                         >  2:663309/1‑151 (MQ=255)  tttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCg                                                                                                                                                         >  1:861277/1‑150 (MQ=255)   ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGa                                                                                                                                                        <  2:705618/150‑1 (MQ=255)   ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGa                                                                                                                                                        <  2:635570/150‑1 (MQ=255)           gTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGaaaaaaaa                                                                                                                                                 >  1:808311/1‑149 (MQ=255)                   aataatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTa                                                                                                                                       >  1:860569/1‑150 (MQ=255)                      aatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAAcc                                                                                                                                   >  2:894179/1‑151 (MQ=255)                      aatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAAcc                                                                                                                                   >  1:266827/1‑151 (MQ=255)                          ttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACcggcg                                                                                                                               >  1:244953/1‑151 (MQ=255)                             ttttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTg                                                                                                                            >  2:395000/1‑151 (MQ=255)                                tctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGaa                                                                                                                          >  1:393344/1‑150 (MQ=255)                                tctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGaa                                                                                                                          >  2:792776/1‑150 (MQ=255)                                  tgttttcctgcCCGCTTTTTCCACCAAAAAATGGGATAACAAATGAGGGG‑AAATACCCCCCCCCCAACCCGCCCCCTTGTTACTGTTGTGA‑AAAGGGGGC‑T‑TTT‑TTTTTAGGACGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGaa                                                                                                                          <  2:574654/145‑1 (MQ=255)                                    gttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGc                                                                                                                       >  1:238330/1‑149 (MQ=255)                                      tttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                                   >  1:370695/1‑151 (MQ=255)                                       ttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCtt                                                                                                                     >  2:737365/1‑148 (MQ=255)                                         ctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATg                                                                                                                >  1:40588/1‑151 (MQ=255)                                           ggcccgttATTTCACCCAAAAAGGGGATGACTAATAAGCGG‑AAATATCCCCTCCCCTACCCGGCCCTTTTTTAAATTTTTGT‑AC‑AGGGGC‑C‑TTATTTTTTAGGAGGGCAAAAAAAAT‑CCGCCAGTAAACCGGCGGTGAATGCTTGCAt                                                                                                                 <  1:200200/141‑1 (MQ=255)                                           gcccGGGTATTTCACCAAAAAGTGGGATGACTAAAGAAGGG‑AAATACCCCCCCCCCAACCCGCCCCTTTTTTAGGGTTTTGT‑AAAGGGGGC‑T‑TTT‑TTTTTAGGACGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCAt                                                                                                                 <  1:175677/146‑1 (MQ=255)                                             ccGCTTATTCCAGCAAAAAAGGGGATAACTAAAGAACGG‑AAATAACCCCTCCCCAACCCGCCCCCTTTTTACAGTTGTTA‑ACAAGGGGC‑T‑TTT‑TTTTTATGCCGACAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATgg                                                                                                               <  2:302596/148‑1 (MQ=255)                                             ccGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                             >  2:337241/1‑150 (MQ=255)                                              cgcgTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATgg                                                                                                               >  2:19198/1‑147 (MQ=255)                                                   aTTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAttt                                                                                                       >  1:770315/1‑150 (MQ=255)                                                    ttccacccACAAAGGGGAGGAAAAATGAAGGG‑AAATACCCCCCCACCAACCCGGCCCTTTGTTACAGTTGTGA‑AAAGGGGGC‑T‑TGT‑TTTTTAGGAGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGa                                                                                                          <  1:354776/140‑1 (MQ=255)                                                    ttccaccaaaaagggggATGACAAAT‑AACGGAAGAAAACCCCCCCCCAACCCGCCCCTTTTTTAAATTTTTGA‑AAAGGGGGC‑G‑TGT‑TTTTAAGGCCGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAg                                                                                                           <  1:329701/131‑1 (MQ=255)                                                     ttCACCCAAAAGTGGAATGACTAATGAACGG‑AGATAACCCCTCACCAAACCGCCCCTTTGTTACAGTTGTGA‑AAAAGGGGC‑T‑TGA‑TTTTTATGAGGAAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt                                                                                                        <  2:364750/147‑1 (MQ=255)                                                            acaAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCt                                                                                             >  2:49151/1‑151 (MQ=255)                                                            aaaaaGGGGATGAATAATGACGGG‑AAAAAACCCCCCCCCTACCCGCCCCTTTTTTACAGTTGTTA‑AAAGGGGGC‑T‑TGT‑TTTTTAGGAGGGCAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGtttt                                                                                                <  1:117012/146‑1 (MQ=255)                                                              aaaGGGGATAAATAATAAAGGG‑AAATAATCCCTC‑CCCACCCCGCCCCTTTTTA‑AGTTTTGTAAAAAGGGGC‑CTTTT‑TTTTTATGAGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCttt                                                                                           <  2:625896/151‑1 (MQ=255)                                                               aaTGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt                                                                                          >  1:743966/1‑151 (MQ=255)                                                                aTGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt                                                                                          >  1:358276/1‑150 (MQ=255)                                                                     aTGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCt                                                                                     >  2:66962/1‑150 (MQ=255)                                                                         cTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                                  >  1:311203/1‑149 (MQ=255)                                                                           aaagAACGG‑AAATACCCCCCCCCCAACCCGCCCCTTTTTTACTTTTTTGT‑AAAGGGGGC‑T‑TTT‑TTTTTAGGGGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGc                                                                              <  2:847847/148‑1 (MQ=255)                                                                               aaCGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGGGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAttt                                                                          >  2:244707/1‑151 (MQ=255)                                                                                   gaaaaaaacccccCACCAACCCGCCCCTTTTTTACTGTTGTAA‑AAAGGGGGC‑T‑TTT‑TTTTTATGGGGGGAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTcc                                                                       <  2:671051/142‑1 (MQ=255)                                                                                         aaTCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACt                                                                 >  2:100003/1‑151 (MQ=255)                                                                                          aTCCGTCAACAAACCGGCCACTTGTTATCTTCTGGA‑ACAAGGGGT‑T‑TG‑‑TTTTTAATAGTGAGAAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATGTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCAc                                                                  <  2:8102/149‑1 (MQ=255)                                                                                            ccccccccaacccGCCCCTTTTTTACATTTGTGT‑ACAGGGGGC‑C‑TTT‑TTTTTAGGAGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                              <  1:248321/140‑1 (MQ=255)                                                                                            cccTCCCCAACCCGCCCCTTTTTTACAGTGGTTT‑A‑AAGGGGCTT‑TTT‑TTTTTGTGAGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                              <  2:296511/151‑1 (MQ=255)                                                                                                  ccaacccGCCCCTTTTTTACGGTTGTGT‑AGAGGGGTC‑T‑TGT‑TTTTTAGGAGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAAc                                                           <  2:463379/143‑1 (MQ=255)                                                                                                       ccGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGtt                                                   >  2:1038199/1‑151 (MQ=255)                                                                                                        cGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCg                                                     >  1:548724/1‑148 (MQ=255)                                                                                                        cGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg                                                  >  1:134125/1‑151 (MQ=255)                                                                                                           cccttttttagtgTTGTGT‑ACAGGGGGC‑T‑TGT‑TTTTTAGGAGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGaca                                               <  1:163870/139‑1 (MQ=255)                                                                                                             ccTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGacac                                              >  1:832574/1‑150 (MQ=255)                                                                                                               ttGTTATAGTGGTGT‑ACAGGGGGC‑T‑TTT‑TTTTTAGGACGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGacac                                              <  2:319043/148‑1 (MQ=255)                                                                                                                        ttGTTT‑AAAAGGGGC‑C‑TGA‑TTTTTAGGCGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCa                                    <  1:121128/149‑1 (MQ=255)                                                                                                                        ttGTGT‑CAAGGGGGT‑T‑TGT‑TTTTTAGGCCGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc                                     <  1:526467/148‑1 (MQ=255)                                                                                                                           gtgaaaggggggC‑C‑TTT‑TTTTTATGCGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAg                                   <  2:469964/140‑1 (MQ=255)                                                                                                                                cAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtctc                            >  1:829775/1‑150 (MQ=255)                                                                                                                                 aaGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa                           >  2:419959/1‑150 (MQ=255)                                                                                                                                      ggggctttttttttATGACGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaat                         <  2:860441/141‑1 (MQ=255)                                                                                                                                      ggggctttttttttAGGAGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa                           <  2:128930/139‑1 (MQ=255)                                                                                                                                           tttttttttAGGCCGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaat                     <  2:861149/145‑1 (MQ=255)                                                                                                                                               tttttATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTaa                >  1:645965/1‑150 (MQ=255)                                                                                                                                               tttttATGACGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt                  <  2:417943/148‑1 (MQ=255)                                                                                                                                               ttttaATGAGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTaa                <  1:624207/150‑1 (MQ=255)                                                                                                                                                      gACGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAg          <  2:249694/149‑1 (MQ=255)                                                                                                                                                           cGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAAc   >  1:662723/1‑151 (MQ=255)                                                                                                                                                            gAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACt  >  2:876009/1‑151 (MQ=255)                                                                                                                                                            |                                                                                                                                                        TTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGG‑AGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGT‑ACAAGGGGC‑C‑TGA‑TTTTTATGACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACT  >  NC_000913/1622666‑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: ‑