Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,153,836	G→A	Q338* (CAA→TAA)	argE ←	acetylornithine deacetylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,153,836	0	G	A	82.9%	98.5 / 13.4	41	Q338* (CAA→TAA)	argE	acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base G (2/5);  new base A (18/16);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCCAGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGCAAACCATTAAGTTCATTGAGTGTCATGCCAGGCAGCGGACGAATATCCATATGCAACTC  >  NC_000913/4153600‑4154076                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 ggATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCtt                                                                                                                                                                                                                                    >  2:225475/1‑251 (MQ=255)   aTAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACCCCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCg                                                                                                                                                                                                                                  >  2:345073/1‑251 (MQ=255)   aTAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCg                                                                                                                                                                                                                                  >  2:143558/1‑251 (MQ=255)           ggcAAGCATCGCAGCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCGGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACagta                                                                                                                                                                                                                          <  2:71286/250‑1 (MQ=255)                  atcgcatcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTcacc                                                                                                                                                                                                                   >  1:26622/1‑251 (MQ=255)                       atcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCAGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTc                                                                                                                                                                                                              >  1:204217/1‑251 (MQ=255)                              attattGCGTCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGttttt                                                                                                                                                                                                       <  1:265434/251‑1 (MQ=255)                                tattGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGACTGATTAATTGAGCCAGGCCCCAACACCAGCGTAGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGc                                                                                                                                                                                                     >  1:303690/1‑251 (MQ=255)                                                         aGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCa                                                                                                                                                                            <  1:102082/251‑1 (MQ=255)                                                              cGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAg                                                                                                                                                                       <  1:87685/251‑1 (MQ=255)                                                                   gcggcgCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTtga                                                                                                                                                                   <  1:1820/250‑1 (MQ=255)                                                                      gcgcGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTgatgat                                                                                                                                                               <  1:241918/251‑1 (MQ=255)                                                                          gAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACATGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCgg                                                                                                                                                           >  1:49003/1‑251 (MQ=255)                                                                                cTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGc                                                                                                                                                      >  2:84375/1‑250 (MQ=255)                                                                                            aCGTTTTAATGCCAGCAAAAATGGTTAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCATCCACTTCAACCAGTTGATGATTCGGTGTGCATTCATAGCCa                                                                                                                                           <  2:454634/249‑1 (MQ=255)                                                                                            aCGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAg                                                                                                                                          <  1:208040/250‑1 (MQ=255)                                                                                            aCGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAg                                                                                                                                          <  2:49003/250‑1 (MQ=255)                                                                                               ttttAATGCCAGCAAAAATGGAGAATTACCTGGGTTATCAGTTCGTGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGAt                                                                                                                                      >  2:91247/1‑251 (MQ=255)                                                                                                      gCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGa                                                                                                                               <  2:369669/251‑1 (MQ=255)                                                                                                              aaaTGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCtc                                                                                                                       <  2:357554/251‑1 (MQ=255)                                                                                                              aaaTGGTGAATTACATGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCtc                                                                                                                       <  1:279062/251‑1 (MQ=255)                                                                                                                 tGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGGTGCAGCtcgct                                                                                                                    >  2:119404/1‑249 (MQ=255)                                                                                                                             tGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGCTGCAGCTCGTCGACCGGCAgacg                                                                                                        >  2:223937/1‑251 (MQ=255)                                                                                                                              gggTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTACGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGGCGATCGTCGACCGTCAgacg                                                                                                        >  2:141398/1‑250 (MQ=255)                                                                                                                                           gcgGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCCAGCGtt                                                                                          <  2:410627/251‑1 (MQ=255)                                                                                                                                                 gggCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCCGCTCGTCGACCGGCAGACGACCCGGCCAGCGTTCGCTc                                                                                     >  1:362398/1‑250 (MQ=255)                                                                                                                                                                     ccAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCCAGCGTTCGCTCACCGGAGCCAATGCATCGTTg                                                                <  2:95534/251‑1 (MQ=255)                                                                                                                                                                             tcatcaGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACAGCCCGGCCAGCGGTCGCTCCCCGGAGCCCATGCCTCGTTTATCAAAcc                                                        >  1:180291/1‑251 (MQ=255)                                                                                                                                                                                     ttGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCCAGCGTTTGCTCCCCGGAGCCAATGCATCGTTGaagaacggaagagacc                                                >  1:265248/1‑236 (MQ=255)                                                                                                                                                                                     ttGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCCAGCGTTCGCTCACCGGAGCCAATGCATCGTTGa                                                               <  2:265248/236‑1 (MQ=255)                                                                                                                                                                                                    aaTTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCCCGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGCAAACCATTAAGTTCATTGAGTGTCATGc                                 >  1:13767/1‑251 (MQ=255)                                                                                                                                                                                                    aaTTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGGCGACCGTCAGACGACCCGGCCAGCGGTCCCCCACCGGAGCCCATGCATCGTTTGGCAAACCATTAAGTTCAtt                                            >  1:356243/1‑240 (MQ=255)                                                                                                                                                                                                    aaTTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGAGGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCCAGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGCAACCCATTAAGTTCAtt                                            <  2:356243/240‑1 (MQ=255)                                                                                                                                                                                                       tGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCCAGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGCAAACCATTAAGTTCATTGAGTGTCATGCCa                               <  1:200142/250‑1 (MQ=255)                                                                                                                                                                                                          gCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATGCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCCAGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGCAAACCATTAAGTTCATTGAGTGTCATGCCAGGCa                           <  1:481622/251‑1 (MQ=255)                                                                                                                                                                                                           ccAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCCAGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGCAAACCCATAAGTTCATTGAGGGTCATGCCAGGc                            >  1:140092/1‑249 (MQ=255)                                                                                                                                                                                                                      cacCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCCAGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGCAAACCATTAAGTTCATTGAGTGTCATGCCAGGCAGCGGACGAatat               <  1:491339/251‑1 (MQ=255)                                                                                                                                                                                                                             gTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGCCCGTCAGACGACCCGGCCAGCGTTCGCTCCACGGAGCCAATGCATCGGTGAGCAAACCATTAAGGTCATTGAgtgg                                      >  1:321543/1‑220 (MQ=255)                                                                                                                                                                                                                             gTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCCAGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGCAAACCATTAAGTTCATTGAgtgt                                      <  2:321543/221‑1 (MQ=255)                                                                                                                                                                                                                                   cATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCCAGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGCAAACCATTAAGTTCATTGATTGTCATGCCAGGCAGCGGACGAAGATCCATATGCAACTc  <  1:141398/251‑1 (MQ=255)                                                                                                                                                                                                                                   cATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGAAGACCCGGCCAGCGTTCGCTCCCCGGAGCCAATGCATCGTTGAGCAAACCATTAAGTTCATTGAGTGTCATGCCAGGCAGCGGACGAATATCCATATGCAACTc  >  1:231747/1‑251 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 GGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCCAGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGCAAACCATTAAGTTCATTGAGTGTCATGCCAGGCAGCGGACGAATATCCATATGCAACTC  >  NC_000913/4153600‑4154076

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