Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	91,949	C→T	I179I (ATC→ATT)	ftsI →	peptidoglycan DD‑transpeptidase FtsI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	91,949	0	C	T	97.7%	114.4 / ‑4.2	44	I179I (ATC→ATT)	ftsI	peptidoglycan DD‑transpeptidase FtsI
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (19/24);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

TGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGG‑CTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTG  >  NC_000913/91705‑92190                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  ggcactccGCTGGCGGTATCAGCGTCGGTGACAGCTGGAAGGCGCTGG‑CTCACTCTCTCAATATCCCGCTGGATCAGCTTTCCGCCCGCATTAACGCCAACCCGAAAGGTCGTTGTATTTATCTGGCGCGTCAGTTGACCCCTGACATTGCTGACTCCATCAAAAAAATGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCAGCCGGCGAAGTGACTGCTCACCTTATTGGCttt                                                                                                                                                                                                                                             <  2:190604/244‑1 (MQ=255)           gtgtgtttgcagagccGGTGACCGCTGGAAGGCG‑TGGTCTCACGCGCTCAATATTCCGCTGGTTCAGCTTTCATACCGTATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGGCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCCTCCGGCGAAGTGACTGCGCTCCTCATTGGCTTTACTAACGTc                                                                                                                                                                                                                                    <  2:221609/235‑1 (MQ=255)                                  gCTGGAAGGCGCTGG‑CTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCAGTAACGCCAACCCGAAAGGGCGCTTTATGTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACGGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGtt                                                                                                                                                                                                            <  2:461516/251‑1 (MQ=255)                                    gggAAGGCGCTTG‑CTAACGCGCTCAATATTCCGCTTGATCAGCTTTCAGCCCTCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTCACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCTTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGTTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTga                                                                                                                                                                                                          <  1:420124/250‑1 (MQ=255)                                       aaGGCGCTGG‑CTAACGAGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGACGTTCTTATCCATCCGGCGAAGTGACAGCTCTCCTCATTGGCTTTACTAACGTCGTTAGTCCAGGGATTGAGGGCGTTgaga                                                                                                                                                                                                        >  1:163090/1‑250 (MQ=255)                                                    aaCGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGAt                                                                                                                                                                                             <  1:337149/249‑1 (MQ=255)                                                          cTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCt                                                                                                                                                                                     >  1:362644/1‑251 (MQ=255)                                                           tCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCAGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATc                                                                                                                                                                                                                                                   >  2:66756/1‑188 (MQ=255)                                                           tCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCAGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATc                                                                                                                                                                                                                                                   <  1:66756/188‑1 (MQ=255)                                                            cAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGCCTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATACATGGCtt                                                                                                                                                                                    >  2:42944/1‑250 (MQ=255)                                                             actatTCCGCTGGAACAGCTTTCAGCCCGCATTAAGGCCAAACCGTCAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCAGTTATTATCCATCCGGCGAAATGACTGCTCACCTCATTGGCTTACCTACCGCCGATAGTCAAGGGGTTGAGGGCGTTCagaagtatatcgataacaggctaa                                                                                                                                                                                   >  2:161372/3‑231 (MQ=255)                                                                   ccGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCa                                                                                                                                                                            <  1:28649/251‑1 (MQ=255)                                                                                   cAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCAtt                                                                                                                                                            >  1:199488/1‑251 (MQ=255)                                                                                   cAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTAATTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTAAATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTAGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGACTTACCGGGCAGCCGGGTGCGGGCAtt                                                                                                                                                            >  2:314244/1‑251 (MQ=255)                                                                                    aGCCCGCATTAACACCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGTGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAATAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATtg                                                                                                                                                           <  1:204362/251‑1 (MQ=255)                                                                                      cccGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGGCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATtgtg                                                                                                                                                         <  1:215398/251‑1 (MQ=255)                                                                                          cATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTa                                                                                                                                                     <  1:176689/251‑1 (MQ=255)                                                                                              aaCGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTaa                                                                                                                                                    >  1:59638/1‑248 (MQ=255)                                                                                                         gAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTc                                                                                                                                      >  1:397548/1‑251 (MQ=255)                                                                                                                    ttatttatCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGa                                                                                                                            <  1:42944/250‑1 (MQ=255)                                                                                                                                     cAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGCCTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTg                                                                                                           >  2:272427/1‑250 (MQ=255)                                                                                                                                         tGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGc                                                                                                      >  1:426016/1‑251 (MQ=255)                                                                                                                                         acgctcttccgatctcGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGAc                                                                                                                                                >  1:361199/16‑209 (MQ=255)                                                                                                                                         acgctcttccgatctcGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGAc                                                                                                                                                <  2:361199/194‑1 (MQ=255)                                                                                                                                            aCCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAg                                                                                                   <  1:195099/251‑1 (MQ=255)                                                                                                                                                 gACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGctg                                                                                               <  1:51037/250‑3 (MQ=255)                                                                                                                                                       gCGGAATACATCAAACAACTGAAACCGCCGGGGTTTCATCTGCGTGAAGAATCTCTCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTTCTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACaa                                                                                        <  2:121433/251‑1 (MQ=255)                                                                                                                                                           aCTACATCAAAAAACGGAAACTGCCGGGGATTCATCTGCGGGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAAGTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTg                                                                                    <  2:397548/251‑1 (MQ=255)                                                                                                                                                                            aaactgCCGGGGATTCATCTGCGTGAAGAATCGCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCCCAACCTGGCGCTGAGTATtgatga                                                                   <  2:355140/251‑1 (MQ=255)                                                                                                                                                                              actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGCGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATtgatgacc                                                                 >  2:22779/1‑249 (MQ=255)                                                                                                                                                                                 gCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATACATGGCTTACCGGGCAACCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTAGCGTAATTGAAGATATTTCTTCAACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGcc                                                              >  2:251837/1‑251 (MQ=255)                                                                                                                                                                                 gCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACg                                                                <  2:603/249‑1 (MQ=255)                                                                                                                                                                                          tCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGgcgc                                                     >  1:73602/1‑251 (MQ=255)                                                                                                                                                                                                gCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGATTATTGATGAACGCCTGCAGGCGCTGGttt                                               <  1:381720/251‑1 (MQ=255)                                                                                                                                                                                                  gTGAAGAATCTGGCCGTTATTATCCATCCGGCGAAGTGACTGCTAGCCGCAATGGCTTTACTAACGGCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGGCAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTAt                                             <  2:425796/251‑1 (MQ=255)                                                                                                                                                                                                    gaagaaTCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATcg                                           >  1:1269/1‑251 (MQ=255)                                                                                                                                                                                                             cGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGaa                                  <  1:241318/251‑1 (MQ=255)                                                                                                                                                                                                               ccGTTATTATCCATCCGGCGACGTGACTGCTCTCCTCAGTGGCTTTTCTAACGTCGATAGTCAAGTGAGTGAGGGCGTTGAGAAGAGGTTCGATAAATGGCTTACCGTGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTTACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCTCTGGTTGTTCGCTAACTGaaca                                <  2:290343/251‑1 (MQ=255)                                                                                                                                                                                                                       atccatccGGCGAAGTGACTGCTCCCCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGAAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGt                         <  2:332031/250‑1 (MQ=255)                                                                                                                                                                                                                           acccgGCGAAGGGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCGTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGGCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCTCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCt                    <  2:124051/249‑1 (MQ=255)                                                                                                                                                                                                                             ccgGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGAGGCCttt                  >  1:407900/1‑251 (MQ=255)                                                                                                                                                                                                                                  gAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTTATGAACGCCTGCAGGCGCTGGTTTATCGCGACCTGAACAATGAGGTGCCCTTTAACa              >  2:217775/1‑250 (MQ=255)                                                                                                                                                                                                                                  gAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACa              >  1:119770/1‑250 (MQ=255)                                                                                                                                                                                                                                   aaGTGACTGCTCACCGCATTGGCTTTACTAACGACGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTGCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATGGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAg            <  1:80955/251‑1 (MQ=255)                                                                                                                                                                                                                                       gACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGCGGGCGTTGAGAAGTGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTg        >  2:267406/1‑251 (MQ=255)                                                                                                                                                                                                                                             ccACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACTGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTg  <  2:82771/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  TGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGG‑CTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTG  >  NC_000913/91705‑92190

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