Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,825,061	0	G	A	35.0%	57.6 / 19.2	40	intergenic (‑124/+79)	ves/spy	HutD family protein Ves/ATP‑independent periplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base G (24/2);  new base A (0/14);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 5.17e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCC  >  NC_000913/1824817‑1825302                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   cAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTAGCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAAcc                                                                                                                                                                                                                                              <  2:202278/251‑1 (MQ=255)     gCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑CTTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGGGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACACCT‑GTTTGAGAA‑GAGAAAAGAAAAccgccg                                                                                                                                                                                                                                          >  2:534361/1‑251 (MQ=255)                                ggcgggTGGAAACGTCAAAATTTCGCGCTTTTCCCCGGCATCATTCCCCCAAAGATTCCCCG‑TATTTTTCCCAATATCAAAGTTTTCCTTCCCCCCCCCCGGCCTTTTTTTTTC‑TTTGTTATCGTTCAGGTACCCG‑TGG‑TTGTTTATCGCCAACG‑CGGGTGATTTACGTT‑AGTGGTGATCAGTTTTTTTTTACACTT‑TTTTGAGAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGt                                                                                                                                                                                                               <  1:167584/247‑1 (MQ=255)                                       gggAACGCGCAAATTTCGCGTTTTCCCCCGGCAGTATTCCCCCACAGTTTCCCCGAAATTTTACGAA‑TAACAAAGTTTCCCATCCCACCTCCCGCCCTTTTCTTTTC‑TTTTTTATCGTTCAGGTACCCG‑TGG‑AGGTTTATCGGCAACT‑CGGGTTTTTTACGTT‑AGTGGTGATCAGTTTTTTTTTACAACT‑TTTTGAGAA‑GAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc                                                                                                                                                                                                        <  1:562478/251‑1 (MQ=255)                                                aaaaTTTCGCGCGTTTCACCGGCAGTTTTTCCCAAAATTTTCCCCGACTTTTTACGCA‑TATAAAAGTATCCCATCCTCGCCCCCGCCCCGTTCTTTTA‑TTTGTTATCGTTCAGGTACCCG‑GGG‑TTGTTTAGCGGAACCT‑GGGGTTTTTTACGTT‑AGTGGTGATCGGTTTTTTTTTACACCT‑GTTTGAGAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                                                                                                                                                                                               <  1:395290/250‑1 (MQ=255)                                                 aaaTTTCGCGGTTTTCACCGGAAGCTTTTCCCCACAGATTCACCGACTTTTTACGCA‑TATAAAAGTTTTCCATACCCCCCCCCGCCCCTTTTTTTTC‑TTTTTTATCTTTC‑GGTAACCG‑GGGGAGGTTTATGGGAAACT‑GGGGTTATTTACGTT‑AGGGGTGATCAGGTTTTTTTTACAACT‑GTTTGAAAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                                                                                                                                                              <  1:396121/251‑1 (MQ=255)                                                  aaTTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgg                                                                                                                                                                                             >  1:253905/1‑251 (MQ=255)                                                         cgcgTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                                                                                                                                    >  1:163749/1‑221 (MQ=255)                                                              ttCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGTTCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgggtttcaggtttt                                                                                                                                                                                 >  2:491178/1‑239 (MQ=255)                                                                     gcagctTTCCGCAACAGATCCCCCGACTTTTTACGCA‑TATAAAAGATTTCAATCCCCCCCCCGGGCCTGTTTTTTTA‑ATTTTTAGCGTTCAGGTATCCG‑TGG‑ATGTTTACCGG‑AACTTGGGGTGATTTCCTTT‑AGTGGTGATCAGGTTTTTTTTAAAACT‑GTGTGAGAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGa                                                                                                                                                                         <  1:612000/251‑1 (MQ=255)                                                                      cagcaTTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGa                                                                                                                                                                         >  1:495536/1‑250 (MQ=255)                                                                       agcaTTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGGTTTTTGGCTGaa                                                                                                                                                                        >  2:4772/1‑250 (MQ=255)                                                                         caTTGCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGATAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAAAAGACCGGCGGA‑Cttcaatagttttgctgaa                                                                                                                                                                        >  2:163293/1‑232 (MQ=255)                                                                                   caGATTCACCGACATTTTCCGCA‑TATCACAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAACGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑TTTTAATTTTTTTGCtttaaaatttatctgcc                                                                                                                                                           >  2:462701/1‑235 (MQ=255)                                                                                   caGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTcagcc                                                                                                                                                           >  1:306897/1‑250 (MQ=255)                                                                                            cGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa                                                                                                                                                                                                                        >  1:388664/1‑182 (MQ=255)                                                                                                aTTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGG‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAttt                                                                                                                                              >  1:576540/1‑251 (MQ=255)                                                                                                      cctaataaaaaagtttcccATACCCTCTCCCGGCCCGTTTTTTTC‑TTTGTTATCGTTCAGG‑AACCGGTGT‑AGTGTTATCGGCAACG‑GGGGTGTTTTACGTT‑AGTGGTAATCGGGTTTTTTTAACACCT‑TTTTGAGAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTAcc                                                                                                                                          <  1:390729/233‑1 (MQ=255)                                                                                                       cacataacaagatATTCCACACTCGCTCCCGGCCGGTT‑TTTTCTATTGTTATTGTTAAGGTAACCG‑TGG‑AGGGTTATCGGCAACT‑TGGGTGATTTCCGTT‑AGTGGTAAACGGTCTTTTTTTACACTT‑GTTGGAGAA‑AAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtt                                                                                                                                        <  1:176367/239‑1 (MQ=255)                                                                                                                 aGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTgg                                                                                                                              >  1:701850/1‑251 (MQ=255)                                                                                                                             tCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTTAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGtt                                                                                                                    >  2:358900/1‑248 (MQ=255)                                                                                                                               ccccccGGCCTGTTTTTTTT‑ATTGTTATCGGTCAGGTAGCCG‑GGG‑TGGGTTATCGGAAACG‑CGGGTGTTTTA‑GTTGAGTGGTGATCGGGTTTTTTTAACAACT‑GTTGGAAAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc                                                                                                                <  2:306897/248‑1 (MQ=255)                                                                                                                                     ggCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAAAT‑CGGGTGTTTTACGTT‑AGTGGTGCTCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAACGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGG‑CTTCAGTTTTTTGGCTGAAAGATTATTCAGCAGTTTCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCCAACGCTtct                                                                                                          >  2:318542/1‑251 (MQ=255)                                                                                                                                        ccGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCa                                                                                                        >  1:111507/1‑250 (MQ=255)                                                                                                                                               tttC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGc                                                                                                >  1:94797/1‑251 (MQ=255)                                                                                                                                                tttcatgtttttcgtcaggttatccG‑TGG‑GTGTTTATCGCCACCT‑GGGGTGATTTACTTT‑GGGGGTGATCATGTTTTTTTTAAACCG‑GTTTGAGAA‑AAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAg                                                                                                 <  1:478225/227‑1 (MQ=255)                                                                                                                                                  c‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCAtt                                                                                             >  1:45009/1‑251 (MQ=255)                                                                                                                                                               taaGGTACCCG‑TGG‑AGGTTTATCGGAATCT‑GGGGTGTTTGACGTA‑AGTGGTGATCGGGTTTTTTTTACA‑GT‑GTTTGAAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCtttt                                                                                  <  1:458847/248‑1 (MQ=255)                                                                                                                                                                    gtaacggggg‑ATGGTTATCGGCATGT‑GGGGTGTTTTCGGTT‑GGTGGTGATCAG‑TTTTTTTTGAAACTTGTTGGAGAA‑AAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGt                                                                             <  1:4772/243‑1 (MQ=255)                                                                                                                                                                                  gTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGCA‑GAGAAAAGAAAAACGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTTGCTGAAAGATTATTGAGCAGGTGCAGGGCTTTTTCCTTTTGCCGCTGGGCGGTCTGTCTAACGCTTATCAAAATTAGCATTAAATTGCTTTTTCTTTTTCGGctgcaggg                                                                >  1:301170/1‑244 (MQ=255)                                                                                                                                                                                        ggCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCCAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGAtgttgt                                                          >  1:554632/1‑251 (MQ=255)                                                                                                                                                                                                      aTTTACGTT‑AGAGGTGATCCGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAAACGCCGATCCTGTCCACCGCATTACTGCCAGGTAGTGGACAAGACCGGCGGT‑CTTACGTTTTTTTGCTGAAAGATTCTTCAGCAGTTGCAGGGATTTTATCTTTTGTCGCTGGATGTTTTGTCAGACGCTTCTCAAAATTAGACTTAACTTGCTTTTTCTGTTCCGGCGTCCGGGTGTTGTAGATTTTTTTCTg                                             >  2:629076/1‑251 (MQ=255)                                                                                                                                                                                                         tACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGAGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGAAGGTGCAGGCATTTTACCTTTTGCCGCTGGACGGTCTGTCAGACGATTCTCAAAATTAGCAATAAATTGCTTTTTTTGTTCCGGCGCCAGGGTGTTGTAAATTTTGTTCTGGGt                                          >  2:172322/1‑251 (MQ=255)                                                                                                                                                                                                             tt‑AGTGGTGCTGAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAACCCGACGATCCTGTCCACCGCATTACTGCAAGGTAGAGGACAAGAACGGCGGAACTAAAGTTTTT‑GGGTGAAAGCTTATTCAACACTTGCAGGCCTTTTACCTTTTTCCCCTGGACGATCTGTCCAACGCTTCCCAAAATTAGCATTTAATTTTTTTTTCTGTTTCGGCGTCAGGGTGTTGTAAATTTTGTTCTggggtttc                                      >  2:287800/1‑246 (MQ=255)                                                                                                                                                                                                                                ttttttacactggtgtGAGAA‑AAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATGCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCTGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCTGCATGt                      <  1:450158/235‑1 (MQ=255)                                                                                                                                                                                                                                    ttacattttttagaaaa‑gaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGc                  <  1:659061/235‑1 (MQ=255)                                                                                                                                                                                                                                       aaCT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGATTCTTGGCAGGAAGAATTTTAAGCAGATACAGGCATTTTAACTTTTGCCGCCCGACGTTCTCTCAGACCCTTCTCAACAATACCATTTCATTGATTTTTCTCTTTCGGGGCCAGGATGTTGTAAATTTTCTTTTGGGTTTCCATGTCCGCCAgacagttagcgtttc             >  1:671777/1‑237 (MQ=255)                                                                                                                                                                                                                                         actttttGAGAA‑GAAAAAAAAAAACCGCCGTTCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTgcg            <  1:491178/247‑1 (MQ=255)                                                                                                                                                                                                                                                 agaa‑gaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGttctt    <  2:253905/251‑1 (MQ=255)                                                                                                                                                                                                                                                   aa‑gagaaAAGAAAACCGCCGATCCTGTCCACCGAATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAAACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTcc  >  1:197267/1‑251 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   CAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCC  >  NC_000913/1824817‑1825302

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