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.9%	46.3 / 68.6	66	intergenic (+75/+144)	yneM/dgcZ	inner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base G (35/2);  new base A (0/29);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.04e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCC‑GCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTA  >  NC_000913/1622675‑1622962                                                                                                                                                    |                                                                                                                                                    cGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCC‑GCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑c                                                                                                                                           <  2:284332/151‑1 (MQ=255)  gTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCC‑GCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑cc                                                                                                                                          >  2:248826/1‑151 (MQ=255)   tACTGGGAATAATTTTATTTTCTGGTTTTCTGGCC‑GCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCg                                                                                                                                         >  1:207945/1‑151 (MQ=255)     cgtggAATAATTTTATTTTCTGGTTTTCTGGGC‑GCGTATCTCCGCCAAAAATTGGACGGCTAATGAGCGGAGACACGCCCTCACCGAGCCGGGCCC‑TTGAT‑ACAGTTG‑GCTTCAAGGGGCC‑TGAGTTATATG‑ACGTCAAAAAAAAA‑CCGc                                                                                                                                        <  2:606113/147‑1 (MQ=255)     cTGGGAATAATTTTATTTTCTGGTTTTCTGGCC‑GCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGcc                                                                                                                                       >  2:90793/1‑151 (MQ=255)           ataatTTTATTTTCTGGTTTTCTGGCC‑GCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTaa                                                                                                                                  >  1:617575/1‑150 (MQ=255)            taatTTTATTTTCTGGTTTTCTGGCC‑GCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAAc                                                                                                                                >  1:370053/1‑151 (MQ=255)                ttttttttttGGTTTTTCGGCCGGTGTTTTTCACCCAAAAAGGGAATGACAAAGGAAGGGAGAT‑ATCCCCCCCCCAACCCGCCCC‑TTTTT‑AAAGTTGTTGTGAAAGGGCGC‑TGTTTTTTAGG‑AGGGGAAAAAAAAA‑CCGCCAGTAAACcg                                                                                                                              <  1:192801/150‑1 (MQ=255)                 ttattttCTGGTTTTCTGGCC‑GCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACcggcg                                                                                                                           >  2:896335/1‑151 (MQ=255)                      ttctggttttctggCC‑GCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGa                                                                                                                       >  1:1053314/1‑150 (MQ=255)                          ggttttctggCC‑GCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCt                                                                                                                  >  1:376263/1‑151 (MQ=255)                          ggttttctggCC‑GCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCt                                                                                                                  >  2:601361/1‑151 (MQ=255)                            ttttctggCC‑GGTTTTTCCAGCCAAAAATGGGATGATTAAGGAAGGGAGATAACCCCCCCCCAACCCGGCCCT‑TTATT‑AAATTTG‑TTTAAAGGGGGCT‑TGATTTTTATG‑ACGGCAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATg                                                                                                                    <  2:810603/147‑1 (MQ=255)                            ttttctggCC‑GCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCtt                                                                                                                 >  1:410669/1‑150 (MQ=255)                            ttttcgggCC‑GTTTTTTTCACCAAAAAAGGGGATGATTAAGGAAGGGAAATATCCCCCCCCCTACCCGCCCCT‑TTTTT‑ACATTTG‑TGTA‑AAGGGGCCTTTTTTTTTATG‑CGGGCAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCt                                                                                                                  <  1:249168/149‑1 (MQ=255)                             tttctggCC‑GCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTg                                                                                                                >  1:1153760/1‑150 (MQ=255)                             tttcggcCC‑GTGTTTTTCACCCAAAAAGGGGATGACTAATGAAGGGAGATAACCCCCCCCCAACCCGCCCCT‑TTTTT‑ACATTTG‑TGAAAAGGGGGCT‑TGTTTTTTATG‑CCGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTg                                                                                                                <  1:67208/150‑1 (MQ=255)                              ttctggCC‑GCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                               >  1:172103/1‑150 (MQ=255)                                 tggCC‑GCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATgg                                                                                                           >  1:744202/1‑151 (MQ=255)                                          tATTTCAGCAAAAAAGGGGATAACTAATGAAGGGAAATAACCCCTCCCCTACCCGCCCCT‑TTGTT‑ACAGTTG‑TG‑AAAAGGGGCCTTGTTTTTTATG‑AGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATa                                                                                                        <  1:16957/146‑1 (MQ=255)                                            ttccaccaaaaagggggATGACAAATAAAGGGGAATACTCCCCCCCCCCCCCGCCCCT‑TTGTT‑AGATTTG‑TGAACAGGGCCGT‑TGATTTTGATG‑GCGTAGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt                                                                                                    <  1:1119461/134‑1 (MQ=255)                                                  ccACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAAACGACCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGCTTTTTATG‑ACGGCGAAACAAAA‑CCGCCAGCAAACCGGCGGTGAATACGTGGATGGCTTGGTTTttgt                                                                                               >  2:445660/1‑147 (MQ=255)                                                      aaaGGGGATGATAAATGAAGGGAAAAAACCCCCCCCCTACCCGCCCCT‑TTTTA‑ACGTTTG‑TGTACAGGGGGCT‑TTTTTTTTAGG‑CGGGAAAAAAAAAC‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGttt                                                                                             <  1:1121415/145‑1 (MQ=255)                                                             aTGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCt                                                                                 >  1:300564/1‑150 (MQ=255)                                                              tGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTa                                                                                >  1:275372/1‑150 (MQ=255)                                                               gACTAATGAACGGAGATAATCCCTCACCTAATCGGCCCC‑TTGTT‑ACAGTTA‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CAGCCAGTAAACCAGCGGCGCATGCATGCACGGGTAGGTTTGGTTTTTGCTTTTATGCt                                                                                 >  2:965641/1‑148 (MQ=255)                                                                aCTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                              >  1:142/1‑150 (MQ=255)                                                                 cTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTaa                                                                               >  2:841039/1‑148 (MQ=255)                                                                    aTGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGc                                                                          >  1:333006/1‑150 (MQ=255)                                                                     tGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAt                                                                        >  2:297962/1‑151 (MQ=255)                                                                       aaCGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGTATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAttt                                                                      >  2:788300/1‑151 (MQ=255)                                                                       aaCGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAttt                                                                      >  2:17493/1‑151 (MQ=255)                                                                         aggaataactCCCCCCCCTACCCGCCCTT‑TTTTT‑ACGGTTT‑TGTATAAGGGCCC‑TGTTTTATATG‑TCGGAAAAAAAAAT‑CCGCCAGTAAACCGGCGGTTAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTc                                                                    <  2:1144464/142‑1 (MQ=255)                                                                          gaaaataccccccCACCTACCCGCCCCC‑TTTTA‑ACAGTTG‑TG‑AAAAGGGGCC‑TTGTTTTTAAG‑GGGGGAAAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtttt                                                                     <  2:294162/140‑1 (MQ=255)                                                                           ggaaaaaacccctcccTAACCCGCCCCTTTTTT‑ACGTTTG‑TGAAAAGGGGCCC‑TTTTTTTTAGG‑ACGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCa                                                                         <  1:19353/131‑1 (MQ=255)                                                                                     ccccccccaCCCGCCCC‑TTTTTAACGTTTG‑TGAAAAGGGGGCG‑TGTTTTTTAGG‑AGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATa                                                         <  2:380924/143‑1 (MQ=255)                                                                                        cccaccccGCCCTC‑TTGTT‑AGATTTG‑TGAAAGGGGGGCT‑TTTTTTTTAGG‑AGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATaa                                                        <  1:588625/142‑1 (MQ=255)                                                                                          cTACCCGGCCCT‑TTGTT‑AAATTTG‑TGTACAGGGGGTT‑TGTTTTTTAGG‑CGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAAc                                                       <  1:561442/147‑1 (MQ=255)                                                                                          cTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGaa                                                    >  2:928974/1‑150 (MQ=255)                                                                                            acccGGCCCT‑TTTTT‑ACAGTTT‑TGAAAAGGGGGCT‑TTATTTTTATG‑ACGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATc                                                  <  1:371327/148‑1 (MQ=255)                                                                                             cccGCCCCT‑TTGTT‑AGTTTTG‑TGTACAGGGGGCT‑TTTTTTTTATG‑ACGAAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCg                                                 <  2:615688/149‑1 (MQ=255)                                                                                               cGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg                                              >  1:241042/1‑151 (MQ=255)                                                                                                  ccct‑TTTTT‑ACATTTG‑TGAAAAGGGGGCT‑TTTTTTTTAGG‑AGGGCAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGaca                                           <  1:90816/151‑1 (MQ=255)                                                                                                  cccc‑TTTTT‑AATGTTG‑GGGAAAAGGGGGC‑TGTTTTTTAAGCAGGAAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGaa                                                    <  2:981859/143‑1 (MQ=255)                                                                                                     c‑TTGTT‑ACCGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGacac                                          >  2:131447/1‑149 (MQ=255)                                                                                                       ttGTT‑AAAGTTG‑TGAAAAGGGGGCT‑TTTTTTTTATG‑CGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGAcaca                                         <  1:681157/149‑1 (MQ=255)                                                                                                         ttgttaggtttttgtACATGGGTCT‑TTTTTTTTAGG‑ACGGAAAAAAAAAT‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGAcaca                                         <  2:20102/137‑1 (MQ=255)                                                                                                             acttttG‑TGTAAAGGGGGCT‑TGTTTTTTATG‑ACGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCa                                   <  1:765529/146‑1 (MQ=255)                                                                                                                gTTG‑TGTAAAGGGGGCT‑TGTTTTTAATG‑AGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc                                 <  1:147095/149‑1 (MQ=255)                                                                                                                   tgtgaaaaggggatc‑TGATTTTAATG‑AGGGGAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGt                              <  2:340207/136‑1 (MQ=255)                                                                                                                   g‑tgtACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGttt                            >  2:593703/1‑151 (MQ=255)                                                                                                                     tgtACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtttt                           >  2:487741/1‑151 (MQ=255)                                                                                                                     tgtACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtttt                           >  1:317604/1‑151 (MQ=255)                                                                                                                      gtACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTAAGTTAAAAGGGATTTTTCTGCACTGATAAGGAATCGTTGAAACAGTAGAATAAGttttt                          >  2:414070/1‑150 (MQ=255)                                                                                                                       aaaagggggCT‑TGTTTTTTAGG‑AGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtc                          <  1:240778/145‑1 (MQ=255)                                                                                                                        aCAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtctc                        >  1:390700/1‑151 (MQ=255)                                                                                                                         aaaGGGGCCATGTTTTTAATG‑GGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtttt                           <  1:542291/147‑1 (MQ=255)                                                                                                                             ggggctttttttttATG‑CCGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatg                    <  1:489839/142‑1 (MQ=255)                                                                                                                             ggggctttttttttAGG‑ACGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaa                      <  2:953364/140‑1 (MQ=255)                                                                                                                              gtcc‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatga                   >  2:631012/3‑150 (MQ=255)                                                                                                                              gggctttttttttAGG‑ACGAAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatga                   <  1:630966/143‑1 (MQ=255)                                                                                                                                   ttaTTTTTATG‑GCGAAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaa                  <  2:1030058/145‑1 (MQ=255)                                                                                                                                      tttttATG‑GGGGAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaatg                <  1:576799/146‑1 (MQ=255)                                                                                                                                      tttttATG‑ACGGCG‑AAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAAc          >  1:153771/1‑151 (MQ=255)                                                                                                                                                 ggCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTa  >  1:833415/1‑150 (MQ=255)                                                                                                                                                 ggCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTa  >  1:1005157/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCC‑GCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTG‑TGTACAAGGGGCC‑TGATTTTTATG‑ACGGCGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTA  >  NC_000913/1622675‑1622962

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