Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,557,375	A→T	8.0%	N128I (AAT→ATT)	iraD →	RpoS stabilzer after DNA damage, anti‑RssB factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,557,375	0	A	T	8.0%	204.8 / 8.1	63	N128I (AAT→ATT)	iraD	RpoS stabilzer after DNA damage, anti‑RssB factor
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base T (2/3);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01

TTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTAAAAGCTCGCTCATCGAAGAGACAAAAAAGGTGCTGAAACTTAACTCACTGTGGCGCGTCAGGGTGCGATTA  >  NC_000913/4557089‑4557660                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                              tttAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATg                                                                                                                                                                                                                                                                                             >  2:118456/1‑289 (MQ=255)   tAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGATTGtg                                                                                                                                                                                                                                                                                           <  2:463382/289‑2 (MQ=255)       tctcCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCt                                                                                                                                                                                                                                                                                       <  1:149501/289‑1 (MQ=255)              aTTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGcaca                                                                                                                                                                                                                                                                                <  2:350691/289‑1 (MQ=255)                             cATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTcgcg                                                                                                                                                                                                                                                                  <  2:368548/288‑1 (MQ=255)                              aTTACTTGTTTCAGGCCTTCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACATCTCCAGCGCTGGTATTGCCGTCAGCTGCGTTCTTCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCTATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGCATGTCAGCTAAGCACACTCTCCCGCTCGCGCa                                                                                                                                                                                                                                                                >  1:350867/1‑289 (MQ=255)                              aTTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTTGAGAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTcgcg                                                                                                                                                                                                                                                                  <  2:378428/289‑1 (MQ=255)                              aTTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGATTGtgc                                                                                                                                                                                                                                                                                          <  2:248895/263‑3 (MQ=255)                              aTTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGATTGtgc                                                                                                                                                                                                                                                                                          >  1:248895/1‑261 (MQ=255)                                 aCTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGtt                                                                                                                                                                                                                                                             <  1:118456/289‑1 (MQ=255)                                                     aGTGGCAAGTGCATAACCCGAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTcgcg                                                                                                                                                                                                                                                                  >  1:266380/1‑264 (MQ=255)                                                     aGTGGCAAGTGCATAACCCGAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTcgcg                                                                                                                                                                                                                                                                  <  2:266380/264‑1 (MQ=255)                                                                  tAACCCGAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCaa                                                                                                                                                                                                                                        <  2:241567/277‑1 (MQ=255)                                                                  tAACCCGAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCaa                                                                                                                                                                                                                                        >  1:241567/1‑277 (MQ=255)                                                                        gAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCCGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTAt                                                                                                                                                                                                                      >  1:230738/1‑289 (MQ=255)                                                                        gAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCCCAATCTCCCGCTCGCGCAGTTCCCGCACGGAGCGTTCGGACACGGGCTTCACGCGGGTAt                                                                                                                                                                                                                      >  2:410023/1‑289 (MQ=255)                                                                        gAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGtt                                                                                                                                                                                                                                              >  1:432318/1‑265 (MQ=255)                                                                        gAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGtt                                                                                                                                                                                                                                              <  2:432318/265‑1 (MQ=255)                                                                        gAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTAt                                                                                                                                                                                                                      >  1:152375/1‑289 (MQ=255)                                                                        gAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTAt                                                                                                                                                                                                                      >  2:488979/1‑289 (MQ=255)                                                                                 aacattcaggagCTGTTATTGCCGTCAGCTGCTTTCTGCACGTTTATTTCTTGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGt                                                                                                                                                                                                             <  2:284037/280‑1 (MQ=255)                                                                                 aCATCTCCAGAGGTGGTATTGCCGTCAGGTGCGTTCTGCACGGTTATTTCATGAGCCACGGATTGCGGCGTTGAAGGTTAATTTGAAAGAAGCGGACTGCAACACTCTGGCGATCAGCATT‑‑GAGATAATGCTGGACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGAGACGCTGGAGAAGGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCAAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGt                                                                                                                                                                                                             <  2:350867/289‑1 (MQ=255)                                                                                      tcCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCa                                                                                                                                                                                                        <  1:401439/289‑1 (MQ=255)                                                                                      tcCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCa                                                                                                                                                                                                        <  1:68345/289‑1 (MQ=255)                                                                                                gTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTCCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACTCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAACCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGaa                                                                                                                                                                                              >  1:224134/1‑289 (MQ=255)                                                                                                          cAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCa                                                                                                                                                                                     <  2:48542/288‑1 (MQ=255)                                                                                                           aGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGACTGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGGCAGAAGCTTTCGCACACGGGATTCATGAGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCAcc                                                                                                                                                                                   >  2:197547/1‑289 (MQ=255)                                                                                                                    ttgataggTTATGTCAGGAGCCCCGGATTGCGGGGTGGAAGGTTAAATTGAGAGAGGCGTGCTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGGTGAACCGCTAACGTTTGATGTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTcc                                                                                                                                                                          <  2:296105/282‑1 (MQ=255)                                                                                                                                    tGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCg                                                                                                                                                                                                                           >  1:43880/1‑224 (MQ=255)                                                                                                                                    tGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCg                                                                                                                                                                                                                           <  2:43880/224‑1 (MQ=255)                                                                                                                                      aGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGa                                                                                                                                                          >  1:20263/1‑287 (MQ=255)                                                                                                                                           cGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGtt                                                                                                                                                                                                                                                             >  1:86822/1‑183 (MQ=255)                                                                                                                                           cGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGtt                                                                                                                                                                                                                                                             <  2:86822/183‑1 (MQ=255)                                                                                                                                                               tAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGa                                                                                                                               <  2:20263/289‑1 (MQ=255)                                                                                                                                                                       aagaagCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGATTGt                                                                                                                                                                                                                                                                                            <  2:446153/124‑1 (MQ=255)                                                                                                                                                                       aagaagCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGATTGt                                                                                                                                                                                                                                                                                            >  1:446153/1‑124 (MQ=255)                                                                                                                                                                                            tGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCAt                                                                                                                                                                                                                  <  2:39950/177‑1 (MQ=255)                                                                                                                                                                                            tGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCAt                                                                                                                                                                                                                  >  1:39950/1‑177 (MQ=255)                                                                                                                                                                                                gATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAg                                                                                               <  1:488979/288‑1 (MQ=255)                                                                                                                                                                                                        tt‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAgg                                                                                      >  2:233436/1‑289 (MQ=255)                                                                                                                                                                                                        tt‑‑GAGACAATGCTGTACCAAGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAgg                                                                                      <  2:196389/289‑1 (MQ=255)                                                                                                                                                                                                                        tACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGTCAGGCAATCCCACAGCCGACGAGGCCAACCTGCTTTaa                                                                        >  1:321600/1‑289 (MQ=255)                                                                                                                                                                                                                        tACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAAcc                                                                                <  1:31218/281‑1 (MQ=255)                                                                                                                                                                                                                        tACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAAcc                                                                                >  2:31218/1‑281 (MQ=255)                                                                                                                                                                                                                        tACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTaa                                                                        >  2:12116/1‑289 (MQ=255)                                                                                                                                                                                                                        tACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTaa                                                                        >  1:289947/1‑289 (MQ=255)                                                                                                                                                                                                                             cgacgaTGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGCTTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTcc                                                                                                                                                                          >  2:360075/1‑186 (MQ=255)                                                                                                                                                                                                                             cgacgaTGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTcc                                                                                                                                                                          <  1:360075/186‑1 (MQ=255)                                                                                                                                                                                                                             cgacgaTGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTcc                                                                                                                                                                          >  2:399797/1‑186 (MQ=255)                                                                                                                                                                                                                             cgacgaTGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTcc                                                                                                                                                                          <  1:399797/186‑1 (MQ=255)                                                                                                                                                                                                                             cgacgaTGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTAAAAgc                                                                    >  1:284814/1‑288 (MQ=255)                                                                                                                                                                                                                                          cTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTAAAAGCTCGCTCATCGCaga                                                      >  2:473478/1‑289 (MQ=255)                                                                                                                                                                                                                                                  tGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTAAAAGCTCGCTCATCGAAGAGACaaaaa                                              >  1:364924/1‑289 (MQ=255)                                                                                                                                                                                                                                                         gtctggGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCca                                                                                                  >  2:304782/1‑230 (MQ=255)                                                                                                                                                                                                                                                         gtctggGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCca                                                                                                  <  1:304782/230‑1 (MQ=255)                                                                                                                                                                                                                                                         gtctggGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTAAAAGCTCGCTCATCGAAGAGACAAAAAAGGt                                          >  1:133843/1‑286 (MQ=255)                                                                                                                                                                                                                                                         gtctggGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTAAAAGCTCGCTCATCGAAGAGACAAAAAAGGTGc                                        <  2:284814/288‑1 (MQ=255)                                                                                                                                                                                                                                                          tctggGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTAAAAGCTCGCTCATCGAAGAGACAAAAAAGGTGCTg                                      <  1:410023/289‑1 (MQ=255)                                                                                                                                                                                                                                                                   aTGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTAAAAGCTCGCTCATCGAAGAGACAAAAAAGGTGCTGAAACTTa                               >  1:178792/1‑287 (MQ=255)                                                                                                                                                                                                                                                                    tGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTAAAAGCTCGCTCATCGAAGAGACAAAAAAGGTGCTGAAACTTAACt                            <  1:473478/289‑1 (MQ=255)                                                                                                                                                                                                                                                                        ggTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTAAAAGCTCGCTCATCGAAGAGACAAAAAAGGTGCTGAAACTTAACTCAc                         <  2:364924/288‑1 (MQ=255)                                                                                                                                                                                                                                                                                           tGGAGAATGTCAGCTAAGCACAATATCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTAAAAGCTCGCTCATCGAAGAGACAAAAAAGGTGCTGAAACTTAACTCACTGTGGCGCGTCAGGGTGCg      <  2:289947/288‑1 (MQ=255)                                                                                                                                                                                                                                                                                               gaATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTAAAAGCTCGCTCATCGAAGAGACAAAAAAGGTGCTGAAACTTAACTCACTGTGGCGCGTCAGGGTGCGATTa  <  2:133843/288‑1 (MQ=255)                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                              TTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACATCTCCAGAGCTGGTATTGCCGTCAGCTGCGTTCTGCACTGTTATTTCATGAGCCACGGATTGCGGCGTTGCAGGTTAATTTGAAAGAAGCGTACTGCCACACTCTGGCGATCAGCCTT‑‑GAGATAATGCTGTACCACGACGATGAACCGCTAACGTTTGATCTGGTCTGGGATAATGGCGGTTGGCGCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTAAAAGCTCGCTCATCGAAGAGACAAAAAAGGTGCTGAAACTTAACTCACTGTGGCGCGTCAGGGTGCGATTA  >  NC_000913/4557089‑4557660

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑