Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	91,904	G→A	E164E (GAG→GAA)	ftsI →	peptidoglycan DD‑transpeptidase FtsI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	91,904	0	G	A	97.6%	139.2 / ‑5.2	41	E164E (GAG→GAA)	ftsI	peptidoglycan DD‑transpeptidase FtsI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/22);  minor base C (1/0);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

GGTGAGCGTGCCGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATC  >  NC_000913/91661‑92148                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     ggTGAGCGTGCCGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGcc                                                                                                                                                                                                                                               <  2:212139/251‑1 (MQ=255)   tGAGCGTGCCGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGt                                                                                                                                                                                                                                             >  1:5504/1‑251 (MQ=255)           ccGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTAtcc                                                                                                                                                                                                                                     <  2:139302/251‑1 (MQ=255)               aaaaaGCGATTTGGTCAGACCCGAAAGAAGTGCATGACGCGGGCTGTATCAGCGTCTGTGACCGCTGGAAGGCGTTGGCTAACGCGCTCAATATTCCGTTGGATCAGCTTTCAGCCCCCATTAACGTCTACCCGAAATGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGAC‑TGTGCAGACTACATCAAAAAACTGAAACTGCCGGGGATTCATTTGCGTGAAGAATCTCGCCGTTATTatccatcc                                                                                                                                                                                                                                 <  1:301031/250‑1 (MQ=255)                       aTTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGt                                                                                                                                                                                                                         >  2:332205/1‑251 (MQ=255)                                 aCCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCAc                                                                                                                                                                                                               <  1:73406/251‑1 (MQ=255)                                   ccGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCt                                                                                                                                                                                                             <  2:343300/251‑1 (MQ=255)                                         gaagTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCAttgg                                                                                                                                                                                                       >  2:63064/1‑248 (MQ=255)                                                              aTCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCa                                                                                                                                                                                                           >  2:198316/1‑226 (MQ=255)                                                              aTCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCa                                                                                                                                                                                                           <  1:198316/226‑1 (MQ=255)                                                                                        gcTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGtt                                                                                                                                                        >  2:268234/1‑251 (MQ=255)                                                                                               aaCGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAAtct                                                                                                                                                                                                                                                   <  1:324598/153‑1 (MQ=255)                                                                                               aaCGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAAtct                                                                                                                                                                                                                                                   >  2:324598/1‑153 (MQ=255)                                                                                                aCGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATaaa                                                                                                                                                <  1:285230/251‑1 (MQ=255)                                                                                                  gcgcTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATg                                                                                                                                              <  2:225381/251‑1 (MQ=255)                                                                                                      tCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCtt                                                                                                                                          >  2:213534/1‑251 (MQ=255)                                                                                                           atTCCGCTGTATTAGCTTTCAGCCCGAATTAACGCCAACCCTAAAGGTCGCTTTATTTATCGGGCCCGTAAGGTGAACCCTTACATG‑GCGGACTACATCAAAAACCTGAAACTGCCTGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCTATAAATGGCTTACCgg                                                                                                                                     <  2:371747/251‑1 (MQ=255)                                                                                                           atTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCgg                                                                                                                                     <  2:24234/251‑1 (MQ=255)                                                                                                                   ggATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGcc                                                                                                                               <  1:164656/249‑1 (MQ=255)                                                                                                                   ggATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCg                                                                                                                              <  1:63064/250‑1 (MQ=255)                                                                                                                    gATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCggg                                                                                                                            >  1:71972/1‑251 (MQ=255)                                                                                                                     aTCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGt                                                                                                                           <  2:275907/251‑1 (MQ=255)                                                                                                                            ttCAGCCCGCCTTTCCGCCACCCCGAACGGGCTCTTTCTTTATCTGGCGCGCCGGGTGAACCCTGACATG‑TCGGACTACATCAAAAACCTGAACCTGCCGTGGATTCCTCTACGTTAAGAATCTCGCCGTTATTATCCATCCGTCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGCGCGTTGAGAAGAGTTTCGATATATGGCTTACCGGGCGGCCGGGTGAgcg                                                                                                                      <  2:288300/249‑1 (MQ=255)                                                                                                                            ttCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGt                                                                                                                                                                                                                         <  1:303918/150‑1 (MQ=255)                                                                                                                            ttCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGt                                                                                                                                                                                                                         >  2:303918/1‑150 (MQ=255)                                                                                                                                              cAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCAGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCt                                                                                                  <  1:164059/251‑1 (MQ=255)                                                                                                                                                aCCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTAt                                                                                                <  2:5504/251‑1 (MQ=255)                                                                                                                                                                   ttatatggcacgccaggagacaccTGACCTG‑GCGGCCTACATCAAAAACCTGCAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATACGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCCGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGat                                                                              >  2:310942/23‑250 (MQ=255)                                                                                                                                                                     atCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGTGTAATTGAAGatat                                                                            <  2:50682/250‑1 (MQ=255)                                                                                                                                                                        tGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGGGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATttct                                                                        >  2:304429/1‑251 (MQ=255)                                                                                                                                                                         ggCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATAttt                                                                          <  1:305816/248‑1 (MQ=255)                                                                                                                                                                            gggTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATCGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTtcttct                                                                     >  1:38856/3‑250 (MQ=255)                                                                                                                                                                                    tGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTAACCTCATTGGCTTTACTAAAGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGAGAGCGCATTGTGCGTCAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGc                                                            >  1:246714/1‑251 (MQ=255)                                                                                                                                                                                       aCCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGcc                                                           >  2:345065/1‑249 (MQ=255)                                                                                                                                                                                           tGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGctg                                                     >  1:86749/1‑249 (MQ=255)                                                                                                                                                                                             tcATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGATATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGCGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTgcg                                                   >  2:196344/2‑251 (MQ=255)                                                                                                                                                                                                                                 gggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCa                <  1:245006/251‑1 (MQ=255)                                                                                                                                                                                                                                          cTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGt       >  2:253824/1‑251 (MQ=255)                                                                                                                                                                                                                                          cTGCGCGAAGACTGTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTggg       >  2:371472/1‑250 (MQ=255)                                                                                                                                                                                                                                            gCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGttt     <  1:169325/251‑1 (MQ=255)                                                                                                                                                                                                                                               tGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGGGACTCCTCAACTCATTAGCTTTACTAACGTCGATAGTCAAGAGATTGAGGGCGTTGAGACGCGTTTCGCTAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGtttctc  >  1:371891/1‑248 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     GGTGAGCGTGCCGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATG‑GCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATC  >  NC_000913/91661‑92148

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