Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	46,976	A→G	L390L (CTA→CTG)	yaaU →	putative transporter YaaU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	46,976	0	A	G	98.4%	191.5 / NA	61	L390L (CTA→CTG)	yaaU	putative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (28/32);  minor base T (1/0);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTG‑GCGGGCCGGGTAATTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGTCGTATTGGCACCATTGTTTCGACCTGGGCACTACCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAGACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGG  >  NC_000913/46729‑47214                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                               ggCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGGTCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTTTCTG‑GCGGGCCGGGTAATTTGCAGTG‑GCTCTATACTTATGAACTCTTTCCGACAGATATCCCCGCCTCTGCCCCGGGCGTGATTATGTCCTTAAGCCGTATTGGCCCCATTGGTTTGACCTGGGCACTTCCg                                                                                                                                                                                                                                             >  1:5035/1‑251 (MQ=255)      ccATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTG‑GCGGGCCGGGTAATTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCtt                                                                                                                                                                                                                                        <  1:367092/251‑1 (MQ=255)                 tGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGGGATGGCCTTTGCGGTGTATGCCTTTTTCTCTG‑GCGGGCCGGGGAATTTGCAGTG‑GGTCTATCCTTATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGGGGGCGTGATTATTTCCTTAAGCCGTATTTGCCCCATTGTTTCGTCCCTGGCACTGCCGATCTTTTTTAATAAtt                                                                                                                                                                                                                             >  1:114375/1‑251 (MQ=255)                          gggcaggatggcggtggcgctggcgGTTTTGGGGCTGATCACGGGTATGGGGATCTGGCTGGTGGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTG‑GCGGGCCGGGTAATTTGCAGTG‑TCTCTATCCTAATGAACTCTTCCTGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCGTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACTGTATc                                                                                                                                                                                                                     <  2:436444/236‑1 (MQ=255)                                gatgaCGCGGGCGTTGGTGGTTTTGGGGCTGAAGCCGGATTTGGGGATTTGGCTGGTAGTGATGGCCTTTGCGTTGTATGCCTTTTTCTCTG‑GCGGGCCGGGTATTTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAaca                                                                                                                                                                                                              <  2:465329/251‑1 (MQ=255)                                     cgctggcgctggCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTG‑GCGGGCCGGGTAATTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATg                                                                                                                                                                                                         <  1:31775/251‑1 (MQ=255)                                               ggcgGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTG‑GCGGGCCGGGTAATTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATgggg                                                                                                                                                                                                <  1:515448/250‑1 (MQ=255)                                                  gtgtttGGGGCTTATTCCGGATATGGGGTTCTGGTTGGGTGTTATGGCCTTTGCTGGGTAGGCCTTTTTCTCTG‑GCGGGCCGGGTAATTTGCAGTG‑GCTCTATCCGAATGAACTCTTCCCTACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGTGGGCgg                                                                                                                                                                                            <  2:404272/248‑1 (MQ=255)                                                        ggggCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTG‑GCGGGCCGGGTAATTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTAtct                                                                                                                                                                                      >  1:176086/1‑251 (MQ=255)                                                           gctgttCCCGGATAGGGGGATCTGGTTGGTAGTGATGGCCGTGGCGGTGTATGCCTTTTTCTCTG‑GCGGGCCGGGTAATTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCTTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGc                                                                                                                                                                                   <  1:632369/246‑1 (MQ=255)                                                           gctgaTCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTG‑GCGGGCCGGGTAATTTGCAGTG‑GCTCTATCCTAATGCACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCCATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGc                                                                                                                                                                                   >  1:287504/5‑251 (MQ=255)                                                               agttctgatctgggtgtCTGGATGGTAGTGGTGGTCTTTGTGGTGTTTGCCCTTTTCTCTG‑GCTGCCCTTATAATTTGCCGTG‑TCGCTATCTTCACGAGCTTTTCCCGTCTAATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTCTTGGTTCTATTGTTTCGACATGGGCACTGCCGATCTTTATGAACAAGTACGATATCAGTAGCACGATGCTAATGTCTGCGGGTAtctaactgt                                                                                                                                                                                <  2:502694/234‑7 (MQ=255)                                                                                   gcgggttgttttggtcttggCGGTGTATGCCTTTTTC‑GTGTGCGGGTCGGGTAATTTGCTGTG‑GATCGATCATAAGGAACTCTTCCAGAAAGATTCCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAATGCGTATTGGCCCCATTGTTTCGACCTGGGCGCTGCCGATCTTTATCAATAATTGCTGTATCAGTATCACGATGCTAATGGGGGAGGGTATCTCGCTGTTTGGCTTGTTGATTTCTGTAg                                                                                                                                                           <  2:213120/232‑1 (MQ=255)                                                                                      ggTAGTGATGGCCTTTTCGGTGTATGCCTTTTCCTCTG‑GCGGGCCGGGTAAGTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGt                                                                                                                                                                                          <  1:494381/217‑1 (MQ=255)                                                                                         aGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTG‑GCGGGCCGGGTAATTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTg                                                                                                                                                     >  1:11697/1‑251 (MQ=255)                                                                                            gctGGCCATTGCGGTGTATGCCTTTTTCTCTG‑GCGGGACGGGTAATTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATAACCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGGATTGACTCCATTGTTTCGACCTGGGCACTGCCGATCCTTATCAATACTTACGGTATCAGTAACACGGTGCTAATGGGGGAggg                                                                                                                                                                                           >  2:12168/3‑210 (MQ=255)                                                                                            gATGGCCTTTGCGGTGTATGCCTTTTTCTCTG‑GCGGGCCGGGTAATTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCggg                                                                                                                                                                                           <  1:12168/210‑1 (MQ=255)                                                                                               gtcgctgccggtggagtccTTGTTCTCTG‑GCGGGCGGGGGAATTGGCA‑TGTGCTCTATCTCAAGGAACTCTTCCCGACAGATATCCGAGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCATCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTCCGGTAGCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCTTAGCGTTTGcgccc                                                                                                                                                <  2:90872/233‑5 (MQ=255)                                                                                                     tGCGGTGTATGCCTTTTTCTCTG‑GCGGGCCGGGTAATTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTc                                                                                                                                         >  1:678129/1‑251 (MQ=255)                                                                                                           gtATGCCTTTTTCTCTG‑GCGGGCCGGGTAATTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAgggg                                                                                                                                   >  1:304951/1‑250 (MQ=255)                                                                                                              tGCCTTTTTCTCTG‑GCGGGCCGGGTAATTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGGATCTCGCTGTTTGGCTTGTTGATTTCCGTAACGTTTGCCCCGGGGGCTCGAGGGATgg                                                                                                                                >  1:142978/1‑250 (MQ=255)                                                                                                               gCCTTTTTCGCTG‑GCGGGCCTGGTAATTTGCAGTG‑GCTCGATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGt                                                                                                                                <  1:632653/250‑1 (MQ=255)                                                                                                                              cGGGCCGGGTAATTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAAccc                                                                                                                 >  2:115025/1‑250 (MQ=255)                                                                                                                               gggCCGGGTAATTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCCCTGGCGCAAAcc                                                                                                                  >  1:670789/1‑249 (MQ=255)                                                                                                                                   cGGGTAATTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCt                                                                                                                                                                                                       <  1:584956/160‑1 (MQ=255)                                                                                                                                   cGGGTAATTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGCTTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGGTGTCACTGGCGCAAACCAGCact                                                                                                            >  2:268390/1‑249 (MQ=255)                                                                                                                                   cGGGTAATTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGAGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCt                                                                                                                                                                                                       >  2:584956/1‑160 (MQ=255)                                                                                                                                         aTTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCtt                                                                                                                                                                                                                                        <  2:587990/121‑1 (MQ=255)                                                                                                                                         aTTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCtt                                                                                                                                                                                                                                        >  1:587990/1‑121 (MQ=255)                                                                                                                                                gTG‑GCGCTATCATAATGATCTTGTCCCGACAGTTATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGCCCTGGGCACTGCTGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGCCGTTCcgcg                                                                                               <  2:243067/251‑1 (MQ=255)                                                                                                                                                 tG‑GCTCTATCCTAATGAACTCTTCCCGACAGAGATCCGAGCCTCTGCCGTGGGCGGGATTATGTCCTTAAGCCGGATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGAATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTTCCCCGGAGACTCGAAGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCgg                                                                                              >  1:120723/1‑251 (MQ=255)                                                                                                                                                      tctATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAatat                                                                                                          >  1:104138/1‑235 (MQ=255)                                                                                                                                                      tctATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAatat                                                                                                          <  2:104138/235‑1 (MQ=255)                                                                                                                                                         aTCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCagag                                                                                        <  1:22015/250‑1 (MQ=255)                                                                                                                                                         aTCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAgaga                                                                                       <  1:445783/251‑1 (MQ=255)                                                                                                                                                              aataatatcttcccaccagagttCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGGATTGGCACCATTGTGTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATggg                                                                                  <  2:520398/229‑1 (MQ=255)                                                                                                                                                              aaTGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTTGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATgg                                                                                   <  1:614875/250‑1 (MQ=255)                                                                                                                                                                tGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATg                                                                                    <  1:633316/247‑1 (MQ=255)                                                                                                                                                                tGAACTCTTCCCAAAAGATATCCGCGCCTCTGCCGTGGGCGCGGTTATGTCCTTCAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGTTTTCCGTAACGTTTGACACGGAGAATCGAGGGATGTAACTGGCGCAAACCAGCACTATGACGACGCGCGGGCAGGGGATgggg                                                                                 >  2:633316/1‑250 (MQ=255)                                                                                                                                                                          ccGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCCGCCATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTc                                                                       >  2:506845/1‑250 (MQ=255)                                                                                                                                                                                  tatCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGATAATGGGGGCGGGt                                                                                                                                                                                          <  1:298522/127‑1 (MQ=255)                                                                                                                                                                                  tatCCGCGCATCTGCCGTGGGCGAGATTATGTCATTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCTATAATTACGGTATCAGTAACCCGATGCTAATGGGGGCGGGt                                                                                                                                                                                          >  2:298522/1‑127 (MQ=255)                                                                                                                                                                                   atCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTtctct                                                             >  1:487794/1‑251 (MQ=255)                                                                                                                                                                                                  gTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAat                                              >  1:239336/1‑251 (MQ=255)                                                                                                                                                                                                  gTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAat                                              >  2:483332/1‑251 (MQ=255)                                                                                                                                                                                                   tGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGCGACTCGAgg                                                                                                                                     >  2:344893/1‑163 (MQ=255)                                                                                                                                                                                                   tGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAgg                                                                                                                                     <  1:344893/163‑1 (MQ=255)                                                                                                                                                                                                    gggCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAata                                             <  1:569073/250‑1 (MQ=255)                                                                                                                                                                                                     ggCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTCTCAATAATTACGGAATCAGTAACACGATGCTAATGTGGGCGGGGATCTCGCAGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGTGTCTCGAGGGATGTCACTGGCGAAACCCAGAATTATGAACATCCGAGGGCTGGGATTGGGGGAATTTTTTCAGGTTTCTCCCTTTTCTTAAtacataat                                           >  2:619842/1‑244 (MQ=255)                                                                                                                                                                                                                     tAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCcg                            <  1:483332/250‑1 (MQ=255)                                                                                                                                                                                                                      aagccttATTGGCCCCTTTGTTTCGACCTGGGCACTGACGTTCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCcgc                           <  2:209545/244‑1 (MQ=255)                                                                                                                                                                                                                            tATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACAGTATCAGTAACACGATGCTAATGGGGGAGGGTATCTCGCTGTTTGACTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACACGAGGGATGTCACTTGCGCAAACCAGCAATATGACGTTCCGCGGGCGGAGAAGTGGGTAAATTGTTCAGATTTCTCTCATTTCTGAATCAATATTATTGACAATAAGCCGCGTGAat                     >  2:197842/1‑250 (MQ=255)                                                                                                                                                                                                                              ttGGCACCACTGTGTCGACCTGGGCTCTGCCGATCTTTATCAATAATGACTGTATCAGGAACACGATGCTAATGGGGGTGGGTATCTCGCGGTTTGTCTTGTTGATTTCCGTAGCGGTTGCCCCGGAGTCTCGAGGGATTTCACTGGCGCATACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAAGTGTTCATATTGCTCTCTTTTCTGAATCAATATTATTGACTACAAGCCGCTTGAatag                   <  2:287504/250‑2 (MQ=255)                                                                                                                                                                                                                                 gCACCATTGTTTCGACCTGGTCACTGCCGATCTTTATCAATAATTACTGTACCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGa               <  1:368067/251‑1 (MQ=255)                                                                                                                                                                                                                                   cccATTGCTTCGACCTGGTCACTGCCGATCTTTATCAATAAGTACGGTAACAGTAACACGATGCGAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGTGCAGAGAATGGTGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACt             <  2:245293/250‑1 (MQ=255)                                                                                                                                                                                                                                   aCCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATAAGACt             >  1:145634/1‑251 (MQ=255)                                                                                                                                                                                                                                      aTTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTGCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTaca          <  1:89972/251‑1 (MQ=255)                                                                                                                                                                                                                                         gTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACtt       >  1:661040/1‑251 (MQ=255)                                                                                                                                                                                                                                         gTTTCGACCTGAGCACTGACGATCTTTATCACTAATTACGGTATCAGTACCACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGATGATTTCCGTAACGTTTGCCCCGGAGACTCGAGGGATGTCATTGGAGCAAACAAGCAATATGACGATCCGCGGGCAGAGAATGGGGTCAATTGTTCAGATTTCTCTATTTTCTGAATCAATTTTTTTGACTATAGGCCGCGTCACTATATTACtcaactt       >  2:276602/1‑245 (MQ=255)                                                                                                                                                                                                                                             cGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTtgtg   <  1:506845/251‑1 (MQ=255)                                                                                                                                                                                                                                              gACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTAAAAGCCGCGTGAATATATGACTACACTTTGTgg  >  2:554708/1‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                               GGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTG‑GCGGGCCGGGTAATTTGCAGTG‑GCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGTCGTATTGGCACCATTGTTTCGACCTGGGCACTACCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAGACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGG  >  NC_000913/46729‑47214

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