Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,580,388	G→A	P251P (CCG→CCA)	tktB →	transketolase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,580,388	0	G	A	95.2%	125.5 / NA	43	P251P (CCG→CCA)	tktB	transketolase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (20/20);  minor base C (2/0);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAAGCTGATTGGTTTTTACGATCACAACGGTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCGAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATGAGAAGTTTGCCGCCTATAAAAAGGCTCATCCGCAACTGGCAGAAGAGTTTACCCGACGGATGAGCGGTGGTT  >  NC_000913/2580159‑2580629                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   gCAAGCTGATTGGTTTTTACGATCACAACGGTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAagag                                                                                                                                                                                                                              <  1:167564/251‑1 (MQ=255)   aaGCTGATTGGTTTTTCCGATCACAACGGTATTCCCATCGACGGTGAAACAGAAGGCTGGTTTCCCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAACCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAATCGTGAAAGATACGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAggc                                                                                                                                                                                                                            <  2:185672/251‑1 (MQ=255)         aTTGGTTTTTACGATCACAACGGTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGCACCTTGAAGCGCAAAGAGTGAAAGATAAGCCGTCGCCGCTTATCTGCCGGACGGTGATTGGCTTTGTATCGCC‑AACAAACGCAGGTGAGGGAGCggcgcacgg                                                                                                                                                                                                                      >  2:231216/1‑251 (MQ=255)          ttGGTTTTTACGATCACAACGGTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGACGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GAAGGTAAGGAAGAggcgcccggc                                                                                                                                                                                                                     >  1:15514/1‑251 (MQ=255)             gTTTTTACGATCACAACGGTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGACGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAggcgcacggcgc                                                                                                                                                                                                                   >  1:450988/1‑250 (MQ=255)                 ttACGATCACAACGGTATTGCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGAACCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGcacca                                                                                                                                                                                                               <  1:431641/250‑1 (MQ=255)                  aacgaacacaacggtaaggtcaacGACGGTGAAACAGAAGGCTGGTTTAACGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGGGAAGGAAGCGATCCGTGAAGCGCAAAGCGTTAAAGATAAGCCGACGCTGATTATCTGCCGTACGGTGATTGGCTGTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACt                                                                                                                                                                                                             <  2:32259/228‑1 (MQ=255)                             cGGTAGTTCCATCGTCGGTGAAACAGAAGGCTGGTTTACCGACGATACGTCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTGGAAGCGCAAAGCTTGAAAGATAAGCCGTCTCTGATTAGCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGTAAGAGGCGCACGGCGCACCACTGGGGtaagaag                                                                                                                                                                                                  <  1:436367/251‑1 (MQ=255)                               gTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGgaagaagaa                                                                                                                                                                                                <  1:426754/251‑1 (MQ=255)                                         cGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTgg                                                                                                                                                                                      >  1:341097/1‑251 (MQ=255)                                          gACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGGAGAAGAAGGGGCGCTGGc                                                                                                                                                                                     >  1:207469/1‑251 (MQ=255)                                                gAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGc                                                                                                                                                                                <  2:15514/250‑1 (MQ=255)                                                         ggCTGGTTTACCGACGACACGGCAAAACGTTTTGTAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGACCCGCAGGAGGTGAAGGAAGCTATCCTTGAAGCGCAAAGCGGGGAGGATAAGGCGACGGTGCTTAACTGCCGTACGGTGAATGGCTTTGGTTAGCCCAATAAA‑GCCGGTAAGGAAGCGGCGAACAGGGCCCCACTCGGGGAAGCAGAAGGGGCCCTGGAACGGCaaaaaacgcg                                                                                                                                                                      >  2:156997/1‑246 (MQ=255)                                                                  aCCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCCCCACTGGGGGAAGAAGAAGTGGCGCTGGCAAGGCCAAAACTGGGCTGGCAcc                                                                                                                                                              >  1:448140/1‑250 (MQ=255)                                                                  aCCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCAcc                                                                                                                                                              <  1:427867/250‑1 (MQ=255)                                                                       cgaTACGGCAAAACGTTGTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGc                                                                                                                                                        <  2:75482/251‑1 (MQ=255)                                                                             ggCAAAAAGTTTTGAAGCCTATCACTGGACTGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAAGAAGCGATCCTTGAAGCGCCAAGCGTGAAAGATAAGCCGTAGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGCGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCCAAAACTGGGATGGCACCCTCCCCCAtttt                                                                                                                                                  >  2:264662/1‑250 (MQ=255)                                                                              gCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATACGCAGGCGGTGAAGGAAGCGCTCCTTGAAGCGCAAAGCGTGAAAGATCAGCCGTCGCTGCTTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGCAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAAGAGGGCTGGCACCAACCGCCATTTgg                                                                                                                                                 >  2:499637/1‑250 (MQ=255)                                                                                          gAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAag                                                                                                                                      <  1:39186/250‑1 (MQ=255)                                                                                          gAAGCCTATCAATGGCATGTGATCCCTGAAATCGCCGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGGAAAGCGTGAAAGAAAAGCAGGCGATAATAAACTGCCGTACGGAGATTGGCTTTGGTCCGCCCAATAAA‑TCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGGAGAAGGAGTGGAGACGGAACGGCAAAAACTGGGCTGGAACCAACCGACATTTGAGAtacattaga                                                                                                                                      >  2:361379/1‑242 (MQ=255)                                                                                                tATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCCGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGCAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCCTCCGCCATTTGAGATACCTAAAAAGAtttt                                                                                                                               >  2:278274/1‑250 (MQ=255)                                                                                                     cTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGGAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAAATGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACg                                                                                                                          >  2:415316/1‑251 (MQ=255)                                                                                                                   aTGAAATCGCCGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGGGGAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGGTTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGGAAGGCGCACCACTGGGGGGAAAAGAAGTGGCTCTGGCACGGCCAAAACTGGGCAGGCCCCATCCGCCAGTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGt                                                                                                            >  2:119931/1‑251 (MQ=255)                                                                                                                    tGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGCAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGACGTAAGGTGATTGGCATTGGATCGCCAAATAAA‑GCAGGGAAGGAAGAGGCGCACGGCGCCACACTGGGGGAAGAAGAAGCGGCGCTGGGACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTg                                                                                                           >  1:48662/1‑251 (MQ=255)                                                                                                                     gAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGGACGGGGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGGAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGACACAGTGGCGCTGGCACGGCCAAAACTGGGCTGGCCCCATCCGCCATTTGAGATCCCTACAGAGACTTATAACGCCTggggagcccgcg                                                                                                           >  2:447403/1‑241 (MQ=255)                                                                                                                       aaTCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAATGAATAGGCGCACGGCGCTCCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCTCGTGaaa                                                                                                        <  1:373831/251‑1 (MQ=255)                                                                                                                                   cGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCGGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCCCTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCAGCATCCGCCAGTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAg                                                                                            <  1:531785/251‑1 (MQ=255)                                                                                                                                    gATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAgg                                                                                           >  1:269020/1‑250 (MQ=255)                                                                                                                                          cAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAGGGCGAAAAAGCgcagca                                                                                     <  1:278274/251‑1 (MQ=255)                                                                                                                                                            aTCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATGAGAAGtt                                                                   >  1:502937/1‑251 (MQ=255)                                                                                                                                                                      cgcAAAGCGTGAAAGATAAGCCGGCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGACAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGGGGCGTTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATGAGAAGTTTGCCGCCtat                                                         <  2:448140/251‑1 (MQ=255)                                                                                                                                                                                 aaaGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAAt                                                                           >  2:150650/1‑222 (MQ=255)                                                                                                                                                                                 aaaGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAAt                                                                           <  1:150650/222‑1 (MQ=255)                                                                                                                                                                                      tAAGCCGTGGCTCATGATCAGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGATGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCGAAAGAGATTTATCACGACTAGGATGCCCGTGAAAAAGACGAAAAAGCGCAGCAGAGCTGGAATGAGAAGTTTGGCGCCTATAAAAAGGCTCATCCGc                                         >  2:373629/1‑251 (MQ=255)                                                                                                                                                                                                   aTTATCTGCCGTACGGTGATAGGCTTTGGTTCGCAAAATAAA‑GAAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATACCTAAAGAGATTTATCACGCCTGGGATGCCCGAGAGAAAGGCGCAAACGCGCAGCAGTGCTGGAAGGAGAAGTTTGCAGCCAAGAAAAAGGCTAATCCGAAACTGGCagc                               >  2:138257/1‑247 (MQ=255)                                                                                                                                                                                                         tGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCAttt                                                                                                                                                   <  2:388234/126‑1 (MQ=255)                                                                                                                                                                                                         tGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCAttt                                                                                                                                                   >  1:388234/1‑126 (MQ=255)                                                                                                                                                                                                          gCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCTCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAATGGCGAAAAAGCGCAGCAGAGCTGGAATGAGAAGTTTGCCGCCTATAAAAAGGCTCATCCGCAACTGGCAGAAGAGttt                        <  1:230250/248‑1 (MQ=255)                                                                                                                                                                                                           cctgacggtagttGGCTTTGGTTCGCCAAATAAA‑GCAGGTGAGTAAGAGGCGCACGGCGCCCCACGGGGGGAAGAAGAAGGTGCGCGGGCACGGCAAAAACTGGTCGGGCACCATCCGCTATTTGAGATCCCTAAAGAGATTTAGCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATGAGAAGTTTGCCGCCTATAAAAAGGCTCATCCGCAACTGGCATAAGAGGTTAccc                    <  2:471976/240‑1 (MQ=255)                                                                                                                                                                                                               aCGGTGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATGAGAAGTTTGCCGCCTATAAAAAGGCTCATCCGCAACTGGCAGAAGAGTTTACCCGACg                <  1:499637/251‑1 (MQ=255)                                                                                                                                                                                                                   tGATTGGCTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATGAGAAGTTTGCCGCCTATAAAAAGGCTCATCCGCAACTGGCAGAAGAGTTTACCCGACGGATg            >  1:405041/1‑251 (MQ=255)                                                                                                                                                                                                                          cTTTGGTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGAAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGACAAAGCGCAGCAGAGCTGGAATGAGAAGTTTGCCGCCTATAAAAAGGCTCATCCGCACATGGCCGCAGAGTTTACACGCAGGATGAGCggtg     >  2:183769/1‑251 (MQ=255)                                                                                                                                                                                                                              ggTTCGCCAAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATGAGAAGTTTGCCGCCTATAAAAAGGCTCATCCGCAACTGGCAGAAGAGTTTACCCGACGGATGAGCGGTGGtt  <  1:264662/250‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   GCAAGCTGATTGGTTTTTACGATCACAACGGTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCGAATAAA‑GCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATGAGAAGTTTGCCGCCTATAAAAAGGCTCATCCGCAACTGGCAGAAGAGTTTACCCGACGGATGAGCGGTGGTT  >  NC_000913/2580159‑2580629

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