Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,402,938	G→A	A301T (GCT→ACT)	hflK →	regulator of FtsH protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,402,938	0	G	A	95.7%	131.5 / ‑1.6	47	A301T (GCT→ACT)	hflK	regulator of FtsH protease
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (17/27);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01

CGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCT  >  NC_000913/4402710‑4403173                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            cGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCtt                                                                                                                                                                                                                       <  1:59690/251‑1 (MQ=255)        aTGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGTCTGCTCGTACGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTaaa                                                                                                                                                                                                                <  1:378859/251‑1 (MQ=255)                 ggtgctcgctgtgggccGTCAACTCGCAGGTTGCTGGGCCGCCGGAAGAAGTAAGAGCGGCGTTTGACGATGCGATAGCCGCGCGTGAAAACGAATAGCAATACATGCGTGAAGCAGAAGCCTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCGTGAAGGGGCGCGTGCGTACATGGCCCAGACCATCCGTGAAACTCAGGGTGAAGTGGCGCGCGTTGCTAAACTTCTGCCg                                                                                                                                                                                                       <  2:66067/235‑1 (MQ=255)                  tATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGACGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAACCGAACAGCAATACATTCGTGCAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTACCCTCGAAGAGACGCGTGCGTACAAGGCCAAGCCCCTCCTGGCACCTCAGGGCGAAGTGGCGCGCTTTGCTAACATTCTGCCgg                                                                                                                                                                                                      >  2:172122/1‑251 (MQ=255)                       cGCTGCTGGAGGTACACTTCCAGGCTGCTCTGCCGCCGGAAGAAGTAAAGGCGGCGTTTGACCATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGGATACCAACGAAGTTCAGGCGCTTGCGAACGGTTAGGCGCAACGTCTCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGAGCATCCGGGAAACTCAGGGTGAAGTGGCGCTCTTTACTTATCTGCTGCCGGa                                                                                                                                                                                                     <  2:223402/247‑1 (MQ=255)                        gctgctGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAatat                                                                                                                                                                                                 >  1:257540/1‑250 (MQ=255)                        ccggctgGACGGAAACTTCCAGGCTGCTCGATCGCTGAAAGAAGGAAAAGCGGCGTTTGTCGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACGCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGACGGAAtata                                                                                                                                                                                                <  2:317621/248‑1 (MQ=255)                           gctgGACGTCAACGTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAg                                                                                                                                                                                             <  1:269956/251‑1 (MQ=255)                               gACGTCAACTTCCAGGCTGCTCGTCAGCCGGCAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACAAACAAAGATCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGCGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCCGGGCGAAGTGGCGCGCTTTGCTAAACTTCTACCGAAATATAAAGccac                                                                                                                                                                                         >  2:159146/1‑249 (MQ=255)                                    catcTTCCAGGCTGCTCGGCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCgg                                                                                                                                                                                    <  1:97341/248‑1 (MQ=255)                                         tCCAGGCTGCTCGGCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACAGTTGTGAAGCAGATGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCTCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAAt                                                                                                                                                                                <  2:272571/250‑1 (MQ=255)                                           cAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTAc                                                                                                                                                                             <  1:347423/251‑1 (MQ=255)                                                        cgcctGAAGAAGTAAAAGCGGCGTATGACGATGCGATTGCCTCGCGTGAAAACGAACAGCAATACAATCGTTAAGCAGAAGCGGATACCAACGAAGTTCAGCCGCTTGCTAACGTTCAGGCGCAACTTATCCTCGAAGAGGCGCGTGCGTACAAGTCCCAGACCATCCTGGAAACTCAGTGTGAAGTGGCGCGCTTTGCTAAACTGCTGCCGTAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGtatg                                                                                                                                                                <  2:128730/251‑4 (MQ=255)                                                          ccgGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGAAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCCGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCAACCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATACAGCCGCGCCGGAAATTACTCACGAGCGACTgg                                                                                                                                                               >  2:329368/1‑249 (MQ=255)                                                              aagaagTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATc                                                                                                                                                          <  1:172122/251‑1 (MQ=255)                                                               agaagTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCg                                                                                                                                                         >  1:219239/1‑251 (MQ=255)                                                                 aagTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCg                                                                                                                                                         >  1:213200/1‑249 (MQ=255)                                                                   gTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACCTTCGTGAAGCAGAAGCGTATACCCACGAAGTTCAGCCGCGTGCCAACGGTCAGGCGCAACGTATCCTCGAACAGGCGCGCGCGTACAAGGCCCAGACCATCATGGAAACTCAGGGTGAAGCGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCAGGAAATTACTCGCGAGCGTCTGTCAATCGAGAc                                                                                                                                                     >  1:31389/1‑251 (MQ=255)                                                                           ggcgTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGTTGAAGTGGCGCGGTTTGCTAAACTTCTGCAGGAATATAAAGCCGCGCCGGAAATTACTCGCTATCGTCTGTATATCGAGACGATGGaaa                                                                                                                                             <  1:371905/251‑1 (MQ=255)                                                                               tttGACGATGCGCTTGACGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGCAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTg                                                                                                                                                                                                                     >  2:331245/1‑175 (MQ=255)                                                                               tttGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTg                                                                                                                                                                                                                     <  1:331245/175‑1 (MQ=255)                                                                                   aCGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTg                                                                                                                                     >  1:334526/1‑251 (MQ=255)                                                                                      aTGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGt                                                                                                                                  <  1:23708/251‑1 (MQ=255)                                                                                                        acaaCGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTAAGCCGCGAGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGAGAAGGGGCGCGCTTTGCAAAACTTCTGCCGGAATATAAAGCCGCGCAGGAAATTCCTCGCTCGCGACTGTAATTAGAGACGATGGAAAAAGTGTTGGGTAAAACCCGCAAAggactg                                                                                                                >  2:211544/3‑246 (MQ=255)                                                                                                                 aGCAATACATTCGTGAAGCAGAAGCATATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGAt                                                                                                       <  1:19966/251‑1 (MQ=255)                                                                                                                   cAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGAGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGTAGTGGCGCGCTTTGCTAAACTTCTGCCGAAATATAAAGCCGCACCGGAAATTCATCGCGAGCGTCTGTATATCGGGAAGGTGGAAAAAGTGTTGGGTAAAACCCGCACAGTGCTGGTTatcaata                                                                                                      >  2:7847/1‑244 (MQ=255)                                                                                                                      aaCATTCATGAGGCACAAGCGGATACCAATGTAGTTTAGGCGCGTGCGAACGGTCATGCGCACCGTATCCTCGAAGAGGCGCGTGCGTACAAGGACCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGATCTTTGCTAAACTTATGCCGTAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTAGATCGAGACGATGGAAAAAGTGTTGGGTAATACCCGCAAAGTTCTGGGTAACGATAAAgg                                                                                                  <  2:278197/250‑1 (MQ=255)                                                                                                                             gTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCg                                                                                                                                                         >  1:332783/1‑189 (MQ=255)                                                                                                                             gTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGGGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCTGAAATTACTCGCGAGCGTCGGTATATCg                                                                                                                                                         <  2:332783/189‑1 (MQ=255)                                                                                                                              tGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTcgc                                                                                                                                                                         >  1:374052/1‑172 (MQ=255)                                                                                                                              ggAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGGATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGTAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCGGGAAATTACTcgc                                                                                                                                                                         <  2:374052/171‑1 (MQ=255)                                                                                                                                            atgccaacgcagtcctgtcgcgTGCGAACTGTCGGGCGCATCGTATCCTCGAAGAGGCGCGGGCGTACAAGGCCCAGCCCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGTATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGTGTAACGCTCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGcc                                                                             <  2:280598/232‑1 (MQ=255)                                                                                                                                              aCCAACGAAGTTCAGCCGCGTGCGAACTGTCAGGCGCAACGTATCCGCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGtt                                                                          <  1:102287/251‑1 (MQ=255)                                                                                                                                                         tCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCga                                                                                                                                                        <  1:308562/162‑1 (MQ=255)                                                                                                                                                           aGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTg                                                             <  1:272247/251‑1 (MQ=255)                                                                                                                                                                     gAACGGTCAGGCGCAACGTATCCTCGCAGAGGCGCGTGCGTAAAAGGCCCAGACCATACTGGAAACTCAGGGTGGAGTGTCGCGCTTTGCTAAACTTATTCCGGAATATAAAGCCGAGCCGGCAATTACTCGCGAGCGTCTGTATATCGAGACGAAGGAAAAAGTGTTGGTTAACACCCGACATGTGATGGTTAACGATAAAGGGGGCAACCTGAGAGTTCATCCGGTAGACCAGATGATGAACggtggt                                                    >  2:145929/1‑250 (MQ=255)                                                                                                                                                                       aCGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACACAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAGCCTGATGGTTCTACCGTTAGACCAGATGCTGAAAGGTGGTAAc                                                 <  1:296510/251‑1 (MQ=255)                                                                                                                                                                               gcgcAACGTATCCTCGCAGAGGCGCGTCAGTACAAGACCAAGCCCACCCTGGAAACTCATGGGGAAGTGGCGCGCTTTACTAAACTTCTGCCGGAAAATAAAGCCGCGCCGGGAATTACTCGCGAGCATCTGTATTTCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAGGTGCTGGGTAACAACAAAGGTGGGAACCTGATGGTTCTGCCGTTAGCCCAGAAGCTGaccgtgggtagagaccctgg                                         >  2:233928/1‑232 (MQ=255)                                                                                                                                                                                  cAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTgcggc                                      <  1:128519/251‑1 (MQ=255)                                                                                                                                                                                             cGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGGAACTCAGGGTGAAGTGACGCGCTTTGCTACACTTCTGCCGGAATATCAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATAGAGACGATGGAAAAAGTGTTGGGTAACACCCGCCAAGTGCTGGTTAACGAGAAAGGTGGCAACCTGATGGTTCTGCAGTTAGACCAGATGCTGACAGGTGGTAACGCCCCTGCGGCGAAGAGCgct                            >  2:61881/1‑248 (MQ=255)                                                                                                                                                                                                  agGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGCTGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACggg                      >  1:24314/1‑250 (MQ=255)                                                                                                                                                                                                       cgtgcgtACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAg                 <  1:34534/251‑1 (MQ=255)                                                                                                                                                                                                        gtgcgtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGATTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGAGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAATACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTACTGCCGTTAGACCAGATGATGAAAGGTGGTAACGCCCCTGCGGCGAAAAGCGATAACGGTGCCAGc                >  2:128897/1‑251 (MQ=255)                                                                                                                                                                                                          gcgtACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTGAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGGGTTGGGTAACACCCGCAAAGTGCTGGGTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGGCCAGATGCTGAAAGGAGGTAACGCCCCTGCGGCGAAAGGCGATAACGGTGCCAGCa               >  1:32340/1‑250 (MQ=255)                                                                                                                                                                                                                aaGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATctgct        <  1:211544/251‑1 (MQ=255)                                                                                                                                                                                                                    cccAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGt    <  1:329368/251‑1 (MQ=255)                                                                                                                                                                                                                      aataccctCCTGGAAACTAAGGGTAAAGTGGAGCGCGTTGCTAAACTTATGCCGGAATATAAAGCCGCGCCGGAAATTAATCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGCGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGAGGCAACCTGATGGATCTGCCGTTAGACCAGATGCTGAAAGGTGGTACCGCACCTGCGGAGCAGAGCGATAACAGAGACAGCAATCAGCTGCGTCt  >  2:341780/8‑251 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            CGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCT  >  NC_000913/4402710‑4403173

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