Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,213,557	C→T	I108I (ATC→ATT)	yjaA →	stress response protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,213,557	0	C	T	97.1%	227.5 / NA	69	I108I (ATC→ATT)	yjaA	stress response protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (39/28);  minor base G (0/2);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATT  >  NC_000913/4213310‑4213796                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                gagaCGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCAt                                                                                                                                                                                                                                              <  2:129855/251‑1 (MQ=255)      gCTGCCTTCAGTACCCAGCGCCTGTTAATCGCCCATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGCACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATgcgc                                                                                                                                                                                                                                          >  2:507236/1‑250 (MQ=255)        tGCCTTAAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTCGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATAGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAAt                                                                                                                                                                                                                                       >  1:554190/1‑251 (MQ=255)           cTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATc                                                                                                                                                                                                                                      >  1:86895/1‑249 (MQ=255)                    ccAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCTCGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAAGGAGGGCGTTTTGTCACAGGTTGAGTAACGCATTCTCCATGAAGTGGTCGCTTGGGCAACGTTAATGAAAGATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCa                                                                                                                                                                                                                           <  2:695479/251‑1 (MQ=255)                     catcgcCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGTAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCTGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTGCCACATTCATGCGCAATCAGCGGTACTCAg                                                                                                                                                                                                                          <  2:403109/248‑1 (MQ=255)                       gcgcCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTg                                                                                                                                                                                                                        >  1:56322/1‑251 (MQ=255)                        cgcCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGa                                                                                                                                                                                                                       <  1:152173/251‑1 (MQ=255)                           cttttAATAGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTgct                                                                                                                                                                                                                      >  2:329832/4‑247 (MQ=255)                            tGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCCGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGt                                                                                                                                                                                                                   >  2:32134/1‑251 (MQ=255)                             gTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAg                                                                                                                                                                                                                    <  1:492178/249‑1 (MQ=255)                                aaTCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCgg                                                                                                                                                                                                               >  2:183454/1‑251 (MQ=255)                                     ccAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATg                                                                                                                                                                                                          >  2:150742/1‑251 (MQ=255)                                               tgttgCAGAAAAAGTTCTGCCAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCtt                                                                                                                                                                                                 >  1:617396/1‑250 (MQ=255)                                                  tgCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAg                                                                                                                                                                                              <  1:28645/250‑1 (MQ=255)                                                  ggCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACTCACGTTCAACCTGTCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGAGTACTGGAAATGAATGAGGGCGGTTTGTCACAGTTTGAGGAACGCAGTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCCGTTCCACATTCATGCGCAATCAGCGGAACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGc                                                                                                                                                                                             >  2:132361/2‑251 (MQ=255)                                                        aaaaGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAaat                                                                                                                                                                                       >  1:26424/1‑251 (MQ=255)                                                           aGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAaataat                                                                                                                                                                                    >  1:211319/1‑251 (MQ=255)                                                                tGCAAGATCTTGTTCTGAAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCgtgt                                                                                                                                                                               >  2:242404/1‑251 (MQ=255)                                                                   aagattgtGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATAGTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTGTGTCACAGGTTGAGGAACGCATTCGCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCAc                                                                                                                                                                            <  2:714/244‑1 (MQ=255)                                                                   aaGATCTTGTTCTGCAACTCCACCCACGTTCAACCGGGCATTCTTTTTTGCGAGCAAATCGTATGGATATTGTTGTGAGCGCACTGGAAAGGAATGAGGGCTGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGTTCGCTGGGGCAACGTTAATGAAATATCTCCAGTTCCGCATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCAc                                                                                                                                                                            <  1:674185/251‑1 (MQ=255)                                                                        cTTGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATtg                                                                                                                                                                       >  1:578451/1‑251 (MQ=255)                                                                          tGTTCGGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATtg                                                                                                                                                                       <  2:26424/249‑1 (MQ=255)                                                                               tGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGt                                                                                                                                                                 >  1:339469/1‑250 (MQ=255)                                                                                gCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGt                                                                                                                                                                 <  1:267552/249‑1 (MQ=255)                                                                                      ccacccacGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATGCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAAcc                                                                                                                                                         <  1:665375/251‑1 (MQ=255)                                                                                        acccacGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGAATGGAGATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCTCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGGACTCAGCGATAGCGCGGTCATGGCAATGCTTAAGCGGAGATAATCGTGTCACCATTGGTGGGTACTAAACCTg                                                                                                                                                       >  1:561245/1‑251 (MQ=255)                                                                                                   aCCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGccc                                                                                                                                            <  1:662589/251‑1 (MQ=255)                                                                                                     cTGGCCTTCTTTTTTGCCAGCAAAACGTATGGATATTATTGTGAGCGCACTGGCAATGCCTGAGGGCGGTTTGTCCCAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCCTGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGCAATAATCGTGTCACCATTGGTGGGTACTAACCCTGAAGTTCAGCacac                                                                                                                                          >  2:460529/1‑251 (MQ=255)                                                                                                             cTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCa                                                                                                                                                                          <  2:314980/211‑1 (MQ=255)                                                                                                             cTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCa                                                                                                                                                                          >  1:314980/1‑211 (MQ=255)                                                                                                             cTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATg                                                                                                                                   >  1:21847/1‑250 (MQ=255)                                                                                                             cTTTTTTGCCAGAAAAACGTATGTATCTTTTTGTGAGCGCACTGTACATGAATGAGGGCGGATAGACACAGGTTGGGGAACGCATTCTCATTGCAGTGGTCGCTGTGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAAACAGCGGTACTCCGTGATAGTGAGGTCTTGCCAATGCTTAAGCAGCAATAATCGTGGCACCATTGGTGGGTAATAACCCTGAAGTTCAGCCAGCCGGGATgt                                                                                                                                  >  2:119838/1‑250 (MQ=255)                                                                                                                ttttGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTACGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACAGGGATgaga                                                                                                                                >  2:581612/1‑250 (MQ=255)                                                                                                                         cAAAACGTATGGATATTGTTGTGAGCGCACTGGAAAGGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGcc                                                                                                                      <  2:617396/251‑1 (MQ=255)                                                                                                                         cAAAACGTATGGAAATTGTTGTGAGCGCACTGGAAATGCATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGcc                                                                                                                      >  2:679935/1‑251 (MQ=255)                                                                                                                           aaaCGTATGGATAGTGTTGTGAGCGCACTGGAAATGAATGATGGCGGTTTGTCACAGGTTGAGTAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCTGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGATAAAAAATCGCCTa                                                                                                                    <  2:632357/251‑1 (MQ=255)                                                                                                                               gTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGGCACAGGTTGAGGAACGCATTCTCCATTAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTg                                                                                                                                                                                                                  <  1:664561/153‑1 (MQ=255)                                                                                                                               gTATGGATATTGTTGTGAGAGAACTGGAAATGAATGACGGCGGTTTGTCAAAGGTTGAGGAACGCATTATCCATGAAGTGGTCGCTGGGGCAAAGTTAATGAAATATCGCAAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTg                                                                                                                                                                                                                  >  2:664561/1‑153 (MQ=255)                                                                                                                                 aTGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGACACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATg                                                                                                                                                                                                          >  1:298537/1‑159 (MQ=255)                                                                                                                                 aTGGATATTGGTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGccc                                                                                                               <  1:544869/250‑1 (MQ=255)                                                                                                                                 aTGGAGATTGTTGTGAGCGCACTGGAAAGGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCAGTCGCCATGAAGTGGTCGCTGGGGCAACGGTAAGGAAATATCGCCAGGTCCACATTCATGCGCAATCAGCGGTGCTCAGTGATAGTGCGGTCATg                                                                                                                                                                                                          <  2:298537/159‑1 (MQ=255)                                                                                                                                    gATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCcac                                                                                                           >  1:212783/1‑251 (MQ=255)                                                                                                                                                       tGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAAc                                                                                        >  2:539286/1‑251 (MQ=255)                                                                                                                                                         gAAATGAATGAGGGCTGTTTGTCCCAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGATGGGGCAACGTTAATGACATATCGCCAGTTCCACATTCATGCGCAATCAGCGATACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAACCCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACgg                                                                                      <  2:162178/251‑1 (MQ=255)                                                                                                                                                             ggaatgaGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATAc                                                                                  <  2:561245/250‑1 (MQ=255)                                                                                                                                                               aatgaGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACgg                                                                                <  1:150742/251‑1 (MQ=255)                                                                                                                                                                aggaGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGt                                                                               <  2:578451/249‑1 (MQ=255)                                                                                                                                                                   aGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTc                                                                            <  1:32134/251‑1 (MQ=255)                                                                                                                                                                    gggCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTcc                                                                           <  1:473016/251‑1 (MQ=255)                                                                                                                                                                       cGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTcccc                                                                         >  2:83919/1‑250 (MQ=255)                                                                                                                                                                        ggTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCa                                                                        >  1:599579/1‑250 (MQ=255)                                                                                                                                                                                    gTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTc                                                           >  1:389606/1‑251 (MQ=255)                                                                                                                                                                                    gTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAAGGAAATAGCGCCAGTTCCACATGCATGCGCAAGCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTGCCCCAAATCGTCCACCTc                                                           <  2:86895/251‑1 (MQ=255)                                                                                                                                                                                         ggAACACATTCGCCATGACGTGGGCGTTGGGGCACGGTTAATGAAACCTTGCAAGTTCAACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGCAGTTCAGCCCACCGGGATGAGACAAAATCGCCTACGCCCACACATACGCCAGATTCAGCAACGAATACGGCTTCCCACAATCGTCCACCTccgtgc                                                      >  2:478369/1‑246 (MQ=255)                                                                                                                                                                                                    tcCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTaccc                                           >  1:397276/1‑248 (MQ=255)                                                                                                                                                                                                     ccATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCt                                          <  1:66678/251‑1 (MQ=255)                                                                                                                                                                                                     ccATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAAGGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAAcc                                           <  1:460529/250‑1 (MQ=255)                                                                                                                                                                                                      cATGAAGTGGGCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACAGTCATGCGCAATCAGCGGTACTCAGTGATAGGGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAAGCGCCTACGCCCCCACATACGCCAGATGCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCtt                                         <  2:339469/251‑1 (MQ=255)                                                                                                                                                                                                         gAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTaa                                       >  1:643122/1‑250 (MQ=255)                                                                                                                                                                                                          aaGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTaaa                                      >  1:489230/1‑250 (MQ=255)                                                                                                                                                                                                                  gCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACctt                             >  1:133735/1‑251 (MQ=255)                                                                                                                                                                                                                     ggggCAACGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACcttctt                          >  1:240023/1‑251 (MQ=255)                                                                                                                                                                                                                          aaCGTTAATGAAATATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTtat                        >  1:391767/1‑248 (MQ=255)                                                                                                                                                                                                                              ttAATGAAATATGGCCAGTTCCACATTCATGCGCAATCAGCGGTACACAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGCAATAATCATGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCATACGCACCCACATACGCCAGATTCAACAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGCAccc                 >  2:345422/1‑251 (MQ=255)                                                                                                                                                                                                                                     aaTATCGCCAGTTCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGcc          <  1:132728/251‑1 (MQ=255)                                                                                                                                                                                                                                             cAGTTCCACATTCATGCGAAATCAGCGGTACTCAGTGATAGGGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGAACTAAACCTGAAGTTCAGCCAACCGGGATTAGAAAAGATCGCCTACGCCCCCACATACGCCGGATTCAGCAACGGAGACGGTTTACCCAAATAGTCCACCTCAGAGCGTCCCGGAACCTTAATACCCACATTCTTATAGCACCCAACCGCCTGCTCAtt  >  2:534142/1‑251 (MQ=255)                                                                                                                                                                                                                                                ttCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCgg                                                                                                                                       >  2:271661/1‑115 (MQ=255)                                                                                                                                                                                                                                                ttCCACATTCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCgg                                                                                                                                       <  1:271661/115‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                GAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATT  >  NC_000913/4213310‑4213796

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