Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,978,164	T→A	21.8%	Y245F (TAT→TTT)	lysA ←	diaminopimelate decarboxylase, PLP‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,978,164	0	T	A	21.8%	104.6 / 6.6	58	Y245F (TAT→TTT)	lysA	diaminopimelate decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base T (24/19);  new base A (3/10);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 5.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGGTGTACC  >  NC_000913/2977930‑2978398                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           tCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCataata                                                                                                                                                                                                                                        <  1:739386/239‑1 (MQ=255)     cTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATAAGCCAGACTGCGCTACCC‑GGG‑AGCGACCCGGTTCAATTT‑CCCGTTTCACAGGGTGGCCCAAA‑A‑AGCGGGAGGATTGTTCACCCCCGGC‑GT‑TTCACAGA‑CCCtaacaaatt                                                                                                                                                                                                                                    >  1:661931/1‑236 (MQ=255)      tGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGtt                                                                                                                                                                                                                                   <  1:179834/239‑1 (MQ=255)         cAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCgg                                                                                                                                                                                                                                <  1:1011610/239‑1 (MQ=255)                 ttatgtggtgtTACCTCCCGTTCTTTGCGGGGCGCATCAGATTTTGGAACCCGGCATAACCCAGAACAAAGGGGGGG‑CCCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGGGCTACCA‑GGA‑AGCGCCCCGGTTCAATTT‑CCAGTTTCACAGGGGGGCCCAAA‑T‑GGCGGGCGTTTTGCTCCCGCGCGGC‑AT‑CCCAAAGA‑CCATAATAAGGTGCGGTATCAAcc                                                                                                                                                                                                                        <  2:376965/234‑1 (MQ=255)                        gtggtAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCC‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTTACAAGGGGGCCCACA‑T‑TGCGGGCGATTTGCTCACCCGCGGC‑CT‑TTCCCCAA‑CCATAATAATGTTCGGGATCCACccccctt                                                                                                                                                                                                                  >  1:172804/1‑235 (MQ=255)                             aaCTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑TGCGGGCGAATTTCTCACCCGCGGC‑AT‑TCCCCCAA‑CCAACATAATATTCGGGATCCAACCCCTCtttcac                                                                                                                                                                                                             >  1:315889/1‑234 (MQ=255)                                   cGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑CGCGACCCCGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTTCTCCCCCGCGGC‑CT‑TTCACAGA‑CCCTAATAATTTTCGGGATCACACGCCCCTTCAccccctttg                                                                                                                                                                                                      >  1:187825/1‑233 (MQ=255)                                     tACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGGGGCCCCAA‑T‑GGCGGGCGAGTTTCTCACGCGCGGC‑AT‑TCCCCAGA‑CCCTAATAAAGTTCGGTGTCAACCCCCCCTTTAcccctttgtat                                                                                                                                                                                                    >  1:678883/1‑231 (MQ=255)                                          tGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCCTTTGTTTTCCGCCCCC‑CCCCTGAG‑TAATTAATACGCCAGACTCCGCTACCA‑GGA‑ACCCACCCGGTTCAATTT‑TCAGTTTCACACGCGGGCCCCAA‑A‑GGCGGGGGATTTGCTCACACGCGGC‑AC‑TTCACAAG‑CCCCAAAAATGTTCGGTAGCAACCCCCCCCTCAaccccttcttaaagca                                                                                                                                                                                               >  2:47479/1‑226 (MQ=255)                                          tGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGtt                                                                                                                                                                                                                                   <  1:72148/203‑1 (MQ=255)                                               ggCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCCCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGCCCCCGGTCAATTT‑CCAAGTTCCCAGGGTGGCCCCAA‑T‑GGCGGGGGAGTGGCTCACCCCCGGCAAT‑CCCCCCGC‑CCAAAAAAAT‑TTCCGTATTAACCGCCCCTTCccccctttttaaaagaaaaaaa                                                                                                                                                                                          >  2:1086993/1‑221 (MQ=255)                                                            ttggcacccgcacccacccaccaAAAAGGGGC‑G‑CCCCCCATTTGTGTTACGCGCCG‑CACCCGAGATAATAAAGCCGCCAGACGGGTGTCCCG‑GGA‑GGGGCCCGGGTTAAATTG‑CAAGTTTCCCAGGGGGGCCAAAA‑T‑GGCGGGCGTTTTGCCCCCCCGCGCC‑TT‑CCCACAGA‑CCATAATAAGGTCCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCAccgcccg                                                                                                                                                                             <  2:238849/222‑1 (MQ=255)                                                                  cccGGCACCAACCAGCACAAAGGGGCCG‑CTCCCAATTTGGTGGACCCCCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGaaa                                                                                                                                                                       <  1:297543/239‑1 (MQ=255)                                                                  caccgCATCACCCA‑CAAAAAGGGGGGGCCTCCCCTTTTGTTTGAGGCTCGG‑ACACTGAT‑TACTTAACCCGCCAGCCGGCGC‑ACCC‑GGAGAGCGACCGGGTTCATTTC‑CGGGTTTAAAGGGGTGGCCTAAA‑T‑GGCGGGCGTTTTGCCCCCGCGCGGC‑TT‑CCCACAGA‑CAATAATAATGTTGGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCCCCGCCCGCAGaaa                                                                                                                                                                       <  2:661931/235‑1 (MQ=255)                                                                     ggcaccccccAGCCCAAAGTGGC‑G‑CCCCCCATTTTGTTGGCGGCCCGCCACATATG‑AAATTAAGCCGCCAGTC‑GGGCCACCAGGGA‑AGCGCCCCGTTTAAATTC‑CCAGTTCCCCGGGGGGCCCAAAA‑T‑GGCGGCCGTTTGCCTCCCGCGCGGC‑TT‑CCCACAGA‑CAATAATAATGTCCGTTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAg                                                                                                                                                                    <  2:644915/231‑1 (MQ=255)                                                                     ggCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCCACG‑CACCTGAG‑TAATTAATAAGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGCTCAATTT‑CCAGTTTTACAGGGGGGCCCCAA‑A‑TGCGGGCGATTTTGCCGCCCGCGGC‑AT‑TCCCCACA‑CCAAAAAAAAGTTCGGTCTCAACACCCCCTTCACCCCGTTGATAAGAAACAAAAAAACCACCCCCCCCACAAATa                                                                                                                                                                     >  2:118239/1‑238 (MQ=255)                                                                                                      ttttttggccccTCCC‑CCCCGGAG‑TAATTAATCCGCCGGGCTGCGCTCCCA‑GGA‑AGCGACCCGGTTCAATTC‑CCAGTTTCACAGGGTGGCCAAAA‑T‑GGGGGGCGATTTCCTCCCGCGGGCA‑AT‑TCCA‑AGACCCAAAATAATGTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGa                                                                                                                                    <  2:1014075/228‑1 (MQ=255)                                                                                                         gtttgACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGGGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCTTAATAATGTTCGGTATAAACACCCCGTTCACCCGATTGGTAAGGAACAGAAATCCCACCGCCCGCAGAAATAGCCCGTTCATACTGACCGAAtgacataacatacgg                                                                                                                                  >  2:892820/1‑224 (MQ=255)                                                                                                          tttgACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGTCGATTTGCTCACCCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGc                                                                                                                                 >  1:761320/1‑238 (MQ=255)                                                                                                            tgACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACCGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCCCCCGCAGAAATAGCCTGTAAATCCTGACCGAAATCAAtaccccgacccccc                                                                                                                              >  2:927671/1‑226 (MQ=255)                                                                                                            tgACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCCTAATAATGTTGGGTATCAACCGCCTCTTCACCCTTTTGTAAAGGAACCGAAAGCCCACCCCCCCCAGAAAAAGCCATTAAATCCTGACCGAATTCGagaaccgaaacccaa                                                                                                                              >  2:520897/1‑225 (MQ=255)                                                                                                             gccGCTCCG‑CCCCTGAG‑TAATTAAGCCGCCAGACTGCCCTACCA‑GGA‑AGCGCCCCGGTTCAATTT‑CCGGTTTCACGGGGTGGCCCAAA‑G‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTCCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCa                                                                                                                             <  1:503129/237‑1 (MQ=255)                                                                                                              aCGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCCCCAt                                                                                                                            >  1:1079525/1‑239 (MQ=255)                                                                                                               cGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGaaa                                                                                                                                                                       >  1:93959/1‑195 (MQ=255)                                                                                                                    cG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCCCCCGCAGAAAAAACCTTTAAATCCTGACCGAATTCGATAAcccgccccccccaaacacc                                                                                                                      >  2:293236/1‑222 (MQ=255)                                                                                                                       caaacgcgtatatttatCCGCCAGCCGGCGCTAGCA‑GGA‑GGCG‑CCCGTGTCAAACT‑CCCCTGTCCCAAGGGGGCCCAAA‑TGGGCGGGCAGTTGGCCCACGC‑CGGCAAT‑TCCCAAGC‑CCAGTATAATGGCCGGGAGCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGcacca                                                                                                                     <  2:591351/228‑1 (MQ=255)                                                                                                                             aG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACcacaca                                                                                                                 >  1:1093246/1‑236 (MQ=255)                                                                                                                                 aattaacCCCCCAGGCTGCGCTCCCA‑GGG‑AGCG‑CCCGTTTCTAATT‑CCCCGTTCCCAGGGTGGCCCAAAGG‑GGCGGGCGTTTTGCTCCCCGGCGCC‑TT‑CCCCAAGC‑CCATAAAATTGTCCGGTAACACCCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGtt                                                                                                           <  2:231284/239‑1 (MQ=255)                                                                                                                                     aatACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCGTCAATAGCCCGTAAATCCTGACCGAATTCGCCAACCCGACGCACCATAACACAACACACCCCTTCCaa                                                                                                       >  2:34279/1‑238 (MQ=255)                                                                                                                                     aatACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAAAAAGCCCACCGCCCGCCGAAACAGCCTGTAAATCCTGACCGAATTCGAGGACCTGAAGCACCAAAACACCACAAACCCGGTCCaa                                                                                                       >  2:365347/1‑238 (MQ=255)                                                                                                                                        gtttaaaagagacgacccA‑GGA‑AGCGCCCGGGTTAAATTT‑CCATTTTCCAGGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCCCCCGCGGGGC‑AT‑TCCACAGA‑CCAAAATAATGTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCATACACCTGTTCCAGa                                                                                                      <  2:660954/221‑1 (MQ=255)                                                                                                                                          gCCAGACTGCGCTACCA‑GGA‑AGGGACCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGA‑AT‑TCCACAGA‑CCATAATAATGATCGGCATCAACCTCCTCTTCACCCTGTTGATAAGGAACAAAAAGACCACCCCCCGCAGAAATAGCCGGTAAATCCTGACCGAATTCCACGACCTTAAGCACCATAGCACCACACACCTCGTCCCGAttg                                                                                                   >  2:164202/1‑236 (MQ=255)                                                                                                                                              ccgcgcgcgctccgagggggccccccGGTTCATTTC‑CGAGTTACACAGGGGGGCAAAAA‑T‑GGCGGGCGATTTTCCCACGCGC‑GC‑ATGCCCACAGC‑CCAAAAAAATGGTGGGTTCCACCCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCa                                                                                                <  2:678883/216‑1 (MQ=255)                                                                                                                                                 gcgccccccA‑GGA‑GGCGCCCCGGTTAATGTA‑CCAGTTCCCCAGTGTGCACAAAA‑G‑GGCGGGCTCTTTGCCCACGGCCGGC‑AT‑TCCACAAA‑CAAAAAAAATTTGCGGATCCAACCGCCTCTTCCCCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTTACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATc                                                                                           <  2:761320/232‑1 (MQ=255)                                                                                                                                                   ggcttcccgggA‑AGCGGCCCGTTTCAATTT‑CCAGTTTCCAGGGGTTGCCAAAA‑T‑GGCGGGCGATTGTCTCACGCGGGGC‑AT‑TCCAAAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCa                                                                                          <  1:100210/230‑1 (MQ=255)                                                                                                                                                         acaaggagccgccccGGTTAATTTT‑CCAGTTTCACGGGGTGCCCAAAA‑G‑GGCGGGCAATTGCCTCACGCGCGGC‑AT‑TCCACAGC‑CAATAATAATGTTGGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGcc                                                                                     <  1:228692/227‑1 (MQ=255)                                                                                                                                                                 gCGACCCGGTTCAATTT‑CCAGTTTCAGAGGGTGCCCAAAA‑T‑GGCGGGCGATTTGCTCACGCGGGGC‑AT‑TCCACAGC‑CCAAAATAATGTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCa                                                                              <  1:161912/238‑1 (MQ=255)                                                                                                                                                                    gccgcgTTTAATTTCCCGGTTGCACAGGGTGCCCAAAA‑T‑GCGGGGCGTTTTGCCCAGGCGGGGC‑AT‑CCACCAGA‑CCAAAAAAAGGTTGGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAAt                                                                            <  2:172804/233‑1 (MQ=255)                                                                                                                                                                    aCCGGGTTCAATTC‑CCATTTTCACAGGGTGGCCCAAA‑G‑GGCGGCCGTTTGGCTCCCGCGCGGC‑AT‑CCCACAGA‑CCAAAAAAATTTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAAtgt                                                                          <  1:1086993/239‑1 (MQ=255)                                                                                                                                                                    aCCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCATACACCTGTTCCAGATGGGCATAAt                                                                                            >  1:660954/1‑221 (MQ=255)                                                                                                                                                                         gTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACCACCCATAGCACCACATACCTGTTCCAGATGGGCATAATCAAACCCAAAACCacagcgcagg                                                                     >  1:1074225/1‑230 (MQ=255)                                                                                                                                                                                                ggggCCAAAA‑T‑GGCGGGCGTTTTGCTCCCCGCCGGC‑AT‑CCCCAAAA‑CCAAAAAAATGTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCa                                                  <  1:402422/234‑1 (MQ=255)                                                                                                                                                                                                  ggCCCAAA‑G‑GGGGGGCGATTGCCACACGGCTGGC‑AT‑CCCCCAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCCCCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAgatga                                             <  1:494481/239‑1 (MQ=255)                                                                                                                                                                                                   gCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAgatgat                                            >  1:220610/1‑239 (MQ=255)                                                                                                                                                                                                    cccAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑CT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTTTAAATCCTGACCGAATTCGCTGACCTGACGCACCATAGCACCCCACCCCTGGTCCAGATGGGCATAATCAACCCCAGAACCAATGAGCATGTGAATGACGACAAGATTCAgatgaaa                                           >  2:256723/1‑237 (MQ=255)                                                                                                                                                                                                     caaaa‑G‑GGCGGGTGATTTGCTCCCGCGCGGC‑AT‑CCCACTGA‑CCATAATAAGGTGCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGatgatg                                           <  1:255281/236‑1 (MQ=255)                                                                                                                                                                                                       aaa‑G‑GGCGGGCGATTTCCCCCCGCGGGGC‑AT‑TCCACAGA‑CCAAAATAAGGTTCGGTAACAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACg                                        <  1:293236/239‑1 (MQ=255)                                                                                                                                                                                                         aattgccgggctttttccccccccgcgGC‑TT‑TCCACGAA‑CCATAAAAATGTCCGGTATCAACCGCCTTTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACg                                        <  2:1079525/216‑1 (MQ=255)                                                                                                                                                                                                         a‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGt                                                                                                                                                               >  2:358408/1‑118 (MQ=255)                                                                                                                                                                                                         a‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGt                                                                                                                                                               <  1:358408/118‑1 (MQ=255)                                                                                                                                                                                                         a‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGCCGtt                                      >  1:434318/1‑239 (MQ=255)                                                                                                                                                                                                         a‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAAAAGCAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCGGACCGAATTCGATGCCCCGACGCACCATAGCCCCCCACACCTGTTCCACACGGGCAAAATCAACGCCAAAACCAAAGCGCATGTGAAAGCCGACAAgtcccaaaagagagcacg                                      >  2:90578/1‑222 (MQ=255)                                                                                                                                                                                                                              cgcgcgGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGa                                                                                                      <  1:99215/156‑1 (MQ=255)                                                                                                                                                                                                                              cgcgcgGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGa                                                                                                      >  2:99215/1‑156 (MQ=255)                                                                                                                                                                                                                                      cttatcccaaGA‑CAATAATAATTTTGGGTATCAACCGCCTCTTCCCCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCa              <  1:520897/229‑1 (MQ=255)                                                                                                                                                                                                                                            cacaGA‑CCATAATAATGTTCGGTATCAACCGCCTATTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt       <  1:365347/239‑1 (MQ=255)                                                                                                                                                                                                                                               acagccAAAAAAATGTTCGGTAACAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt       <  1:327942/233‑1 (MQ=255)                                                                                                                                                                                                                                                 a‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGATCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATCATCAACTCCAGAACCAATGTTCTCATCAATCCCGACCAGATGCAGATAATGACGTTGTATCACGTCCAGAGCGGCGGGCAGATCGgtgtccg  >  1:510731/1‑236 (MQ=255)                                                                                                                                                                                                                                                 a‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCACCGCCAGATGATAACGTAGAATACCATCCAGTGCGGCGGGCACATCgggggaca  >  2:104129/1‑233 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           TCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCCGGTTCAATTT‑CCAGTTTCACAGGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGGTGTACC  >  NC_000913/2977930‑2978398

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