Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,792,833	C→A	R284I (AGA→ATA)	htrL ←	protein HtrL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,792,833	0	C	A	94.2%	148.0 / ‑2.0	52	R284I (AGA→ATA)	htrL	protein HtrL
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (21/28);  total (22/30)
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.43e-01

ATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCC  >  NC_000913/3792586‑3793083                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           atTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTg                                                                                                                                                                                                                                                         <  1:178218/251‑1 (MQ=255)      ccAAACTTATAGGTCGCGACGAGATAACGCGTTAACTCTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGTCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCATTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAgg                                                                                                                                                                                                                                                    <  1:537422/251‑1 (MQ=255)        aaaCTTATAGGGCGCGACGAGATCCCGCGTTAACACTTCTGCAAAATTCAGGACAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATAATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGa                                                                                                                                                                                                                                                   >  2:64077/1‑250 (MQ=255)                 ggTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATtctc                                                                                                                                                                                                                                           >  1:273883/1‑249 (MQ=255)                 ggTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGAGGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATGTTTCCTATCTCAGTTGGTATCAAAACGCTACCCAGTAATGGAGGGAAACCGATTTTATTGAAAGGAAAATTCTCAg                                                                                                                                                                                                                                         <  2:613136/251‑1 (MQ=255)                  gTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAgc                                                                                                                                                                                                                                        <  2:273883/251‑1 (MQ=255)                     gcgacacgatcaCGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTCCAGATAGTCGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGACAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAACTTCTCAGCGcc                                                                                                                                                                                                                                     >  2:275144/12‑251 (MQ=255)                      cgacgaGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCt                                                                                                                                                                                                                                    >  1:216870/1‑251 (MQ=255)                       gacgaGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGCTTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTg                                                                                                                                                                                                                                   >  2:97804/1‑251 (MQ=255)                        acgaGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGc                                                                                                                                                                                                                                  >  2:215442/1‑251 (MQ=255)                         cgaGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCa                                                                                                                                                                                                                                 <  2:546030/251‑1 (MQ=255)                               aaCGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTAGTATCAAAACGCGACCCACTAATGGAGGGAAACCTATTGTATTGAAAGGAAAATGCTCAGCGTCTGCAttttt                                                                                                                                                                                                                            <  2:600644/250‑1 (MQ=255)                                           cTTCTGAAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTCATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtgt                                                                                                                                                                                                               >  2:336047/1‑251 (MQ=255)                                              cTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTAc                                                                                                                                                                                                            >  2:95071/1‑251 (MQ=255)                                                       ttaGGAGAAATAACGGGTAACTTGAGGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTGTTGGCTGGTCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAAGGATGCCAAAGCATCAGTTTTCCTATCTCAGTTGGTATCAGAACCCGTCCCACTAATGAAGGGAAATCTATTGTATTGAAATGAAAAGTCTCAGCGCCTGCACTTTTGccccgtgaggttgactcctgtat                                                                                                                                                                                                    <  2:580648/248‑20 (MQ=255)                                                        cAGGATAAATAACGTGTAACGTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACGGCAGTGATATAAAAAGAAAAAAGGATGCCAAAGCATCATTTTTCCTATCTCAGTTTGTATCAAAACGCTACCCACTAATGGATGGAAACCTATTTTATGGAAAGTAAAATGCGCAGCGCCTGCATTATTGCCCCGAAAGAGTTACTCCGGAagca                                                                                                                                                                                                  <  2:549314/251‑1 (MQ=255)                                                                              gATGTCTCCTAACACCTACAGATAGTAGAAAAACTTTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCGCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTGTTGCCCCTAAAGTGTTTCTCCTGAAGCAGCGAAACAAATCGAACCAttt                                                                                                                                                                             <  1:579618/250‑1 (MQ=255)                                                                              gATGGCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCAttt                                                                                                                                                                             <  2:584539/250‑1 (MQ=255)                                                                               aTGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAACAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCAATAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGTAACAAATCGAACCATTTTc                                                                                                                                                                           >  2:229976/1‑251 (MQ=255)                                                                                       cAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTcc                                                                                                                                                                     <  2:214407/249‑1 (MQ=255)                                                                                         acacCTACAGCTAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTAATGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCCGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTCTTTTCTTGAACGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGACGCAGCGAAACAACTCGAACCATTTTCCTCTcccc                                                                                                                                                                   >  2:268695/1‑249 (MQ=255)                                                                                         acacCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTcccc                                                                                                                                                                   >  1:332043/1‑249 (MQ=255)                                                                                                                                 gCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGTGGCGTTTGTAATAGcac                                                                                                                         <  1:203480/251‑1 (MQ=255)                                                                                                                                     aaTGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATg                                                                                                                      >  2:295450/1‑250 (MQ=255)                                                                                                                                                 ttttACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCCTTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTa                                                                                                                                                                                                                     >  1:524681/1‑143 (MQ=255)                                                                                                                                                 ttttACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTa                                                                                                                                                                                                                     <  2:524681/143‑1 (MQ=255)                                                                                                                                                 ttttACTGCATTGAGATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTAAGTTGGTATAAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGCTAGTTTAGGTTGAAGAGATCGGGGCGTTTGTCATAGCACATGACAAAAATGCCCt                                                                                                         >  2:29541/1‑251 (MQ=255)                                                                                                                                                     ccTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACTCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTAGTc                                                                                                     <  1:9133/250‑1 (MQ=255)                                                                                                                                                        gCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCAt                                                                                                   <  1:268695/250‑1 (MQ=255)                                                                                                                                                        gCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTGGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATc                                                                                                  <  1:275144/251‑1 (MQ=255)                                                                                                                                                        gCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCGCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATc                                                                                                  <  2:597028/251‑1 (MQ=255)                                                                                                                                                              atatAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCcac                                                                                            >  1:456724/1‑251 (MQ=255)                                                                                                                                                              atatAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCcac                                                                                            <  2:359045/251‑1 (MQ=255)                                                                                                                                                                   aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCGCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGt                                                                                       <  1:249111/251‑1 (MQ=255)                                                                                                                                                                    aaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGtt                                                                                      <  2:332043/251‑1 (MQ=255)                                                                                                                                                                         aaaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGCAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTtattat                                                                                 >  1:389859/1‑251 (MQ=255)                                                                                                                                                                                   ccAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCCCTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAAGCGAACCATTGTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTAt                                                                       <  2:343194/251‑1 (MQ=255)                                                                                                                                                                                   ccAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAGCCATTTTCCTCTCCCCATATAGTTTAGTTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTAt                                                                       <  1:597954/251‑1 (MQ=255)                                                                                                                                                                                       aGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTAttttt                                                                   >  1:352635/1‑251 (MQ=255)                                                                                                                                                                                       aGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACGAATGGAGGGAAACCTATTTTAGTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTAttttt                                                                   <  1:428450/251‑1 (MQ=255)                                                                                                                                                                                          atcatTTTTCCTATCTCAGTTGGTATCAACACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTg                                                                 >  1:199667/1‑250 (MQ=255)                                                                                                                                                                                               ttttCCTATATAAGATGGTATCAAAACGCTACCCACTAATGGAGGGACACCTATTTTATTGAAAGGACAATTCTCAGCGCCAGCATTTATGACACAAAAGGGTTACTCATGACGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTCAAGAGATCGGCGCGTTTGTAATAGCACATGACAAAAATACACTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTc                                                           >  2:630647/1‑251 (MQ=255)                                                                                                                                                                                                   ccTATCTCAGTTGGTATCAAAACGCTACCCACTAAAGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTTCTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGAACGGTGCGTTTGTAATAGCACATGACAAAAATGCCCTTGGCATCGTCCACAATGTTATTATTTAATGTTATTTTTTGACTTTCCAGa                                                       <  1:231504/251‑1 (MQ=255)                                                                                                                                                                                                     tATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATCGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAAc                                                     >  1:258906/1‑251 (MQ=255)                                                                                                                                                                                                        ctcAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTa                                                   >  1:588330/1‑250 (MQ=255)                                                                                                                                                                                                                      aaaCACTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGAAGTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTGCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCtt                                    <  1:117144/251‑1 (MQ=255)                                                                                                                                                                                                                           caaCCCACTAATGGAAGGAAGCGTATGTTATTGAAAGGGAAGATCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCAGTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTTCACAATGTTATTATTTAGTGTTATTTGGTGACTTTCCAGATCTAATTTATAGAATTCTTTCCAt                               <  2:609390/249‑1 (MQ=255)                                                                                                                                                                                                                                     aTGGAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCtgtg                     <  1:297522/251‑1 (MQ=255)                                                                                                                                                                                                                                     aTGGAGGGAAACCTATTTTATTGAAAGGAAAATTATAAGCGCATGCATTTTTGCCCCTAAAGTGTTACTCATGAAGCAGCGAAACAAATAGAACCATTTTCCTCTCCCCAGATAGTTTAGGATGAAGAGATCGGGGCGTTTGTAATAGCACATGACAACAATGCAATGGTCATCGTACACCATGTTTTTATTTAGCGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCtgtg                     >  2:503582/1‑251 (MQ=255)                                                                                                                                                                                                                                        gAGGGAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAAc                  <  1:95071/251‑1 (MQ=255)                                                                                                                                                                                                                                           ggAAACCTATTTTATTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACcaa               <  1:595098/251‑1 (MQ=255)                                                                                                                                                                                                                                                   aTTTTATTGAAAGGAAAATTCCCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAAGCGAACCATTTTCCTCTCCCCAGATATTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCAAATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAAAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACATTTgc        <  2:199667/250‑1 (MQ=255)                                                                                                                                                                                                                                                        cTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACAAACAATTGCGccccc  >  2:46421/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           ATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCC  >  NC_000913/3792586‑3793083

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