Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,978,144	G→C	15.0%	R252G (CGT→GGT)	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,144	0	G	C	15.0%	147.9 / 2.0	70	R252G (CGT→GGT)	lysA	diaminopimelate decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base G (19/38);  new base C (9/3);  total (29/41)
Fisher's exact test for biased strand distribution p-value = 1.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTT  >  NC_000913/2977909‑2978362                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                              gtgttttCCATGGAACGACCCTCGGGTGCCAGGGGCATGATATGGTGGTAAGTATCGTGACTTGGCGGGCGCAT‑CAGATCGGTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTCTTTGACGCTCCG‑CCCCTGAG‑TAATTAATACGCCCGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTAACA‑GGGTGGCGCAAA‑T‑GGCGGGCGATTTGCTCCcgcg                                                                                                                                                                                                                             <  2:217450/233‑1 (MQ=255) gtgtgtGCCAGAGAACGACCATCAGCCGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCAcgcg                                                                                                                                                                                                                             <  1:964569/236‑1 (MQ=255) gcgTGTTTCATAGACCCCCCATCCGCTGCCAGGGAACTGAATAGGTGGTAACTACCGTACATTGCCGGGCGCAT‑CAGGTCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCAcgcg                                                                                                                                                                                                                             <  2:510803/238‑1 (MQ=255)     tttcccaacaaCGACCATCCGCGGCCAGGGCACTGTTTTGGGGGTAACTACCGTACATTGCCGGGCGCAT‑CAGATCATTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATGTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑a                                                                                                                                                                                                                       <  2:386494/230‑1 (MQ=255)     gtTCCAGCAAACGACCATCCGCTCGCAGGGCATGGATATGGTGGTAACTACCGTTCATTGCCGGGCGCAT‑CAGAGCGTTGAACCCGGCATCAACCGGCACAAAGTGGCGG‑CTCCCCATTGCGTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGTTACCA‑GGA‑AGCGAGCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCTAAA‑T‑GGCGGGCGATTGGCCCACGCGCGGc                                                                                                                                                                                                                         <  2:949755/238‑1 (MQ=255)     gtTCCAGAGCACGACCCTCGCGTGGCAGGGCCATGATATGGCGGTAACTACCGTACAGTGCCGGG‑GCATACAGATCGTTGACCCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCC‑CACCTGAG‑TAATTAATACGCCAGCCTGCGCGACCA‑GGA‑AGCGACCC‑GGTTCAATGT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑a                                                                                                                                                                                                                       <  2:1103137/239‑1 (MQ=255)          ggagaACGCCCCGCAGCGGCCAGGGCCTGGATAGGGTGGTAACTACCGTACATTGCCGGGCGCAT‑CAGATCGGTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCca                                                                                                                                                                                                                 <  1:193984/238‑1 (MQ=255)            agatctgccATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGGATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCcaca                                                                                                                                                                                                               <  1:938068/232‑1 (MQ=255)             gaACGGCCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCcaca                                                                                                                                                                                                               <  1:917076/238‑1 (MQ=255)             gaACGACCATCAGCTGCCAGGGCCCTGATATGGTGGTAACTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCcaca                                                                                                                                                                                                               <  1:870672/238‑1 (MQ=255)             gaACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGGAG‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCcaca                                                                                                                                                                                                               <  1:5037/238‑1 (MQ=255)             gaACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCcaca                                                                                                                                                                                                               <  1:718160/238‑1 (MQ=255)             gaACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGAATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTACGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCcaca                                                                                                                                                                                                               <  1:953709/238‑1 (MQ=255)               cccgccATCCGCTGCCAGGGCACTGATATGGTGGGAACTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCcaca                                                                                                                                                                                                               <  2:61345/232‑1 (MQ=255)               aCGATCATCAGCTGCCAGGGACCTGTTATGGTGGTAACTACCGTACATGGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCcaca                                                                                                                                                                                                               <  2:244327/236‑1 (MQ=255)                cGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGGGGCGG‑CTCCCCACTTTTTTTACGCTCCC‑CACCCGGG‑TAATTTATAAGCCCGACTGCGCTTCCC‑GAG‑AACGACCCGGGTTCAATTT‑‑CATCTTCACA‑AGGTGGCCCCAA‑A‑GGCGGGCGCGTTGCCCCCCCGCGGC‑aagccccaaaca                                                                                                                                                                                                            >  2:179834/1‑228 (MQ=255)                 gACCATGCGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCa                                                                                                                                                                                                         <  1:426008/239‑1 (MQ=255)                      tCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCataata                                                                                                                                                                                                    <  1:739386/239‑1 (MQ=255)                          cTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATAAGCCAGACTGCGCTACCC‑GGG‑AGCGACCC‑GGTTCAATTT‑CCCGTTTCACA‑GGGTGGCCCAAA‑A‑AGCGGGAGGATTGTTCACCCCCGGC‑GT‑TTCACAGA‑CCCtaacaaatt                                                                                                                                                                                                >  1:661931/1‑236 (MQ=255)                           tGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGtt                                                                                                                                                                                               <  1:179834/239‑1 (MQ=255)                             ccAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCCTTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTTCGCTACCA‑GGA‑ACCGACCC‑CGGTCAATTT‑CCAGTTTCACAGGGGTGCCCCAAA‑T‑‑GCGGGCGGTTTGGTCACCCCCGGC‑AT‑TCCACAGA‑ccccaaaaaatttcg                                                                                                                                                                                             >  1:238849/1‑226 (MQ=255)                              cAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCgg                                                                                                                                                                                            <  1:1011610/239‑1 (MQ=255)                                      ttatgtggtgtTACCTCCCGTTCTTTGCGGGGCGCAT‑CAGATTTTGGAACCCGGCATAACCCAGAACAAAGGGGGGG‑CCCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGGGCTACCA‑GGA‑AGCGCCCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGGGGCCCAAA‑T‑GGCGGGCGTTTTGCTCCCGCGCGGC‑AT‑CCCAAAGA‑CCATAATAAGGTGCGGTATCAAcc                                                                                                                                                                                    <  2:376965/234‑1 (MQ=255)                                      tGATATGGTGGTAACTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCcaca                                                                                                                                                                                                               <  1:1098808/213‑1 (MQ=255)                                             gtggtAACTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCC‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTTACA‑AGGGGGCCCACA‑T‑TGCGGGCGATTTGCTCACCCGCGGC‑CT‑TTCCCCAA‑CCATAATAATGTTCGGGATCCACccccctt                                                                                                                                                                              >  1:172804/1‑235 (MQ=255)                                                  aaCTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑TGCGGGCGAATTTCTCACCCGCGGC‑AT‑TCCCCCAA‑CCAACATAATATTCGGGATCCAACCCCTCtttcac                                                                                                                                                                         >  1:315889/1‑234 (MQ=255)                                                        cGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑CGCGACCC‑CGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTTCTCCCCCGCGGC‑CT‑TTCACAGA‑CCCTAATAATTTTCGGGATCACACGCCCCTTCAccccctttg                                                                                                                                                                  >  1:187825/1‑233 (MQ=255)                                                          tACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGGGGCCCCAA‑T‑GGCGGGCGAGTTTCTCACGCGCGGC‑AT‑TCCCCAGA‑CCCTAATAAAGTTCGGTGTCAACCCCCCCTTTAcccctttgtat                                                                                                                                                                >  1:678883/1‑231 (MQ=255)                                                               tGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCCTTTGTTTTCCGCCCCC‑CCCCTGAG‑TAATTAATACGCCAGACTCCGCTACCA‑GGA‑ACCCACCC‑GGTTCAATTT‑TCAGTTTCACA‑CGCGGGCCCCAA‑A‑GGCGGGGGATTTGCTCACACGCGGC‑AC‑TTCACAAG‑CCCCAAAAATGTTCGGTAGCAACCCCCCCCTCAaccccttcttaaagca                                                                                                                                                           >  2:47479/1‑226 (MQ=255)                                                               tGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGtt                                                                                                                                                                                               <  1:72148/203‑1 (MQ=255)                                                                    ggCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCCCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGCCCC‑CGGTCAATTT‑CCAAGTTCCCA‑GGGTGGCCCCAA‑T‑GGCGGGGGAGTGGCTCACCCCCGGCAAT‑CCCCCCGC‑CCAAAAAAAT‑TTCCGTATTAACCGCCCCTTCccccctttttaaaagaaaaaaa                                                                                                                                                      >  2:1086993/1‑221 (MQ=255)                                                                                  ttggcacccgcacccacccaccaAAAAGGGGC‑G‑CCCCCCATTTGTGTTACGCGCCG‑CACCCGAGATAATAAAGCCGCCAGACGGGTGTCCCG‑GGA‑GGGGCCCG‑GGTTAAATTG‑CAAGTTTCCCA‑GGGGGGCCAAAA‑T‑GGCGGGCGTTTTGCCCCCCCGCGCC‑TT‑CCCACAGA‑CCATAATAAGGTCCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCAccgcccg                                                                                                                                         <  2:238849/222‑1 (MQ=255)                                                                                        cccGGCACCAACCAGCACAAAGGGGCCG‑CTCCCAATTTGGTGGACCCCCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGaaa                                                                                                                                   <  1:297543/239‑1 (MQ=255)                                                                                        caccgCATCACCCA‑CAAAAAGGGGGGGCCTCCCCTTTTGTTTGAGGCTCGG‑ACACTGAT‑TACTTAACCCGCCAGCCGGCGC‑ACCC‑GGAGAGCGACCG‑GGTTCATTTC‑CGGGTTTAAAG‑GGGTGGCCTAAA‑T‑GGCGGGCGTTTTGCCCCCGCGCGGC‑TT‑CCCACAGA‑CAATAATAATGTTGGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCCCCGCCCGCAGaaa                                                                                                                                   <  2:661931/235‑1 (MQ=255)                                                                                           ggcaccccccAGCCCAAAGTGGC‑G‑CCCCCCATTTTGTTGGCGGCCCGCCACATATG‑AAATTAAGCCGCCAGTC‑GGGCCACCAGGGA‑AGCGCCCC‑GTTTAAATTC‑CCAGTTCCCCG‑GGGGGCCCAAAA‑T‑GGCGGCCGTTTGCCTCCCGCGCGGC‑TT‑CCCACAGA‑CAATAATAATGTCCGTTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAg                                                                                                                                <  2:644915/231‑1 (MQ=255)                                                                                           ggCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCCACG‑CACCTGAG‑TAATTAATAAGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGCTCAATTT‑CCAGTTTTACA‑GGGGGGCCCCAA‑A‑TGCGGGCGATTTTGCCGCCCGCGGC‑AT‑TCCCCACA‑CCAAAAAAAAGTTCGGTCTCAACACCCCCTTCACCCCGTTGATAAGAAACAAAAAAACCACCCCCCCCACAAATa                                                                                                                                 >  2:118239/1‑238 (MQ=255)                                                                                                                            ttttttggccccTCCC‑CCCCGGAG‑TAATTAATCCGCCGGGCTGCGCTCCCA‑GGA‑AGCGACCC‑GGTTCAATTC‑CCAGTTTCACA‑GGGTGGCCAAAA‑T‑GGGGGGCGATTTCCTCCCGCGGGCA‑AT‑TCCA‑AGACCCAAAATAATGTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGa                                                                                                <  2:1014075/228‑1 (MQ=255)                                                                                                                               gtttgACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGGGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCTTAATAATGTTCGGTATAAACACCCCGTTCACCCGATTGGTAAGGAACAGAAATCCCACCGCCCGCAGAAATAGCCCGTTCATACTGACCGAAtgacataacatacgg                                                                                              >  2:892820/1‑224 (MQ=255)                                                                                                                                tttgACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGTCGATTTGCTCACCCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGc                                                                                             >  1:761320/1‑238 (MQ=255)                                                                                                                                  tgACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACCGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCCCCCGCAGAAATAGCCTGTAAATCCTGACCGAAATCAAtaccccgacccccc                                                                                          >  2:927671/1‑226 (MQ=255)                                                                                                                                  tgACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCCTAATAATGTTGGGTATCAACCGCCTCTTCACCCTTTTGTAAAGGAACCGAAAGCCCACCCCCCCCAGAAAAAGCCATTAAATCCTGACCGAATTCGagaaccgaaacccaa                                                                                          >  2:520897/1‑225 (MQ=255)                                                                                                                                   gccGCTCCG‑CCCCTGAG‑TAATTAAGCCGCCAGACTGCCCTACCA‑GGA‑AGCGCCCC‑GGTTCAATTT‑CCGGTTTCACG‑GGGTGGCCCAAA‑G‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTCCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCa                                                                                         <  1:503129/237‑1 (MQ=255)                                                                                                                                    aCGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCCCCAt                                                                                        >  1:1079525/1‑239 (MQ=255)                                                                                                                                     cGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGaaa                                                                                                                                   >  1:93959/1‑195 (MQ=255)                                                                                                                                          cG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCCCCCGCAGAAAAAACCTTTAAATCCTGACCGAATTCGATAAcccgccccccccaaacacc                                                                                  >  2:293236/1‑222 (MQ=255)                                                                                                                                             caaacgcgtatatttatCCGCCAGCCGGCGCTAGCA‑GGA‑GGCG‑CCC‑GTGTCAAACT‑CCCCTGTCCCA‑AGGGGGCCCAAA‑TGGGCGGGCAGTTGGCCCACGC‑CGGCAAT‑TCCCAAGC‑CCAGTATAATGGCCGGGAGCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGcacca                                                                                 <  2:591351/228‑1 (MQ=255)                                                                                                                                                   aG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACcacaca                                                                             >  1:1093246/1‑236 (MQ=255)                                                                                                                                                       aattaacCCCCCAGGCTGCGCTCCCA‑GGG‑AGCG‑CCC‑GTTTCTAATT‑CCCCGTTCCCA‑GGGTGGCCCAAAGG‑GGCGGGCGTTTTGCTCCCCGGCGCC‑TT‑CCCCAAGC‑CCATAAAATTGTCCGGTAACACCCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGtt                                                                       <  2:231284/239‑1 (MQ=255)                                                                                                                                                           aatACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCca                                                                                                                                                                                                                 <  2:11426/98‑1 (MQ=255)                                                                                                                                                           aatACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCca                                                                                                                                                                                                                 >  1:11426/1‑98 (MQ=255)                                                                                                                                                           aatACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCGTCAATAGCCCGTAAATCCTGACCGAATTCGCCAACCCGACGCACCATAACACAACACACCCCTTCCaa                                                                   >  2:34279/1‑238 (MQ=255)                                                                                                                                                           aatACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAAAAAGCCCACCGCCCGCCGAAACAGCCTGTAAATCCTGACCGAATTCGAGGACCTGAAGCACCAAAACACCACAAACCCGGTCCaa                                                                   >  2:365347/1‑238 (MQ=255)                                                                                                                                                              gtttaaaagagacgacccA‑GGA‑AGCGCCCG‑GGTTAAATTT‑CCATTTTCCAG‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCCCCCGCGGGGC‑AT‑TCCACAGA‑CCAAAATAATGTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCATACACCTGTTCCAGa                                                                  <  2:660954/221‑1 (MQ=255)                                                                                                                                                                gCCAGACTGCGCTACCA‑GGA‑AGGGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGA‑AT‑TCCACAGA‑CCATAATAATGATCGGCATCAACCTCCTCTTCACCCTGTTGATAAGGAACAAAAAGACCACCCCCCGCAGAAATAGCCGGTAAATCCTGACCGAATTCCACGACCTTAAGCACCATAGCACCACACACCTCGTCCCGAttg                                                               >  2:164202/1‑236 (MQ=255)                                                                                                                                                                    ccgcgcgcgctccgagggggcccccc‑GGTTCATTTC‑CGAGTTACACA‑GGGGGGCAAAAA‑T‑GGCGGGCGATTTTCCCACGCGC‑GC‑ATGCCCACAGC‑CCAAAAAAATGGTGGGTTCCACCCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCa                                                            <  2:678883/216‑1 (MQ=255)                                                                                                                                                                       gcgccccccA‑GGA‑GGCGCCCC‑GGTTAATGTA‑CCAGTTCCCCA‑GTGTGCACAAAA‑G‑GGCGGGCTCTTTGCCCACGGCCGGC‑AT‑TCCACAAA‑CAAAAAAAATTTGCGGATCCAACCGCCTCTTCCCCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTTACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATc                                                       <  2:761320/232‑1 (MQ=255)                                                                                                                                                                         ggcttcccgggA‑AGCGGCCC‑GTTTCAATTT‑CCAGTTTCCAG‑GGGTTGCCAAAA‑T‑GGCGGGCGATTGTCTCACGCGGGGC‑AT‑TCCAAAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCa                                                      <  1:100210/230‑1 (MQ=255)                                                                                                                                                                               acaaggagccgcccc‑GGTTAATTTT‑CCAGTTTCACG‑GGGTGCCCAAAA‑G‑GGCGGGCAATTGCCTCACGCGCGGC‑AT‑TCCACAGC‑CAATAATAATGTTGGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGcc                                                 <  1:228692/227‑1 (MQ=255)                                                                                                                                                                                       gCGACCC‑GGTTCAATTT‑CCAGTTTCAGA‑GGGTGCCCAAAA‑T‑GGCGGGCGATTTGCTCACGCGGGGC‑AT‑TCCACAGC‑CCAAAATAATGTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCa                                          <  1:161912/238‑1 (MQ=255)                                                                                                                                                                                          aCCG‑GGTTCAATTC‑CCATTTTCACA‑GGGTGGCCCAAA‑G‑GGCGGCCGTTTGGCTCCCGCGCGGC‑AT‑CCCACAGA‑CCAAAAAAATTTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAAtgt                                      <  1:1086993/239‑1 (MQ=255)                                                                                                                                                                                          aCCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCATACACCTGTTCCAGATGGGCATAAt                                                        >  1:660954/1‑221 (MQ=255)                                                                                                                                                                                           gccgcgTTTAATTTCCCGGTTGCACA‑GGGTGCCCAAAA‑T‑GCGGGGCGTTTTGCCCAGGCGGGGC‑AT‑CCACCAGA‑CCAAAAAAAGGTTGGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAAt                                        <  2:172804/233‑1 (MQ=255)                                                                                                                                                                                                gTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑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)                                                                                                                                                                                                                                 a‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGt                                                                                                                           <  1:358408/118‑1 (MQ=255)                                                                                                                                                                                                                                 a‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGt                                                                                                                           >  2:358408/1‑118 (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                                                                  >  2:99215/1‑156 (MQ=255)                                                                                                                                                                                                                                                      cgcgcgGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGa                                                                  <  1:99215/156‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                              GCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCAT‑CAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGG‑CTCCCCATTTGTTTGACGCTCCG‑CACCTGAG‑TAATTAATACGCCAGACTGCGCTACCA‑GGA‑AGCGACCC‑GGTTCAATTT‑CCAGTTTCACA‑GGGTGGCCCAAA‑T‑GGCGGGCGATTTGCTCACGCGCGGC‑AT‑TCCACAGA‑CCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTT  >  NC_000913/2977909‑2978362

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: ‑