Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,636,925	G→C	R77P (CGC→CCC)	creC →	sensory histidine kinase CreC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,636,925	0	G	C	98.5%	196.4 / NA	67	R77P (CGC→CCC)	creC	sensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (26/40);  minor base A (1/0);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

AGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGC  >  NC_000913/4636685‑4637170                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      aGGGGCCTGTAATGCGTCTCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTTGCGGAGCTGGCGCGGCCCGATTTGCTATCTGGGGACCCACCGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTcgc                                                                                                                                                                                                                                             >  2:593944/1‑251 (MQ=255)         gtaaggctgaTCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGTGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAat                                                                                                                                                                                                                                       <  2:225034/240‑1 (MQ=255)            aTGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCCGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCgg                                                                                                                                                                                                                                  >  1:310843/1‑251 (MQ=255)               cGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTGCGCGCCAATATCGGTgg                                                                                                                                                                                                                               <  1:465352/251‑1 (MQ=255)               cGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGAAGCCTGGTTCGTACTGGCCATTTTTGTCCAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCCACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATCTCGGTgg                                                                                                                                                                                                                               >  1:257182/1‑251 (MQ=255)                  aTCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGAGCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGTAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCAt                                                                                                                                                                                                                            <  2:396825/251‑1 (MQ=255)                    cGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGGTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGGCTGGGCCACTGGCGCAGGCGGTTAATACGCTACAACATCCCCCGTTTCGCGCCCATAACGGTGGCAtt                                                                                                                                                                                                                           >  2:175176/1‑250 (MQ=255)                            ggTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTg                                                                                                                                                                                                                  <  1:412570/251‑1 (MQ=255)                                      gggCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGa                                                                                                                                                                                                          <  1:230568/249‑1 (MQ=255)                                      gggCTATTTTTTACCGGTGGCGTGGGCAGACTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTTGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAat                                                                                                                                                                                                        <  2:451004/251‑1 (MQ=255)                                               tttACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCt                                                                                                                                                                                               >  1:450536/1‑251 (MQ=255)                                                     ggtggcggtggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTGAAACCGGGCATTCGAAGAGCAACGGAGGGGACGTTTATCGACACCGCAACGTTGCTGGCGGAGCTGGGGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGTAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATTAATATCATGTCTATATg                                                                                                                                                                                          <  2:554661/250‑1 (MQ=255)                                                       tggcggtggcAGCCTGGGTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTTGGGACCCAACGCATGGGCAACTGGCGCAGGGGTTTAATCAGCTACAACATCCCCCGTTTTGCGCCAATATCGGTGGCATTAACCAAGTGCGCAATGAATATCACGTCTATATGAcc                                                                                                                                                                                       >  2:179052/1‑251 (MQ=255)                                                        ggcggtggcAGCCTGGTTCGTACTGGCCATTTTGGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGAcc                                                                                                                                                                                       <  1:423372/250‑1 (MQ=255)                                                        ggcggtggcAGCCTGGGTCGTACTGGCCATTTTTTTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCg                                                                                                                                                                                      <  2:324027/251‑1 (MQ=255)                                                           ggtggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGAACGATg                                                                                                                                                                                   >  1:195239/1‑251 (MQ=255)                                                                     tGGTTCGGACTGGCCCTTTTTGTCAAAGAAGTTAAACCGGGCGTGCGTCGAGCAACCGAGGGGTCGTTGATCGACACCGCAACGCTGCTGGCGGAGCTGGCGCGTCACGATTTGCTCTCTGGGGACCCAACGCATGGGCCACTGGCGCAGGCGTTTAATCAGCTACAACCTCCCCCGTTTCGCGCCAATATAGGTGGCATTAACAAACTGCGCAATGAATATCATGGCTCTATGACCGATGCGCAGGGct                                                                                                                                                                          >  2:26724/1‑249 (MQ=255)                                                                      ggTTCGTACTGGCAATTTTTGTCCAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTCATCAGCTACAACATCCCCCGCTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCACTGAATATCATGTCTATATGACCGATGCGCAGGGCaa                                                                                                                                                                         >  2:343608/1‑250 (MQ=255)                                                                          cgtccTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGCTCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCGGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACAGCCCCCGTTTCGCGCCAATATCGGTGGCATTATCAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTAt                                                                                                                                                                    <  1:46781/247‑1 (MQ=255)                                                                          cGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAAGGGAGGGGACGTTGATCGCCACCGCACCGTTGCGGGCGGAGCTGGCGCGTCACGATGTGGTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTTGCAGTAACAAAGTTCGCAATGAATATCATGTCTATATTACCGATGCGCAGGTCAAAGTa                                                                                                                                                                     <  2:583276/250‑1 (MQ=255)                                                                                ggCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTc                                                                                                                                                               <  1:468086/250‑1 (MQ=255)                                                                                    aTTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTc                                                                                                                                                          >  1:8418/1‑251 (MQ=255)                                                                                    aTTTTTGTAAAAGAAGTTAAATCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGTTTTGCTCTCTGGGGCCCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCCGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACACAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGCATTGTTCGATTc                                                                                                                                                          >  2:281229/1‑251 (MQ=255)                                                                                     tttttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGTCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACGGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTGCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAAGGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGAtt                                                                                                                                                           <  1:86552/249‑1 (MQ=255)                                                                                         tGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCa                                                                                                                                                      >  1:205785/1‑250 (MQ=255)                                                                                          gTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCaaa                                                                                                                                                    <  1:212517/251‑1 (MQ=255)                                                                                          gTCAAAGAAGTGAAACCGGGCGTGCGAAGAGGAACGGAGGGGACGTTGAGCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCACTCGGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCAGTAACAAAGTGCTCAATGAATATCATGTCTATATGACCGATGCGCTGGGCAAAGTATTGTTCGATTCGGCaaa                                                                                                                                                    <  2:34772/251‑1 (MQ=255)                                                                                             aaaGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATaa                                                                                                                                                 <  1:593944/251‑1 (MQ=255)                                                                                                gaagTTAAACCGGGCGTACGAAGAGCAACGGAGGGGACGGTGAGCGACACCGCACCGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGTGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGc                                                                                                                                              <  2:154703/251‑1 (MQ=255)                                                                                                gaagTTAAACCGGGCGGGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGc                                                                                                                                              <  2:67040/251‑1 (MQ=255)                                                                                                       aaCCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTgg                                                                                                                                        >  1:204029/1‑250 (MQ=255)                                                                                                        aCCGGGGGTTCGAGGAGCAACGGAGGGGACGTTGTTCGACACCGCAGCGTTGCTGGCGGGGCTGGCGCGTCCCGATTTGCTCTCTGGGTACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGAGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTATATATGCCCGATGCGCAGGGCAAATTATTGTTCGAATCGGCAAATAAAGCCGTTGGAc                                                                                                                                      <  2:416251/251‑1 (MQ=255)                                                                                                        aCCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGAc                                                                                                                                      >  1:201477/1‑251 (MQ=255)                                                                                                         ccGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAATGTGGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGTGACCCAACGCATGGGCAACTGGCGCAGGGGTGTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAGTGAATATCATGTCTAGATGACCGTTGCGCAGGGCAAAGTATTTTGCGATTCgacgaataggtgcgttggaca                                                                                                                                     <  2:30038/251‑21 (MQ=255)                                                                                                              cGTGCGAAGAGCAACGGAGGGGACGTTGGTCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTACCGTGTTGCTCTCTGGGGACCCAACGCATGGGCAAGTGGCGCAGGCGTTTATTCAGCTACAACATCCCCCGTTTCGCGCCACGATCGGTGGCATTAACGAAGTGCGCAATGATTTTCAGGTCAATATGACCGATGCGAAGGGCCAAGTTTTGTTCGATTCGACAAAACAACCCGTTGGACAggttt                                                                                                                                >  2:422106/1‑248 (MQ=255)                                                                                                              cGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGatt                                                                                                                                <  1:175176/251‑1 (MQ=255)                                                                                                               gTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACCTCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGGATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGatt                                                                                                                                >  2:147293/1‑250 (MQ=255)                                                                                                                tGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGCCCTGTCGCGTCCCCATTTGCTCTATGGGGCCCCAACGCATGGGCACCTGGCGCAGGCGTTTAATCAGCTACAACATACCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGCGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAACGTATTGTTCGATTCGGCAAATCACGCCGTTGGACAGGattat                                                                                                                              >  1:345976/1‑251 (MQ=255)                                                                                                                   gAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTcg                                                                                                                           <  1:86712/251‑1 (MQ=255)                                                                                                                     atagCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTcgcg                                                                                                                         >  1:361174/3‑251 (MQ=255)                                                                                                                         cAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTgg                                                                                                                     <  1:324023/251‑1 (MQ=255)                                                                                                                                ggggACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGGCCCAACGCATGGGCACCTGGCGCAGGCGTTTAATCAGCTACAACATCACCCGTTTCGCGCCAATATCGGGGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGACTGACg                                                                                                              >  2:285565/1‑251 (MQ=255)                                                                                                                                   gACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCt                                                                                                           >  1:146797/1‑251 (MQ=255)                                                                                                                                          aTCGACCCCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTTGAATGACGTCTGGCTAAc                                                                                                    <  2:390225/251‑1 (MQ=255)                                                                                                                                                        gTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGt                                                                                      <  1:107765/251‑1 (MQ=255)                                                                                                                                                           gCTGGCGGAGCTGGCGCGTCTCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTTGTCAGTATg                                                                                   <  1:105287/251‑1 (MQ=255)                                                                                                                                                            cTGGCGAAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTGGGACAGGATTATTCGCGCTGGAAGGACGTATGGATAACGTTGCGTGGTCCGTATgg                                                                                  >  2:22324/1‑251 (MQ=255)                                                                                                                                                                       ggCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCAAGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCgcagca                                                                       <  1:422106/251‑1 (MQ=255)                                                                                                                                                                       ggCGCGGCCTGATTTGCTCTCTGGGTACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCTTTTCGGGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGAGGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCTGTTGGACAGGATTATTCGCGCGGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCgcagc                                                                        <  1:586390/250‑1 (MQ=255)                                                                                                                                                                        gcgcGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCAc                                                                      <  1:285565/251‑1 (MQ=255)                                                                                                                                                                              cccGATTTGCGCTCTGGGGACCCAACGCATGGGCAACTGGCGCCGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCAGTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGAGGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAGAGCCGTTGGACAGGATCATTCGCGCTGGAATGACGTCTGGTTAACGTTGCGTGGTCAGTATGGTGCTCGCAGCCCGTTGCa                                                                <  2:484570/251‑1 (MQ=255)                                                                                                                                                                                  attcgatcgctGGGGACCCAACGCATGGGCAACACGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGCATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATTGTGCGCGCAGCACGTTGCAAAAt                                                            <  2:28211/242‑1 (MQ=255)                                                                                                                                                                                                 aCCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGaa                                             <  1:137638/251‑1 (MQ=255)                                                                                                                                                                                                    cAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGt                                          <  1:320018/251‑1 (MQ=255)                                                                                                                                                                                                         cATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCtgt                                     >  1:460208/1‑251 (MQ=255)                                                                                                                                                                                                          aTGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCtgtg                                    <  1:281469/251‑1 (MQ=255)                                                                                                                                                                                                            gggCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACCAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAACGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCg                                                     >  2:355254/1‑232 (MQ=255)                                                                                                                                                                                                            gggCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCg                                                     <  1:355254/232‑1 (MQ=255)                                                                                                                                                                                                                          ggCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTGGCCGCAc                    <  1:147293/251‑1 (MQ=255)                                                                                                                                                                                                                           gCGTTTAATCAGCTACAACATCCCCCGGTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCAcc                   <  2:205785/251‑1 (MQ=255)                                                                                                                                                                                                                            cGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATCATTCGCGCTGGAATGACGACTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCg                  >  2:1864/1‑251 (MQ=255)                                                                                                                                                                                                                                tAATCGGCTACAACAGCCTCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCATGGAAGATATTGTCTATATGACCGATGCGCAGGGCAATGTGTTGTTCGATTCGGCAAATAAAGCAGGTGGACAGGAGTATTCGCGGTGGAATGGCGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATACTGCCGATCCCGAAAGGTCTTTGATTTATGTTGCCGCATCGAtt               <  2:394357/250‑1 (MQ=255)                                                                                                                                                                                                                                   tCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATgg           >  1:80821/1‑251 (MQ=255)                                                                                                                                                                                                                                       cTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACgg       >  1:580844/1‑251 (MQ=255)                                                                                                                                                                                                                                          caacaTCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGTATTATTCGCGCTGGAATGACGTCTGGCTAACTTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGatgtatgt                            <  1:530166/227‑1 (MQ=255)                                                                                                                                                                                                                                          caacaTCCCCCGTTTCGCGCCAATATAGGTGGCCTTACCAAAGTGCGCACTGAATATCCTGTATATATGAACGATGCGCCGGGCAAAGTCTCGTTCGATTCCGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGCCTGACGTCTGGCTAACGTTGCGTGGTTAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGCTCCCGCAAGTTCTGTGctgtctgtcgatcggacgcgcgtcgtgtaggg    >  2:530166/1‑227 (MQ=255)                                                                                                                                                                                                                                            acaTCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTcgc  <  2:195239/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      AGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGC  >  NC_000913/4636685‑4637170

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