Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,444,905	C→T	G109S (GGT→AGT)	rpsE ←	30S ribosomal subunit protein S5

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,444,905	0	C	T	97.4%	132.7 / ‑4.0	39	G109S (GGT→AGT)	rpsE	30S ribosomal subunit protein S5
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (16/22);  total (16/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

CAGACGACCGATTGCACTGCGGGTTTGAGTAATTTTAATAGTCTTTGCCATGGTTTATTTCCCCAGAATTTCTTCAACGGATTTACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACCACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAATCATATTGCGACGGGCTTTTTCCATCGCTTTCTGGATCGCTGCTGGAACTTCACGCGCTTTACCGTAACCAAAACCAACGCGACCGTTACCATCGCCAACTACAGTCAGAGCTGTGAAGGAGA  >  NC_000913/3444672‑3445138                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                          cAGACGACCGATTGCACTGCGGGTTTGAGTAATTTTAATAGTCTTTGCCATGGTTTATTTCCCCAGAATTTCTTCAACGGATTTACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCgg                                                                                                                                                                                                                          <  2:341457/251‑1 (MQ=255)  aGACGACCGATTGCACTGCGGGTTTGAGTAATTTTAATAGTCTTTGCCATGGTTTATTTCCCCAGAATTTCTTCAACGGATTTACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGt                                                                                                                                                                                                                         <  1:303675/251‑1 (MQ=255)                    cGGGTTTGAGTAATTTTAATAGTCTTTGCCATGGTTTATTTCCCCAGAATTTCTTCAACGGATTTACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTg                                                                                                                                                                                                       <  2:60169/251‑1 (MQ=255)                             gTAATTTTAATAGTCTTTGCCATGGTTTATTTCCCCAGAATTTCTTCAACGGATTTACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTGCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATAc                                                                                                                                                                                              <  2:53496/251‑1 (MQ=255)                              tAATTTTAATAGTCTTTGCCATGGTTTATTTCCCCAGAATTTCTTCAACGGATTTACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGAAACCGGTACCTTCGGAAGCCGGCTGCATGAATAcg                                                                                                                                                                                             >  2:305895/1‑251 (MQ=255)                               aaTTTTAATAGTCTTTGCCATGGTTTATTTCCCCAGAATTTCTTCAACGGATTTACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATAcg                                                                                                                                                                                             <  1:346271/250‑1 (MQ=255)                                aTTTTAATAGTCTTTGCCATGGTTTATTTCCCCAGAATTTCTTCAACGGATTTACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATAcgcg                                                                                                                                                                                           <  2:252748/251‑1 (MQ=255)                                         gTCTTTGCCATGGTTTATTTCCCCAGAATTTCTTCAACGGATTTACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCgt                                                                                                                                                                                  <  1:47402/251‑1 (MQ=255)                                          tCTTTGCCATGGTTTCTTTCCCCAGAATTTCTTCAACGGATTTACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGCACCCgtg                                                                                                                                                                                 >  1:106618/1‑251 (MQ=255)                                          tCTTTGCCATGGTTTATTTCCCCAGAATTTCTTCAACGGATTTACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCgtg                                                                                                                                                                                 >  1:267587/1‑251 (MQ=255)                                                     gTTTATTTCCCCAGAATTTCTTCAACGGATTTACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCttt                                                                                                                                                                      <  1:207653/251‑1 (MQ=255)                                                        tATTTCCCCAGAATTTCTTCAACGGATTTACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTaa                                                                                                                                                                    <  1:169926/250‑1 (MQ=255)                                                               ccAGAATTTCTTCAACGGATTTACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAgggggt                                                                                                                                                            >  2:208786/1‑248 (MQ=255)                                                                     tttCTTCAACGGATTTACCACGCTTGGCCGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCCCTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATAcgc                                                                                                                                                                                            >  2:33559/1‑213 (MQ=255)                                                                     tttCTTCAACGGATTTACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATAcgc                                                                                                                                                                                            <  1:33559/213‑1 (MQ=255)                                                                               ggATTTACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTAt                                                                                                                                            <  2:323711/251‑1 (MQ=255)                                                                                    tACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATAcgcg                                                                                                                                                                                           >  1:69288/1‑199 (MQ=255)                                                                                    tACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATAcgcg                                                                                                                                                                                           <  2:69288/199‑1 (MQ=255)                                                                                              ggCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACCACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTaa                                                                                                                             <  1:309329/251‑1 (MQ=255)                                                                                               gCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTGTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAAt                                                                                                                            <  2:38738/251‑1 (MQ=255)                                                                                                  gCGACCATTTCTGGAGAATTCAGATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAATCat                                                                                                                         <  2:134224/251‑1 (MQ=255)                                                                                                      ccATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAATCATATTg                                                                                                                     <  1:81521/251‑1 (MQ=255)                                                                                                              ggAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAATCATATTGCGACGGGc                                                                                                             <  1:363835/251‑1 (MQ=255)                                                                                                                agaATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAATCATATTGCGACGGGCtt                                                                                                           >  2:100041/1‑251 (MQ=255)                                                                                                                                    cATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAATCATATTGCGACGGGCTTTTTCCATCGCTTTCTGGATc                                                                                       >  2:102116/1‑251 (MQ=255)                                                                                                                                     aTCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAATCATATTGCGACGGGCTTTTTCCATCGCTTTCTGGATCg                                                                                      <  1:83547/251‑1 (MQ=255)                                                                                                                                       cAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCTACTTCCATAACGGCGCGCAGTGCACTACCGGCGATTATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACTCCTTTAACAGGGTGTTGCAGAGTGCCTTTATTCAGCGCGACGTTAATCATATTGCGACTGGCTTTTGCCATCGCTTTCTGGATCgct                                                                                    <  2:306263/251‑1 (MQ=255)                                                                                                                                        aaTAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAATCATATTGCGACGGGCTTTTTCCATCGCTTTCTGGATCgctg                                                                                   >  1:278215/1‑251 (MQ=255)                                                                                                                                                cACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAATCATATTGCGACGGGCTTTTTCCATCGCTTTCTGGATCGCTGCTGGAACt                                                                           >  1:212153/1‑251 (MQ=255)                                                                                                                                                                      ggtggAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTATGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAATCATATTGCGACGGGCTTTTTCCATCGCTTTCTGGATCGCTGCTGGAACTTCACGCACTTTACCGTaaccaa                                                     <  2:255886/251‑1 (MQ=255)                                                                                                                                                                                tAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAATCATATTGCGACGGGCTTTTTCCATCGCTTTCTGGATCGCTGCTGGAACTTCACGCGCTTTACCGTAACCAAAACCAAcgcg                                           >  1:343533/1‑251 (MQ=255)                                                                                                                                                                                      ttAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAATCATATTGCGACGGGCTTTTTCCATCGCTTTCTGGATCGCTGCTGGAACTTCACGCGCGTTACCGTAACCAAAACCAACGCGACCGt                                      <  1:305895/250‑1 (MQ=255)                                                                                                                                                                                            aGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCAATACCGGCGATGATACAGGTACCTTAGGAAGCCGGCTGCATGAATCCGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGCCGATAATCATATTGTGACGGGCTTTTTCCATCGCTTTATGGATAGCTGCTGGAACTTCACGCGCTTTACCGTAACAAAAACCAACGCGACCGTTACCAt                                >  1:351990/1‑250 (MQ=255)                                                                                                                                                                                                 gTTATGCACCCCAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAATCATATTGCGACGGGCTTTTTCCATCGCTTTCTGGATCGCTGCTGGAACTTCACGCGCTTTACCGTAACCAAAACCAACGCGACCGTTACCATCGCCa                           >  1:117846/1‑250 (MQ=255)                                                                                                                                                                                                          cccAGCGACTTCCAGAACGGCGCGCATTGCACTACCTGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTACTCATATTGCGACGGGCTTTTTCCATCGCTTTCTGGATCGCTGCTGTAACTTCACGCGCTTTACCGTAACCAACACCAACGCGACCGTTCCCATCGCCAACTACCGTcc                 >  2:341836/1‑250 (MQ=255)                                                                                                                                                                                                          cccAGCGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAATCATAGTGCGACGGGCTTTTTCCATCGCTTTCTGGATCGCTGCTGGAACTTCACGCGCTTTACCGTAACCAAAACCAACGCGACCGTTACCATCGCCAACTACAGTCa                 <  2:322505/251‑1 (MQ=255)                                                                                                                                                                                                               cGACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAATCATATTGCGACGGGCTTTTTCCATCGCTTTCTGGATCGCTGCTGGAACTTCACGCGCTTTACCGTAACCAAAACCAACGCGACCGTTACCATCGCCAACTACAGTCAGAGCt            >  1:107667/1‑251 (MQ=255)                                                                                                                                                                                                                gACTTCCAGAACGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAATCATATTGCGACGGGCTTTTTCCATCGCTTTCTGGATCGCTGCTGGAACTTCACGCGCTTTACCGTAACCAAAACCAACGCGACCGTTACCATCGCCAACTACAGTCAGAGCtg           >  1:183186/1‑251 (MQ=255)                                                                                                                                                                                                                         aaCGGCGCGCATTGCACTACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAATCATATTGCGACGGGCTTTTTCCATCGCTTTCTGGATCGCTGCTGGAACTTCACGCGCTTTACCGTAACCAAAACCAACGCGACCGTTACCATCGCCAACTACAGTCAGAGCTGTGAAGgaga  <  2:69154/251‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                          CAGACGACCGATTGCACTGCGGGTTTGAGTAATTTTAATAGTCTTTGCCATGGTTTATTTCCCCAGAATTTCTTCAACGGATTTACCACGCTTGGCAGCGACCATTTCTGGAGAATTCATATTTTCCAGGCCATCAATAGTTGCACGAACCACGTTGATCGGGTTGGTGGAACCATAGGCTTTAGCCAGAACGTTATGAACCCCAGCGACTTCCAGAACGGCGCGCATTGCACCACCGGCGATGATACCGGTACCTTCGGAAGCCGGCTGCATGAATACGCGAGAACCCGTGTGAACACCTTTAACAGGGTGTTGCAGAGTGCCGTTATTCAGCGCGACGTTAATCATATTGCGACGGGCTTTTTCCATCGCTTTCTGGATCGCTGCTGGAACTTCACGCGCTTTACCGTAACCAAAACCAACGCGACCGTTACCATCGCCAACTACAGTCAGAGCTGTGAAGGAGA  >  NC_000913/3444672‑3445138

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