Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,936,823	A→G	5.1%	M58V (ATG→GTG)	yebK →	putative DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,936,823	0	A	G	5.1%	261.0 / 3.4	79	M58V (ATG→GTG)	yebK	putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (34/41);  new base G (2/2);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

TGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAGGGAAAATATTTGAGTCCGCAATGGCAACGCTT  >  NC_000913/1936683‑1936969                                                                                                                                             |                                                                                                                                                   tGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACAcg                                                                                                                                            >  2:91018/1‑149 (MQ=255)       aTTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGt                                                                                                                                      >  2:582184/1‑149 (MQ=255)           gAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTc                                                                                                                                  <  4:279283/149‑1 (MQ=255)           gAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTc                                                                                                                                  >  1:57398/1‑149 (MQ=255)            aGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTcc                                                                                                                                 <  2:180207/149‑1 (MQ=255)            aGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTcc                                                                                                                                 <  3:416381/149‑1 (MQ=255)             gCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCt                                                                                                                                >  1:533513/1‑149 (MQ=255)              cAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTg                                                                                                                               >  4:155271/1‑149 (MQ=255)              cAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTg                                                                                                                               >  3:424503/1‑149 (MQ=255)               aaaTCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGa                                                                                                                              >  2:291746/1‑149 (MQ=255)               aaaTCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGa                                                                                                                              <  1:354948/149‑1 (MQ=255)               aaaTCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGa                                                                                                                              <  2:228581/149‑1 (MQ=255)                aagcatagCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATGCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGAt                                                                                                                             >  4:351125/7‑149 (MQ=255)                   cAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGAtttt                                                                                                                          >  1:481853/1‑149 (MQ=255)                      agCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTaaa                                                                                                                       <  3:255169/149‑1 (MQ=255)                      agCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTaaa                                                                                                                       <  2:403310/149‑1 (MQ=255)                      agCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTaaa                                                                                                                       <  2:121019/149‑1 (MQ=255)                      agCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTaaa                                                                                                                       <  3:300773/149‑1 (MQ=255)                      agCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTaaa                                                                                                                       <  1:17775/149‑1 (MQ=255)                                tGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTgg                                                                                                             <  1:321740/149‑1 (MQ=255)                                tGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTgg                                                                                                             <  1:291746/149‑1 (MQ=255)                                   cGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGcac                                                                                                          <  4:65907/149‑1 (MQ=255)                                          aTTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACaga                                                                                                       >  1:121486/1‑145 (MQ=255)                                          aTTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACaga                                                                                                       <  2:121486/145‑1 (MQ=255)                                           ttCTGGCTTCGCCCGTTAACGCGAGCCATTTGAGTCTTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGGCGGTGAATCGCTTCTGTCGCAGCATGGACACGCGCGGCTTTCCTAATAGTAAACTTCATCTGGCACAGAGTCTg                                                                                                  >  4:351117/1‑149 (MQ=255)                                           ttCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTg                                                                                                  >  2:585884/1‑149 (MQ=255)                                                 cTTCGCCCGATCACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAAt                                                                                            <  3:367923/149‑1 (MQ=255)                                                  ttCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATg                                                                                           <  1:214019/149‑1 (MQ=255)                                                     gCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGACa                                                                                        >  2:262621/1‑149 (MQ=255)                                                       ccGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCGTGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACt                                                                                      >  3:296425/1‑149 (MQ=255)                                                       ccGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCGTGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACt                                                                                      >  3:56986/1‑149 (MQ=255)                                                       ccGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACt                                                                                      <  1:438750/149‑1 (MQ=255)                                                       ccGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACt                                                                                      >  1:19714/1‑149 (MQ=255)                                                         gATAACGCTATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTATGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTcc                                                                                    <  3:351117/149‑1 (MQ=255)                                                         gATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTcc                                                                                    >  1:241729/1‑149 (MQ=255)                                                         gATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTcc                                                                                    >  4:564274/1‑149 (MQ=255)                                                         gATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTcc                                                                                    >  1:497401/1‑149 (MQ=255)                                                          aTAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTccc                                                                                   >  4:286529/1‑149 (MQ=255)                                                          aTAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTccc                                                                                   >  4:232806/1‑149 (MQ=255)                                                            aaCGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTa                                                                                 <  3:426812/149‑1 (MQ=255)                                                               gcgATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGt                                                                              <  2:481853/149‑1 (MQ=255)                                                                            gTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTc                                                                 >  4:272298/1‑149 (MQ=255)                                                                            gTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTc                                                                 >  3:321467/1‑149 (MQ=255)                                                                            gTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTc                                                                 >  3:500287/1‑149 (MQ=255)                                                                            gTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTc                                                                 <  2:10321/149‑1 (MQ=255)                                                                               ttGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAAt                                                              <  2:237098/149‑1 (MQ=255)                                                                               ttGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAAt                                                              <  4:282090/149‑1 (MQ=255)                                                                                    gctATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGa                                                         <  4:451713/149‑1 (MQ=255)                                                                                      tATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCaa                                                                     >  2:100139/1‑135 (MQ=255)                                                                                      tATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCaa                                                                     <  1:100139/135‑1 (MQ=255)                                                                                         ggCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACa                                                    >  3:284813/1‑149 (MQ=255)                                                                                            aCTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCg                                                 >  3:163185/1‑149 (MQ=255)                                                                                             cTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGt                                                >  2:100024/1‑149 (MQ=255)                                                                                              tGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGtt                                               <  4:284813/149‑1 (MQ=255)                                                                                              tGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGtt                                               <  2:433909/149‑1 (MQ=255)                                                                                                 aaGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGaa                                            <  1:100024/149‑1 (MQ=255)                                                                                                   gCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATc                                          >  4:541963/1‑149 (MQ=255)                                                                                                   gCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATc                                          <  2:370786/149‑1 (MQ=255)                                                                                                        tGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATaca                                     >  2:238984/1‑149 (MQ=255)                                                                                                           tAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAg                                  <  4:335496/149‑1 (MQ=255)                                                                                                           tAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAg                                  <  2:258222/149‑1 (MQ=255)                                                                                                            aGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTc                                                                 >  3:358919/1‑117 (MQ=255)                                                                                                            aGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTc                                                                 <  4:358919/117‑1 (MQ=255)                                                                                                             gCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAggg                                <  4:366985/149‑1 (MQ=255)                                                                                                             gCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAggg                                >  2:480521/1‑149 (MQ=255)                                                                                                             gCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAggg                                >  1:251543/1‑149 (MQ=255)                                                                                                                     aCGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAGGGAAAATAtt                        <  3:155271/149‑1 (MQ=255)                                                                                                                      cGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAGGGAAAATAttt                       <  3:212570/149‑1 (MQ=255)                                                                                                                        gTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAGGGAAAATATTTGa                     <  4:163185/149‑1 (MQ=255)                                                                                                                        gTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAGGGAAAATATTTGa                     <  3:286529/149‑1 (MQ=255)                                                                                                                        gTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAGGGAAAATATTTGa                     <  3:232806/149‑1 (MQ=255)                                                                                                                         tGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATa                                       <  3:551956/130‑1 (MQ=255)                                                                                                                         tGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATa                                       >  4:551956/1‑130 (MQ=255)                                                                                                                             tCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAGGGAAAATATTTGAGTCCg                >  3:47173/1‑149 (MQ=255)                                                                                                                             tCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAGGGAAAATATTTGAGTCCg                >  1:471998/1‑149 (MQ=255)                                                                                                                                  tCTGTCGCAGCGTGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAGGGAAAATATTTGAGTCCGCAATg           <  4:296425/149‑1 (MQ=255)                                                                                                                                  tCTGTCGCAGCGTGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAGGGAAAATATTTGAGTCCGCAATg           <  4:56986/149‑1 (MQ=255)                                                                                                                                       cGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAGGGAAAATATTTGAGTCCGCAATGGCAAc      <  2:533513/149‑1 (MQ=255)                                                                                                                                       cGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAGGGAAAATATTTGAGTCCGCAATGGCAAc      <  1:480521/149‑1 (MQ=255)                                                                                                                                         cagcaTGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAGGGAAAATATTTGAGTCCGCAATGGCAACGc    <  4:75368/149‑1 (MQ=255)                                                                                                                                           gcaTGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAGGGAAAATATTTGAGTCCGCAATGGCAACGCtt  >  2:371404/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                   TGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCATACACAGGGAAAATATTTGAGTCCGCAATGGCAACGCTT  >  NC_000913/1936683‑1936969

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