Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,967,770	C→T	intergenic (+94/‑146)	rho → / → rfe	transcription termination factor Rho/UDP‑N‑acetylglucosamine‑‑undecaprenyl‑phosphate N‑acetylglucosaminephosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,967,770	0	C	T	97.7%	143.7 / NA	44	intergenic (+94/‑146)	rho/rfe	transcription termination factor Rho/UDP‑N‑acetylglucosamine‑‑undecaprenyl‑phosphate N‑acetylglucosaminephosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (16/27);  minor base G (0/1);  total (16/28)
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.99e-01

CTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGG  >  NC_000913/3967524‑3967994                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                 ctgctCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCAGAAGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGATCATAAATTTGACTTATGCCAAAAACACCACGTGTTTAAGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGGTGTTAGGAATATTCTTgg                                                                                                                                                                                                                               >  2:390127/1‑250 (MQ=255)  tgctCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGa                                                                                                                                                                                                                              >  2:77433/1‑250 (MQ=255)   gctCACGACTCAGGAAGAACTTCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTTCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAAt                                                                                                                                                                                                                            <  2:436907/251‑1 (MQ=255)   gctCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGaa                                                                                                                                                                                                                             <  1:9780/250‑1 (MQ=255)      cACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTg                                                                                                                                                                                                                         <  1:40820/251‑1 (MQ=255)      cACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTg                                                                                                                                                                                                                         <  1:77433/251‑1 (MQ=255)             cAGGGAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGTGCTGGCGCTGTTAGGAAAATTCTTGGAATTTGCTGGCa                                                                                                                                                                                                                   <  2:288557/250‑1 (MQ=255)             cAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCAt                                                                                                                                                                                                                  <  1:56090/251‑1 (MQ=255)                gaagaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGtt                                                                                                                                                                                                               <  1:250285/251‑1 (MQ=255)                  agaaCGGCAGAAAATGTTGATCCTGCGCAAAAGCAAGCACCCGATGGGCGAAAGCGATGCAATGGAATTCCTCATTAAGAAACTGGCAAGGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGGGTTTACGTGGCGTATTGCTTTTATAGCTTTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCAAGttct                                                                                                                                                                                                             <  2:119056/251‑3 (MQ=255)                               aTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCatat                                                                                                                                                                                                >  1:421753/1‑251 (MQ=255)                                     aTCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATGCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATg                                                                                                                                                                                          <  1:150212/251‑1 (MQ=255)                                          gcgcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTaa                                                                                                                                                                                     <  2:147330/251‑1 (MQ=255)                                              aaaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAAtttt                                                                                                                                                                                 <  1:99588/251‑1 (MQ=255)                                                  tCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGc                                                                                                                                                                              >  1:88233/1‑250 (MQ=255)                                                     ttCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGcac                                                                                                                                                                            >  1:435245/1‑249 (MQ=255)                                                       cACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACTCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGTAATTTGCTGGCATGTTATGCAATTTGCATATCaaa                                                                                                                                                                                            <  2:493526/231‑1 (MQ=255)                                                       cACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCaaa                                                                                                                                                                                            >  1:493526/1‑231 (MQ=255)                                                                                       ccTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTg                                                                                                                                         <  1:344991/250‑1 (MQ=255)                                                                                              aataaACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAg                                                                                                                                  <  1:454536/250‑1 (MQ=255)                                                                                                      ggCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGt                                                                                                                          <  1:174238/250‑1 (MQ=255)                                                                                                         aaTGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGtt                                                                                                                      <  2:374665/251‑1 (MQ=255)                                                                                                                 agaccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCt                                                                                                              <  1:366587/251‑1 (MQ=255)                                                                                                                       aTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCtaat                                                                                                        >  1:481553/1‑251 (MQ=255)                                                                                                                         gtcGATTTATTCGAAATGATGAAACGCTCAGAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATTTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCtaata                                                                                                       <  2:395269/248‑1 (MQ=255)                                                                                                                            gATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAAttt                                                                                                   >  1:225096/1‑251 (MQ=255)                                                                                                                                 cttcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTACTTTTTGCACAGGACTGGTGGGTTTGGCACGGACTTTcc                                                                                                                                               >  2:37820/1‑202 (MQ=255)                                                                                                                                 cttcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTcc                                                                                                                                               <  1:37820/202‑1 (MQ=255)                                                                                                                                     gAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGACCGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTg                                                                                             >  2:320187/1‑248 (MQ=255)                                                                                                                                      aaaTGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATTTCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTAGGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGGTTGGAACGGACTTTCCCTTCGTATTAAAGTTCTTCGTGGTTATACTTCTGCTAATTATTTTCTCTgtg                                                                                           >  2:361273/1‑247 (MQ=255)                                                                                                                                                        aaaTTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTAc                                                                       <  1:396013/251‑1 (MQ=255)                                                                                                                                                            ttGTCTTATGCCAAACACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAAtt                                                                                                                                                                                   >  2:426275/1‑139 (MQ=255)                                                                                                                                                            ttGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAAtt                                                                                                                                                                                   <  1:426275/139‑1 (MQ=255)                                                                                                                                                                   atgtcAAAAACGCCACGTGTTTACGTGTCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATGTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCAGCTAATAAGATTCTCTGAGAGCATTCATTGTGAATTTCCGGACAGTGAGt                                                            <  2:114432/247‑1 (MQ=255)                                                                                                                                                                   aTGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCATAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTTGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGt                                                            <  1:58733/251‑1 (MQ=255)                                                                                                                                                                    tGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATg                                                                                                                                                                                                                 >  1:241761/1‑101 (MQ=255)                                                                                                                                                                    tGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAAGTTGCTGGCATg                                                                                                                                                                                                                 <  2:241761/101‑1 (MQ=255)                                                                                                                                                                       cAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTg                                                        <  1:418489/251‑1 (MQ=255)                                                                                                                                                                          aaaCGCCACGTGATTACGGGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGCAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGAGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAAAGCATGCATTGTGAATTTACTGGCAGTGAGTACTgtt                                                      >  2:358760/1‑248 (MQ=255)                                                                                                                                                                                                                    aTCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTg           >  1:339476/1‑251 (MQ=255)                                                                                                                                                                                                                         aaTGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTa                                    >  1:59608/1‑221 (MQ=255)                                                                                                                                                                                                                         aaTGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTa                                    <  2:59608/221‑1 (MQ=255)                                                                                                                                                                                                                         aaTGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGt      <  2:433709/251‑1 (MQ=255)                                                                                                                                                                                                                             ccccgttGGCGATGTTAGGAAAATTCGTTGAATTTGCTGGCATGTTATTCAATTTGCATATCAAATGGTTAATTTTTGCACTTGACTGGTGGGTATGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTTTACTTCTGCTAATAATTTTCTCTGAGATCATGCATTGTGAATTTACTGACAGTGAGTTCTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTCTGCCCGTAAgt  <  2:429593/245‑2 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                 CTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGG  >  NC_000913/3967524‑3967994

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