Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,915,009	G→A	Q176* (CAG→TAG)	proQ ←	RNA chaperone ProQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,915,009	0	G	A	97.7%	120.1 / ‑4.5	43	Q176* (CAG→TAG)	proQ	RNA chaperone ProQ
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (22/20);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

CTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGT  >  NC_000913/1914765‑1915253                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     cTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTAttcttc                                                                                                                                                                                                                                                >  1:155099/1‑251 (MQ=255)       ttGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAgg                                                                                                                                                                                                                                          >  2:279296/1‑251 (MQ=255)          cAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGGGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGc                                                                                                                                                                                                                                       <  1:242617/251‑1 (MQ=255)                    aGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTCCTGtt                                                                                                                                                                                                                             >  2:351025/1‑251 (MQ=255)                    aGCTAACGCGGTAAGCCTAAAAAACATGTTCAGGCCTAGCCCGGACTCCGCTTCCGACCACCAGGCGTTCTGCGCGCACACTCAACGACCTACCCGCTTTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGATGAAATGTCAGAAACCGGGGAGTGCTATTCTTCGCGGGGTGCTTTTAATGtt                                                                                                                                                                                                                             >  1:165424/1‑251 (MQ=255)                         aCGCGGTAAGCCTAAAAAACATGTTCATGCCTCGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTgg                                                                                                                                                                                                                        >  1:21985/1‑251 (MQ=255)                                          aaCATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCt                                                                                                                                                                                                       >  1:151702/1‑251 (MQ=255)                                              tGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTg                                                                                                                                                                                                   >  1:7428/1‑251 (MQ=255)                                                           cccGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCCCAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTATTTGGTGATTTCTAATCCGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGttt                                                                                                                                                                                       >  1:52505/1‑250 (MQ=255)                                                                    cgtgtCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCGTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAgcg                                                                                                                                                                             <  2:281754/247‑1 (MQ=255)                                                                      tttCAGAACACCAGGTGTTCTGCGCGCACAATAAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTATTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGCCAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTATTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGAGCACGAGGTTTACGTTCAGCGcc                                                                                                                                                                           >  2:81046/1‑251 (MQ=255)                                                                       ttCAGAACACCAGGTGTTCTGCGCGCACAATCCAAGACATACCCGAATTAAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCct                                                                                                                                                                          >  1:157993/1‑251 (MQ=255)                                                                                          tGCGCGCACAATCAAAGACATACCAGAATTCAGCTGGACGCGGAAGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTAAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGGGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCggcgtggcgt                                                                                                                                                         >  2:299391/1‑249 (MQ=255)                                                                                            cgcgcACAATCAAAGACATACCCGAATTCAGCTGGACGCGGTCGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGTCGTAGTc                                                                                                                                                     <  1:260481/251‑1 (MQ=255)                                                                                                         aGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGc                                                                                                                                        >  1:298553/1‑251 (MQ=255)                                                                                                               aCCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgc                                                                                                                                  >  2:254172/1‑251 (MQ=255)                                                                                                                   gAATTCAACTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACACGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgcggcg                                                                                                                              >  1:266109/1‑251 (MQ=255)                                                                                                                    aaTTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgcggcgc                                                                                                                             >  2:32065/1‑251 (MQ=255)                                                                                                                     aTTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgctgcgcg                                                                                                                            <  1:55402/251‑1 (MQ=255)                                                                                                                        cAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTg                                                                                                                         >  2:92088/1‑251 (MQ=255)                                                                                                                        cAGCGGGACGCGGACGCCGGCTTTGGTTATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGt                                                                                                                          <  1:68246/250‑1 (MQ=255)                                                                                                                              gacgcggccgcCGTCTTTGGTGATTTCTAATCCGGTGGCATCCCTCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGACACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGtctt                                                                                                                   >  2:326032/1‑251 (MQ=255)                                                                                                                                 gcggacgcCGTCTTTGGTTATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCTCGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCTCGGCTCTGTGCGTcttctt                                                                                                                <  2:238822/251‑1 (MQ=255)                                                                                                                                 gcggacgcAGTCTTTGCTGATTTCTAATACGGCGGCATCCCTCGCGTTTTGCCCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGCAATGTCAGCGACAGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCCCGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTTGCGTAGTCGGACGTGGCTTGCGGTCGCGGCCCGGTGCGTcttctt                                                                                                                >  2:183824/1‑251 (MQ=255)                                                                                                                                       gcCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCAc                                                                                                          <  1:92088/251‑1 (MQ=255)                                                                                                                                            cgttGGTGATTTCTAATGCGGTGGCATCCAGCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGGCAGAGCTTAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGGTTTTGGCGCTTTCTCTACCGGCTTTTGCTCACGATGCTTACGTCCAGCGCCTGCTATGCGGCGTGGCTTATTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACcagc                                                                                                      <  2:353989/248‑1 (MQ=255)                                                                                                                                              ttGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTGCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTGTTGCGCTCGAGGTTTACGTTCATCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCGCACcagc                                                                                                      <  2:210241/248‑1 (MQ=255)                                                                                                                                                               ggTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGACGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTTCGGCAGCCTCGCGtt                                                                                  >  2:204083/1‑251 (MQ=255)                                                                                                                                                                  ggCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGttttt                                                                               >  1:79115/1‑251 (MQ=255)                                                                                                                                                                         atcGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTg                                                                        <  2:354977/251‑1 (MQ=255)                                                                                                                                                                                     cccGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACg                                                            <  1:32065/251‑1 (MQ=255)                                                                                                                                                                                       cGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTGTTTCGCTTGCTGTTCATCACGCt                                                          <  2:151702/251‑1 (MQ=255)                                                                                                                                                                                             cACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGcc                                                     <  1:120674/250‑1 (MQ=255)                                                                                                                                                                                                                aCAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGTTGCTTTTTCTGTTTTTGGCGCTTTCTCGACCGGCGTTTGCGCACGAGGTTGACGGTCAGCGCCTTCTTTGCGGCGTGGCGTATTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTGCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGcttc                                 <  2:157993/251‑1 (MQ=255)                                                                                                                                                                                                                       gagCTGAAAGGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAgc                          <  2:6792/251‑1 (MQ=255)                                                                                                                                                                                                                         gCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAgctg                        <  1:279296/251‑1 (MQ=255)                                                                                                                                                                                                                         gCTGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAgctg                        <  1:86635/251‑1 (MQ=255)                                                                                                                                                                                                                           tGAAATGTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAgctgct                      >  1:280155/1‑251 (MQ=255)                                                                                                                                                                                                                                 gTCAGAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAg                <  2:7428/251‑1 (MQ=255)                                                                                                                                                                                                                                     gaatccGGGGAGTGCTATTCTTCGCGGGGTGCTTTTACTGTTTTTGGCGCTTTCTCTCCCGGCTTTTGCGCACGAGGGTTACGTTCAGCGCCTGCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGGCTGCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTGCAAGCTGCTTGCGAGCATg            <  2:10291/247‑1 (MQ=255)                                                                                                                                                                                                                                     gAAACCGGGGTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATg            >  2:289224/1‑251 (MQ=255)                                                                                                                                                                                                                                             ggTGTGCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACAt    <  1:260602/251‑1 (MQ=255)                                                                                                                                                                                                                                                gtgCTATTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGt  <  1:351025/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     CTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGT  >  NC_000913/1914765‑1915253

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