Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,143,323	C→T	D349N (GAC→AAC)	rne ←	ribonuclease E

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,143,323	0	C	T	97.8%	158.0 / ‑4.2	46	D349N (GAC→AAC)	rne	ribonuclease E
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (22/23);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

CTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTCGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGATCTGTGACTCGATCTGGTAGTGGCT  >  NC_000913/1143080‑1143569                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       cgtcttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGTCAGTACCAGAACAACGCGTACAACCGTGATGACGGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAg                                                                                                                                                                                                                                                 <  2:373246/249‑1 (MQ=255)  ttcttcAATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCCCC‑GCTTCACGCCGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCCTGTTGATgtagg                                                                                                                                                                                                                                                >  2:345159/1‑247 (MQ=255)        gATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGCTGACGTCGATa                                                                                                                                                                                                                                           >  1:141534/1‑250 (MQ=255)           cAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACaa                                                                                                                                                                                                                                       <  2:8487/251‑1 (MQ=255)           cAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTCGATGAAGTCGATAACa                                                                                                                                                                                                                                        >  1:313421/1‑250 (MQ=255)                   gAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGCCGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGAGGCGTACTGGCGTCATGTTGATGAAGGCGATAACAATCCGgccg                                                                                                                                                                                                                               >  2:10152/1‑251 (MQ=255)                                            gTCACGCACGGGGCCAGGACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCGGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAAGCAGGCCGCCGAGGTCACGCAGGCGCAGCTGAc                                                                                                                                                                                                      <  2:201595/251‑1 (MQ=255)                                                   ccgGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAAt                                                                                                                                                                                               <  2:379724/250‑1 (MQ=255)                                                     ggTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAAtct                                                                                                                                                                                             >  1:243587/1‑251 (MQ=255)                                                          cAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCg                                                                                                                                                                                        <  2:434140/251‑1 (MQ=255)                                                                        cgcgGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTg                                                                                                                                                                                                        >  1:176854/1‑221 (MQ=255)                                                                        cgcgGACAAACGTGATGACTGGATTCACCCAGGGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTg                                                                                                                                                                                                        <  2:176854/221‑1 (MQ=255)                                                                           ggACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAgt                                                                                                                                                                       <  1:307843/251‑1 (MQ=255)                                                                              cAAACGTGATGATTGGATTCCCCCATGGATGG‑CTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTa                                                                                                                                                                   <  1:139312/251‑1 (MQ=255)                                                                               aaaCGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTa                                                                                                                                                                   <  1:10152/251‑1 (MQ=255)                                                                                aaCGTGATGCCTGGATTTAAACAGAGAAGGGCCCAGGCGCTGACGGGCCATTTCCAGCAGGCCAAAGCGAGCAATCTGGCTGATTTGAATACGCGCCCGGACCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGCGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCCGCTTCGAGGTTAGTGTTa                                                                                                                                                                   >  2:458224/1‑250 (MQ=255)                                                                                          cTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCa                                                                                                                                                                                                                                    >  2:223775/1‑175 (MQ=255)                                                                                          cTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCa                                                                                                                                                                                                                                    <  1:223775/175‑1 (MQ=255)                                                                                           tGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGCGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTttct                                                                                                                                                       >  2:336964/1‑251 (MQ=255)                                                                                             gATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTttcttc                                                                                                                                                     >  1:268752/1‑251 (MQ=255)                                                                                                      aGTGATGGGCTCCGGCGCTGTCGGGACATGTCCAGCAGTCCAAAGCGATAAATATGGCTGATGTTAATACGCCCACGGTCCTGACGCACT‑GCTTCACGCATACGGTTTTCTACCGCACGCTGGTGGCGTACTGGTGGCATGTTGATGAAGTCGATAGCAAGCAGTCCGCCTAGGGCACGCAGGCGCAGCTGAGGTGCAATCTCATCGGCGGCTTCGAGGTGAGTGTTAAACGCGGTTTCTTCGAGATcgcc                                                                                                                                            <  2:201646/251‑1 (MQ=255)                                                                                                         gATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTTTCTTCGATATcgccgcc                                                                                                                                         <  1:336964/251‑1 (MQ=255)                                                                                                         gATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTTTCTTCGATATcgccgcc                                                                                                                                         <  1:406954/251‑1 (MQ=255)                                                                                                                     cacTGACGTGACATTTCCAGCAGGCCAAAGCGAGATATATGGCTGATTTGGATACGCGCACGGGCCTCAAG‑ACCAGCTTAAAGCAGCCGGTTTTCTACCGCACGCTGGTGGCGTACTGGCATCATGTTGCTGGAATCGATAACAATCAGGCCGCCGAAGTTTCGCAGGCGAAGCTGAAAAGCAATCGTATCGACCGCTACGGGGTTATGGTTACACGCGGTTTCTTAGCTATCGCCGACGCGGGAcgcgcg                                                                                                                             >  2:432768/3‑251 (MQ=255)                                                                                                                              gACATTTCCAGCAGGCCAATGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTGTACGCAGACTGTTTTCTACCGCACGCTGGGGGCGGACTGTCGTCAGGTTGATGAAGGCAATAACAATCAGGCCGCCGGGGTCACGCTTGCGCAGCTGACGAGCAATCGCCGCGGCAGCTTCGTGGTTAGGGTTAAACGCGGTTTCGTCGATAGCGCCTCCTCGGCTCGCACGTGCGGattt                                                                                                                    <  2:38195/251‑4 (MQ=255)                                                                                                                               aCATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCTCGCGTGCGGAGTTg                                                                                                                   <  2:141534/251‑1 (MQ=255)                                                                                                                                     ccAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGAGCGTGCGGAGTTGATGTCg                                                                                                             >  1:103524/1‑251 (MQ=255)                                                                                                                                       agcagGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGAt                                                                                                           >  1:379126/1‑251 (MQ=255)                                                                                                                                                 aGCGAGAAATATTGCTGATTTGAATGCGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCTATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAAc                                                                                                 <  1:106641/251‑1 (MQ=255)                                                                                                                                                      gaAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCTTc                                                                                            <  1:205721/251‑1 (MQ=255)                                                                                                                                                       aaaTATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCTTCg                                                                                           <  1:345159/251‑1 (MQ=255)                                                                                                                                                                 gATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCTTCGGTGCTGTCg                                                                                  >  2:158470/1‑250 (MQ=255)                                                                                                                                                                     tGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATaa                                                                              <  1:27827/250‑1 (MQ=255)                                                                                                                                                                      gAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGt                                                                                                                                                           >  1:288191/1‑172 (MQ=255)                                                                                                                                                                      gAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGt                                                                                                                                                           <  2:288191/172‑1 (MQ=255)                                                                                                                                                                        aTACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAAt                                                                          >  1:322946/1‑251 (MQ=255)                                                                                                                                                                             cgcACGGACCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAac                                                                     >  1:329739/1‑251 (MQ=255)                                                                                                                                                                               cACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGCGGTTAGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGCTGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAaccc                                                                   >  2:255667/1‑250 (MQ=255)                                                                                                                                                                                  ggaccTGACGCACC‑GCTTAACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCAGAGGTCACGCAGGCGCAGCTGACGCGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTCGACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGCTGGCCGTTAGCGCTTCGGTGCTGTCGATAACAATGGAacaacc                                                                 >  2:101157/4‑250 (MQ=255)                                                                                                                                                                                  ggTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAaccacc                                                                 >  1:78864/1‑250 (MQ=255)                                                                                                                                                                                             aCC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTCGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAACCACCAGACGGCCGACg                                                     >  2:163980/1‑251 (MQ=255)                                                                                                                                                                                                                 tttCTACCGCACGCTGGTGGCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAAATTCACGCTGGAAGGCg                                  >  1:325397/1‑251 (MQ=255)                                                                                                                                                                                                                      acagCACGCTGGTGGCGGCCTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGTCGAGCAAGCTCATCGGCAGCTTCGAGTTTAGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGTGCGCGCGTGCGGAGTTGATGTGGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGGACCACCAGACGGCATATGAACTTCACGCTGGAAGGCGGAc                               <  2:442446/246‑1 (MQ=255)                                                                                                                                                                                                                                   ggCGTACTGGCGTCATGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGGTGATATCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGATCTGTGACTc                 <  1:372335/250‑1 (MQ=255)                                                                                                                                                                                                                                   ggCGTACTGGCGTCATGTTGATGAAGGCGATAACAATCAGGCCGCCGATGTCACGCAGGCGCAGCTGACGAGCAATCGCATCGTCAGCTTCGAGGTGAGTGTTAAACGCGTTTTCTTCAATATTGCCGCCGCGGGTCGCGCGTGCGGAGTTGGTGTCGATGGCCGTTAACGTTTCTGTGCTGTCGATAACAATGGAACCACCAGACGGCAGATGAACTTCACGCTGGAAGGCGGACGCGATCGGTGACTCg                <  2:267816/251‑1 (MQ=255)                                                                                                                                                                                                                                                 aTGTTGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGATCTGTGACTCGATCTGGTAGTGGCt  >  2:89877/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       CTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCACGGTCCTGACGCACC‑GCTTCACGCAGACGGTTTTCTACCGCACGCTGGTGGCGTACTGGCGTCATGTCGATGAAGTCGATAACAATCAGGCCGCCGAGGTCACGCAGGCGCAGCTGACGAGCAATCTCATCGGCAGCTTCGAGGTTAGTGTTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGATCTGTGACTCGATCTGGTAGTGGCT  >  NC_000913/1143080‑1143569

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