Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	13,029	C→T	D289D (GAC→GAT)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	13,029	0	C	T	98.5%	208.9 / ‑6.5	66	D289D (GAC→GAT)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (39/26);  minor base G (0/1);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

GACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCT  >  NC_000913/12784‑13272                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    gacgaAGTTGACGGCGAAAAAACCTTCGACGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGAc                                                                                                                                                                                                                                                >  2:185390/1‑251 (MQ=255)   cgaATTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTTCCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTACCCTGCCGTACATCGCTGCAGATGCGACCg                                                                                                                                                                                                                                              <  1:482719/251‑1 (MQ=255)        ttGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGATGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGACTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCg                                                                                                                                                                                                                                         >  1:375631/1‑251 (MQ=255)             ggCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGACGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCGCTGCAGATGCGACCGGTCCGAAaca                                                                                                                                                                                                                                    <  2:349468/251‑1 (MQ=255)              gCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTACCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAacac                                                                                                                                                                                                                                   <  2:162937/251‑1 (MQ=255)                           cGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAACGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGCTGCGACCGGTCCGAAACACATGAACATCaa                                                                                                                                                                                                                        >  1:359086/1‑249 (MQ=255)                               gTTCTGTCAACCAACGGTGAGACCCACTTGGGTGGTGAAG‑ACTGCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATGGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACGGTCTTCCGCGCAGCAGACCGACGTTAACCTGCCGTACATCACTGGAGATGCGACCGGTCCGAAACACAGGAACATCAAAGttac                                                                                                                                                                                                                  <  2:62050/251‑4 (MQ=255)                               gTTCTGGCAACCAACGGTGATACCCCCCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGGTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGAc                                                                                                                                                                                                                  <  2:53247/251‑1 (MQ=255)                                    ggCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCg                                                                                                                                                                                                               <  1:6800/249‑1 (MQ=255)                                           aacGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACt                                                                                                                                                                                                      >  1:79777/1‑251 (MQ=255)                                                 gatccccACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCGGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACAGGAACATCAAAGTGACTCGTGCGAAACTGGAAAg                                                                                                                                                                                                <  2:359086/247‑1 (MQ=255)                                                 aaaaCCCACCTGGGTGGTGAAG‑ACTTCGAATGCCGTCTTAGCAACTTCCTGGGTGAAAAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAATTGACTCGTGCGAAACTGGAAAg                                                                                                                                                                                                <  2:80329/248‑1 (MQ=255)                                                     cccACCTGGGGGGAGAAGCACT‑CGCCACCAGTCTGATCACCTACCTGGTTGAAGAATTCAGGAACGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGACAAAGCGAAAATCGAACTGTCTTCCGCTAAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACCTGAACATCAAAGTGCCTCGTGCGAAACTGGGAAGCCTg                                                                                                                                                                                            >  2:621188/1‑251 (MQ=255)                                                          cTGGGTGGTGAAG‑ACTTCGCCAGCCGTCTGATCACCTACCTGGTTGAAGAATTCAAGAAAGATAAGGGTCTTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTg                                                                                                                                                                                        >  2:62082/1‑250 (MQ=255)                                                             ggtggTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGa                                                                                                                                                                                    <  1:41494/248‑1 (MQ=255)                                                                   gAAG‑ACTTCGACAGCCGTCTGATTAACTACCTGGTTGAAGAATTCAATAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTTCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTgg                                                                                                                                                                               <  1:393122/250‑1 (MQ=255)                                                                     aG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTaa                                                                                                                                                                            <  1:394471/251‑1 (MQ=255)                                                                           tCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAAcc                                                                                                                                                                         <  1:352166/249‑1 (MQ=255)                                                                              aCAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGCAGATCTGGTAAACCGTTcc                                                                                                                                                                    >  1:605826/1‑251 (MQ=255)                                                                                aGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCAt                                                                                                                                                                  >  2:5916/1‑251 (MQ=255)                                                                                    gTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTg                                                                                                                                                                >  1:65167/1‑249 (MQ=255)                                                                                      cTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTG‑CAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGcc                                                                                                                                                            >  2:418949/1‑250 (MQ=255)                                                                                        gATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTG‑CAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGcc                                                                                                                                                            <  1:418949/248‑1 (MQ=255)                                                                                           cAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTg                                                                                                                                                        >  1:171869/1‑250 (MQ=255)                                                                                                  cTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGAAACGATCCGCTGGAAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGCCCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGc                                                                                                                                                >  2:412961/1‑251 (MQ=255)                                                                                                    tgTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGAGCCGCTGGCAATGCAGGGCCTGAAAGAAGCGTCCGAAAAATCGAAAATCGAACTGTGTTCCGCTCAGCAGTCCGATGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACAGCAAAGTGACGCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCa                                                                                                                                               <  1:581646/249‑1 (MQ=255)                                                                                                           gaaTTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAgg                                                                                                                                        >  1:35673/1‑250 (MQ=255)                                                                                                            aatgcAAGAAAGATCAGGGTATTGCCCTGCGCAACGATCCGCTGGCAATGCAGCTCCTGAAAGAAGCGGCAGAAAAAGCGAAACTCTAACTGTCTTCCGCTCATCAGACCGACGTTAACCTGCCGTACATCACTGCAGAGGCGACCGGTCCGAAACACATGAACAGCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGTGCCGCTGAAAGTTTCCCTGCAggcc                                                                                                                                      <  2:222136/247‑3 (MQ=255)                                                                                                              ttCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGc                                                                                                                                    >  1:366800/1‑251 (MQ=255)                                                                                                                  agaaagaTCCGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGc                                                                                                                                >  1:517080/1‑251 (MQ=255)                                                                                                                    aaagaTCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCt                                                                                                                              <  1:116228/251‑1 (MQ=255)                                                                                                                     aagaTCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGGAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTg                                                                                                                             >  2:216500/1‑251 (MQ=255)                                                                                                                       gaTCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGACTGTc                                                                                                                           >  2:182391/1‑251 (MQ=255)                                                                                                                                   gACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGatat                                                                                                               >  1:431767/1‑251 (MQ=255)                                                                                                                                     ccTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATcg                                                                                                             >  2:379071/1‑251 (MQ=255)                                                                                                                                              cGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGtt                                                                                                     >  1:594652/1‑250 (MQ=255)                                                                                                                                               gATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTACGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGAATCTGATATCGACGACGTTAt                                                                                                   >  1:578417/1‑251 (MQ=255)                                                                                                                                               gATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGAAGAAAAAGCGAAAAGCGAACTGTCTTCCGCTCAGCAGGCCGACGTTAACCTGCCGTACATCACGGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTTAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTAt                                                                                                   <  2:300814/251‑1 (MQ=255)                                                                                                                                                   cgctggacaTGCAGCGCATGAAAGAAGCGGCAGACAAAGCGAATATCGAACTGTATTCCGCTAAGCAGACAGAAGTTAACCTGCCGTACATCACTGCAGATGAGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGACAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGCTATCGAAGACGTTGTCCTc                                                                                               >  2:408606/9‑251 (MQ=255)                                                                                                                                                       ggCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTCCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGt                                                                                             <  1:443930/249‑1 (MQ=255)                                                                                                                                                        gCAATGCCGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACGGGAAAGCCGGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTtgg                                                                                          <  2:605826/251‑1 (MQ=255)                                                                                                                                                             gCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTc                                                                                     <  1:185390/251‑1 (MQ=255)                                                                                                                                                             gCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTAGCCTCGTTggtggt                                                                                      <  1:137802/250‑1 (MQ=255)                                                                                                                                                                   ccTGAAAGAAGCTGCAGAAAAAGCGAAAATCGAAGTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTc                                                                               <  2:265496/251‑1 (MQ=255)                                                                                                                                                                                gCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGt                                                                  >  2:15055/1‑251 (MQ=255)                                                                                                                                                                                 cAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTACCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGGAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTAGCTGATATCGACGAAGTTATCCTCGTTGGTGGTCAGACTCGTAGGCCAATGGtt                                                                 >  2:145847/1‑251 (MQ=255)                                                                                                                                                                                  tgAAAAAGCGAAAGTCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCAATGCAGATGCGCCCGGTCCGAAACACATGAACATCATAGTGACTCGTGCGAACCTGGAACGCCTGATTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTGGAACTGCAGGACGTTGGCCTGTCCGTATCTGATTTAGACAACGTTATACTCGTTGGTGGTCaagcacggctgccaaatgtac                                                                >  2:589289/2‑231 (MQ=255)                                                                                                                                                                                      aaaGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCaga                                                             >  1:558666/1‑250 (MQ=255)                                                                                                                                                                                         gCGAAAATCGAACTGTCTGCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAgaagaa                                                         <  2:366800/251‑1 (MQ=255)                                                                                                                                                                                             acaTCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACAGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGCAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCATGTCCGTATCTGATATCGACGACGTTATCCTAGTTGGGGGTCCGGCTCGTATGCCAAAGGTTCAGAAGaaattt                                                     >  2:203574/3‑248 (MQ=255)                                                                                                                                                                                                tCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGACAGTTgat                                                  >  1:554252/1‑249 (MQ=255)                                                                                                                                                                                                 cGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCt                                                  <  1:286925/250‑1 (MQ=255)                                                                                                                                                                                                  gAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTg                                                 >  1:162301/1‑250 (MQ=255)                                                                                                                                                                                                   aaCTGTCTTCCGCTCAGCAGACAGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGCAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGACAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTg                                                 >  2:559805/1‑249 (MQ=255)                                                                                                                                                                                                          ttCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTg                                        >  1:420784/1‑251 (MQ=255)                                                                                                                                                                                                             cGCTCAGCAGACCGCCGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGt                                      >  1:469150/1‑250 (MQ=255)                                                                                                                                                                                                             cGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGt                                      >  1:431298/1‑250 (MQ=255)                                                                                                                                                                                                                tCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAAACGTTCCATTg                                                                                                                                                                >  1:608348/1‑125 (MQ=255)                                                                                                                                                                                                                tCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACAGGAACATCAAAGTTACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCGGGTAAACCGTTCCATTg                                                                                                                                                                <  2:608348/125‑1 (MQ=255)                                                                                                                                                                                                                  agcagACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGGTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAag                                  >  2:396952/1‑249 (MQ=255)                                                                                                                                                                                                                      gACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAACGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCcgc                            >  2:220883/1‑251 (MQ=255)                                                                                                                                                                                                                        ccGACGTTAACATGCAGTACATCACTGCAGATGCGACCGGTCCGCAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCcgcg                           >  2:514502/1‑250 (MQ=255)                                                                                                                                                                                                                          gACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTa                         >  2:581724/1‑250 (MQ=255)                                                                                                                                                                                                                                     gccgtACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAAcc             <  1:274912/247‑1 (MQ=255)                                                                                                                                                                                                                                       cgtACATCCCTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGGTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCg           <  2:536695/249‑1 (MQ=255)                                                                                                                                                                                                                                                cTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCt  >  1:128911/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    GACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCT  >  NC_000913/12784‑13272

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