Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,339,162	C→T	D87N (GAC→AAC)	gyrA ←	DNA gyrase subunit A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,339,162	0	C	T	98.2%	187.6 / NA	55	D87N (GAC→AAC)	gyrA	DNA gyrase subunit A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (24/30);  minor base G (0/1);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

ACCAGCAGGTTAGGAATTTTGGTTGGCATGACGTCCGGAATTTTTTCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGACCGACATCGCATAATCCAGATAGGAGCTCTTCAGCTCTTCCTCAATGTTGACCGGTGTA  >  NC_000913/2338918‑2339397                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             aCCAGCAGGTTAGGAATTTTGGTTGGCATGACGTCCGGAATTTTTTCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGAc                                                                                                                                                                                                                                        <  1:128153/251‑1 (MQ=255)  ccAGCAGGTTAGGAATTTTGGTTGGCATGACGTCCGGAATTTTTTCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCAGGCGGATGATCGTGTTATAGAcc                                                                                                                                                                                                                                       <  1:403236/251‑1 (MQ=255)      cagGTTAGGAATTTTGGTTGGCATGACGTCCGGAATTTTTTCCGGGCCGTCATATTTATCAACGAAATCGACCGGCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGAccgccg                                                                                                                                                                                                                                   <  2:469551/251‑1 (MQ=255)        ggTTAGGAATTTTGGTTGGCATGACGTCCGGAATTTTTTCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCAAATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGAccgccg                                                                                                                                                                                                                                   >  2:349717/1‑249 (MQ=255)          ttAGGAATTTTGGTTGGCATGACGTCCGGAATTTTTTCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTc                                                                                                                                                                                                                               >  2:8495/1‑251 (MQ=255)             ggAATTTTGGTTGGCATGACGTCCGGAATTTTTTCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATAAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGAGAGTGTTATAGACGGCCGAGTAAcc                                                                                                                                                                                                                            >  1:191007/1‑251 (MQ=255)                     gtttggCATGACGGCCGTAATTTGTTCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTGGTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAAGGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGAt                                                                                                                                                                                                                    <  2:96967/249‑1 (MQ=255)                     ggttggCATGACGTCCGGAATTTTT‑CCGTGCCGTCAGTGTTATCAACGAAATCGCCCGTCTCTTTTGCGAGAGCGGCCATCAGTGCATGGTCAATTTTCGCCAGACGGATTTCCGTGTAACGCATTGCCGCCGCAGTTTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGGTATAGACCGCCTAGTCACCATGCGGATg                                                                                                                                                                                                                   <  2:561151/251‑1 (MQ=255)                                           ttttCCGTGCCGTCATAGTTATCAACGAACTCGACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCAcc                                                                                                                                                                                              >  1:456605/1‑251 (MQ=255)                                               ccGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAAc                                                                                                                                                                                           <  1:203681/250‑1 (MQ=255)                                                  ggCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTGACCGATTACGTCACCAACGaca                                                                                                                                                                                       <  2:115417/250‑1 (MQ=255)                                                      gTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑ggg                                                                                                                                                                                  >  2:232614/1‑251 (MQ=255)                                                        cATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCGGTTCATGGGCAATTTTCGCCAGACGGATTTACGTATAACGCATTGCAGCCGCAGAGTCGCAGTCTGA‑AGAACCGAAGTTACACTGACAGTCTAACAGCATAAAACGCAGCGAGAATGGCTGCGCCATGCGGACGCACGAGTTATAGACCGCCGAGGCACCATGGGGATGGGATTTACCGATTCCGTCACCAACGACAA‑GGGc                                                                                                                                                                                 >  1:87125/1‑250 (MQ=255)                                                           atttATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCAGGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCCACGACAC‑GGGCAGAt                                                                                                                                                                             <  2:95829/249‑1 (MQ=255)                                                           aGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACTTCACCAACGACAC‑GGGCAGAt                                                                                                                                                                             <  2:126727/251‑1 (MQ=255)                                                               aTCACCGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTGCGCCAGACGGATTTCCGTATAACGCATGGCCGCCGCAGAGTCGCCGAC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGAttttt                                                                                                                                                                         <  2:155847/251‑1 (MQ=255)                                                                 cAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTta                                                                                                                                                                       <  1:268064/251‑1 (MQ=255)                                                                  aaCGAAATCGACCGGCTCTTTTTCGAGATCGGCCATCAGTTCAGGGGCAATTTGCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGGTACCCGGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCTGACGATCGTGGGATAGACCGCCGAGTCACCGTGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTtat                                                                                                                                                                      <  1:583670/251‑1 (MQ=255)                                                                     gAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTCCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGGC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGGCGCCATGGGGATGGTATTTACCGATTCCGTCACCAACGACAC‑GGGCAGATTTTTTATAgg                                                                                                                                                                   <  2:332039/251‑1 (MQ=255)                                                                        aTCGACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCATACGGATTTCCTTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCtt                                                                                                                                                                <  1:238489/251‑1 (MQ=255)                                                                         tCGACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCttt                                                                                                                                                               >  1:157072/1‑251 (MQ=255)                                                                           gACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGt                                                                                                                                                             >  1:263893/1‑251 (MQ=255)                                                                           gACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGt                                                                                                                                                             <  2:165882/251‑1 (MQ=255)                                                                               gTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGCTTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTccc                                                                                                                                                         >  2:291379/1‑250 (MQ=255)                                                                                tctcTTTTTCGCGATCGGCCATCAGTTCATGGGCACTTTTCGCCATACGGATTTCCGTATAACGCATTACCGCCGCCGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATCCGGCCGATCGTGTTATAGACCGCAGAGTCACCATGGGGATGGTATTTACCGATTCCGTCACCAACG‑CACCGGGCAGATTTTTTATAGGCTTTGTTcc                                                                                                                                                          >  1:283483/1‑249 (MQ=255)                                                                                      tttCGCGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTg                                                                                                                                                  >  1:180409/1‑251 (MQ=255)                                                                                      tttCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACTGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGACTGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTg                                                                                                                                                  >  1:587114/1‑251 (MQ=255)                                                                                        tCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGcc                                                                                                                                                >  1:580072/1‑251 (MQ=255)                                                                                           ggagcGGCCATCAGTTCATGGGCATTTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGGTCCAGTCATTGCCTAg                                                                                                                                             <  2:481332/247‑1 (MQ=255)                                                                                            gaTCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGt                                                                                                                                            >  1:507843/1‑251 (MQ=255)                                                                                                           cATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTa                                                                                                                              <  1:8495/250‑1 (MQ=255)                                                                                                                gcatagttCGCCAGACGGATTTCCGTAAAACTCATTGCCGCCGAATAGTGGCCGTC‑GATAGAACCGAAGTTGCCCTGACCGTCTCCCAGCATATAACGCAGCGAGAATGGCTGCGCCAGGCTGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCTCCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGCCATTGCCTAGTACGTCCATGGCGTAAAGTAc                                                                                                                        <  2:510719/245‑1 (MQ=255)                                                                                                                gCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCGAGTACGTTCATGGCGTAAAGTAc                                                                                                                        <  2:581872/251‑1 (MQ=255)                                                                                                                gCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCAGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTAc                                                                                                                        >  1:205716/1‑251 (MQ=255)                                                                                                                 cAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTAcg                                                                                                                       >  1:99958/1‑251 (MQ=255)                                                                                                                           cAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTa                                                                                                             <  1:546280/251‑1 (MQ=255)                                                                                                                            cgACGGATTTCCGTATAACGCATTGCCGCCGCAGCGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGCTTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACTCGACGGTGTAc                                                                                                            >  2:62318/2‑251 (MQ=255)                                                                                                                                ggATTTCCGTATAACGCCTTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACTACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGc                                                                                                        >  1:430221/1‑251 (MQ=255)                                                                                                                                  aTTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCGGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGGTCATGGCGTAAAGTACGCGACGGGGTACCGGCtt                                                                                                      <  2:262265/251‑1 (MQ=255)                                                                                                                                        gTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGcc                                                                                                >  2:30783/1‑251 (MQ=255)                                                                                                                                                      cgccgcAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTg                                                                                  <  1:312676/251‑1 (MQ=255)                                                                                                                                                       gccgcAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTgg                                                                                 <  2:202300/251‑1 (MQ=255)                                                                                                                                                        ccgcAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGc                                                                                <  2:41027/251‑1 (MQ=255)                                                                                                                                                                   ccGTC‑GATAGAACCGAAGTTCCCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGTTCGTGTTATAGCCCGCCGAGTCCCCTTGGGGATGGTATTTACCGATTACGTCACCAACGACTC‑GGGCATATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTTCGCGACGGTGTACCGGCTTCAGGCCATCTCGTACATCTGGCATCGCACGGcc                                                                     <  2:191007/251‑1 (MQ=255)                                                                                                                                                                    cGTC‑GATAGACACGAAGTTACCCCAAACGTCTACCAGCATATAACGCACCGCGAATGGCAGAGCCATACGGACGAACGTGATATAGACCGACGAGTCACCATGGGGATGGTATTTAACGATTACGTAACCACCGCCAC‑GGGCAGATTTTTTATAGGCTTTGTTCCTGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGAAGGTGTACCGGATTCAGGCCATCTCGGACATCTGGCAGCGAACgtcc                                                                     >  2:375260/1‑247 (MQ=255)                                                                                                                                                                            aGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATg                                                              >  2:290729/1‑250 (MQ=255)                                                                                                                                                                              aaCCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGAcc                                                           <  1:232614/251‑1 (MQ=255)                                                                                                                                                                                     gTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACTACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGACCGACATCg                                                    <  2:430221/251‑1 (MQ=255)                                                                                                                                                                                                        aGCATCTAACACAGCGAGAATGGCTGCGCCTTGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACTGCTTACGTCCACAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACACTCACCGACATCGCCTAATCCAGAAAGGagct                                 >  2:267472/1‑251 (MQ=255)                                                                                                                                                                                                                        agaATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAg                                                                                                                                                        <  2:317071/116‑1 (MQ=255)                                                                                                                                                                                                                        agaATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAg                                                                                                                                                        >  1:317071/1‑116 (MQ=255)                                                                                                                                                                                                                         gaATGGCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGACCGACATCGCATAATCCAGATAGGAGCTCTTCAGCTCTTCCTCa                 >  1:426991/1‑250 (MQ=255)                                                                                                                                                                                                                             ggCTGCGCCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGACCGACATCGCATAATCCAGATAGGAGCTCTTCAGCTCTTCCTCAATGtt            <  1:290729/251‑1 (MQ=255)                                                                                                                                                                                                                                   gcCATGCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCAATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGACCGACATCGCATAATCCAGATAGGAGCTCTTCAGCTCTTCCTCAATGTTGACCgg      <  1:464899/251‑1 (MQ=255)                                                                                                                                                                                                                                        gCGGACGATCGTGTTATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGACCGACATCGCATAATCCAGATAGGAGCTCTTCAGCTCTTCCTCAATGTTGACCGGTGTa  <  2:185864/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             ACCAGCAGGTTAGGAATTTTGGTTGGCATGACGTCCGGAATTTTTTCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGTTCATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTC‑GATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACAC‑GGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGACCGACATCGCATAATCCAGATAGGAGCTCTTCAGCTCTTCCTCAATGTTGACCGGTGTA  >  NC_000913/2338918‑2339397

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