Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,270,523	C→A	L96M (CTG→ATG)	yjbQ →	thiamine phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,270,523	0	C	A	98.9%	314.1 / NA	91	L96M (CTG→ATG)	yjbQ	thiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (47/43);  total (47/44)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCTCTGACACTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGAGTAAAAAATGACCATTTC  >  NC_000913/4270377‑4270668                                                                                                                                                   |                                                                                                                                                  ccTCTCTGACACTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAatgatg                                                                                                                                                 <  2:234843/149‑1 (MQ=255)  ctctctGACACTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATgg                                                                                                                                                <  4:233995/149‑1 (MQ=255)      ctGACACTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGaa                                                                                                                                             >  1:282137/1‑148 (MQ=255)       tGACACTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGa                                                                                                                                              <  2:21241/146‑1 (MQ=255)       tGACACTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACa                                                                                                                                           >  1:45468/1‑149 (MQ=255)           acTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGc                                                                                                                                       <  8:201773/149‑1 (MQ=255)           acTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGc                                                                                                                                       >  5:126784/1‑149 (MQ=255)           acTTAATGAGAACTGCGATCCCACAGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGCGGGAGCAGACCCTATGCCTTCTCATAGCAAATCCGCCATGATGGGAACagctc                                                                                                                                       >  8:326889/1‑145 (MQ=255)               aaTGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCtt                                                                                                                                     <  6:126784/147‑1 (MQ=255)               aaTGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTg                                                                                                                                    <  1:90493/148‑1 (MQ=255)               aaTGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTg                                                                                                                                    <  2:38657/148‑1 (MQ=255)                aTGAGAACTCCGATCCCACCGTACGCCACTAGATGGAGCGTTTTTTCCTCCTCAAGGTTCCCGACAACGGAATTTATGATCATGACTATGAGTGTGCACACTATATGCCTTCTCATTTCAAATCCTCAATGATTTGAACATAGCTTTTa                                                                                                                                  <  2:29254/149‑1 (MQ=255)                 tGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACa                                                                                                                                           <  4:119785/139‑1 (MQ=255)                 tGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACa                                                                                                                                           >  3:119785/1‑139 (MQ=255)                   agaACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTg                                                                                                                               >  7:238993/1‑149 (MQ=255)                       cTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCgg                                                                                                                           <  6:11056/149‑1 (MQ=255)                        tGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGt                                                                                                                          >  3:152955/1‑149 (MQ=255)                          cGATCCCACCGTACGCCACTACATGGAGCGTTTTTTCCTCCGCCCCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATGCTCAATGATGGGAACATCGCTTGTATTGCCGGTGc                                                                                                                        <  5:106421/149‑1 (MQ=255)                          cGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTg                                                                                                                         >  7:304666/1‑148 (MQ=255)                          cGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGc                                                                                                                        <  8:61124/149‑1 (MQ=255)                            aTCCCACCGTACGCCACTACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGTAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCAt                                                                                                                      <  3:42377/149‑1 (MQ=255)                                cACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAg                                                                                                                  <  4:268227/149‑1 (MQ=255)                                 aCCGTACGCCCCGACATGGAGCGTTTTTTCCTCCGCACTGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACCTCGCTTGTATTGCCGGGGCATAAAgg                                                                                                                 >  2:67808/1‑149 (MQ=255)                                 aCCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACTGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAgg                                                                                                                 >  2:72516/1‑149 (MQ=255)                                 aCCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAgg                                                                                                                 <  5:7258/149‑1 (MQ=255)                                   cGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGc                                                                                                               <  1:226005/149‑1 (MQ=255)                                   cGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGc                                                                                                               <  5:61228/149‑1 (MQ=255)                                   cGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGc                                                                                                               <  3:255557/149‑1 (MQ=255)                                   cGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTGCCCTACAACGTAAATTATTAGCATGACTATGAGGGAGCACACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGc                                                                                                               <  7:234029/149‑1 (MQ=255)                                      aCGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTa                                                                                                            >  6:177239/1‑149 (MQ=255)                                       cGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAATGGCGTAt                                                                                                           <  2:45468/149‑1 (MQ=255)                                        gCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGTTATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTAtt                                                                                                          <  8:228426/149‑1 (MQ=255)                                            cGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGt                                                                                                                                   >  8:278704/1‑120 (MQ=255)                                            cGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGt                                                                                                                                   <  7:278704/120‑1 (MQ=255)                                             gACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGAc                                                                                                     >  7:114950/1‑149 (MQ=255)                                                 tGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATTATGTGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGc                                                                                                 <  1:292300/149‑1 (MQ=255)                                                 tGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGATCATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCgg                                                                                                  >  3:22140/1‑148 (MQ=255)                                                   gAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATggg                                                                                                                                               >  3:234263/1‑101 (MQ=255)                                                   gAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATggg                                                                                                                                               <  4:234263/101‑1 (MQ=255)                                                   gAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCAc                                                                                               >  3:317741/1‑149 (MQ=255)                                                   gAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCAc                                                                                               >  5:155630/1‑149 (MQ=255)                                                      cGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATCTGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTg                                                                                            <  5:2425/149‑1 (MQ=255)                                                          ttttCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGc                                                                                          >  7:194593/1‑147 (MQ=255)                                                           tttCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCaa                                                                                        <  8:114950/148‑1 (MQ=255)                                                            ttCCTCCGCACCGTTCCCGACAACTGAAATTATGAGCATGACTTTGAGGGAGAAGACGATATGCCTTCTCATATCAAATCATCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAg                                                                                       <  6:284819/148‑1 (MQ=255)                                                              cctccGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGc                                                                                     <  6:123761/148‑1 (MQ=255)                                                                tccGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGTGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGc                                                                                                               <  4:193764/120‑1 (MQ=255)                                                                tccGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGTGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGc                                                                                                               >  3:193764/1‑120 (MQ=255)                                                                  cGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCAttt                                                                                 <  8:103371/148‑1 (MQ=255)                                                                   gCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTgg                                                                               >  7:85148/1‑149 (MQ=255)                                                                           cccGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTggggg                                                                        <  2:282137/148‑1 (MQ=255)                                                                               aCAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACAt                                                                   >  2:195970/1‑149 (MQ=255)                                                                                cAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATc                                                                  <  1:195970/149‑1 (MQ=255)                                                                                 aaCGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCg                                                                 >  7:88558/1‑149 (MQ=255)                                                                                   cGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGc                                                                >  8:275885/1‑148 (MQ=255)                                                                                    ggAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCAt                                                              >  1:72892/1‑149 (MQ=255)                                                                                    ggAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCAt                                                              >  6:122810/1‑149 (MQ=255)                                                                                     gAAATTATGAGCATGACTATGAGGGAGGAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGATCCTCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATc                                                             <  1:330180/149‑1 (MQ=255)                                                                                              aGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGACCATCGCTTGTATTGCCGGTGCATAAAGGGCGGATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGa                                                    >  5:248825/1‑149 (MQ=255)                                                                                                cATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATc                                                  >  7:172715/1‑149 (MQ=255)                                                                                                    aCTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGt                                              >  3:38624/1‑149 (MQ=255)                                                                                                        tGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGca                                          >  8:127299/1‑149 (MQ=255)                                                                                                           gggAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGcatc                                        >  8:163087/1‑148 (MQ=255)                                                                                                            ggAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGAGTATTCAGACCGGCACCTGGCAAGGCACTTGGCTGGGGGACCATCGCATCCACGGCGGATCGCGTCGcatcat                                      >  3:12949/1‑149 (MQ=255)                                                                                                             gAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCatcatc                                     <  1:272081/149‑1 (MQ=255)                                                                                                              aGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATcg                                    >  2:333908/1‑149 (MQ=255)                                                                                                               gCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATcgc                                   <  6:155630/149‑1 (MQ=255)                                                                                                               gCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATcgc                                   <  4:317741/149‑1 (MQ=255)                                                                                                                cAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGCCCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATcgcg                                  >  4:255166/1‑149 (MQ=255)                                                                                                                 aGACGCTATGCCTTCTCATATCAAATCCTCAATGATGGGCACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGCCCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa                                 <  1:67808/149‑1 (MQ=255)                                                                                                                 aGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa                                 <  1:72516/149‑1 (MQ=255)                                                                                                                   aCGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGaca                               >  8:33546/1‑149 (MQ=255)                                                                                                                       tatGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGCCCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGacac                              >  6:132540/1‑146 (MQ=255)                                                                                                                       tatGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGacac                              <  5:132540/146‑1 (MQ=255)                                                                                                                       tatGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTa                            >  3:52137/1‑148 (MQ=255)                                                                                                                          gCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAg                        >  6:130527/1‑149 (MQ=255)                                                                                                                          gCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAg                        >  3:59342/1‑149 (MQ=255)                                                                                                                           ccTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGCCACTACAAgg                       >  5:106912/1‑149 (MQ=255)                                                                                                                           ccTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAgg                       >  1:6137/1‑149 (MQ=255)                                                                                                                           ccTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAg                        <  6:311678/148‑1 (MQ=255)                                                                                                                            cTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAggg                      >  2:7907/1‑149 (MQ=255)                                                                                                                            cTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAggg                      >  8:322976/1‑149 (MQ=255)                                                                                                                               ctcATATCAAATCCTCAATGCTGGGAACATCGCTTGTATTGCCGGGGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGa                    <  1:17256/148‑1 (MQ=255)                                                                                                                                  atatCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGAGTa                 >  1:72192/1‑148 (MQ=255)                                                                                                                                  atatCAAATCCTCAATGATGGGAACATCGCATGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGAATTTGGCTGGGGGAACATCGCATCCACGGCCGATCGCGTCTCATCATCGCGACACTACAACGGGAGTa                 >  1:123892/1‑148 (MQ=255)                                                                                                                                   tatCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTTGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAggg                      >  7:82176/1‑142 (MQ=255)                                                                                                                                   tatCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTTGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAggg                      <  8:82176/142‑1 (MQ=255)                                                                                                                                        aaTCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGAGTAAAAAAt           >  6:42551/1‑148 (MQ=255)                                                                                                                                         aTCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGAGTAAAAAATGa         <  7:144276/149‑1 (MQ=255)                                                                                                                                            cTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGAGTAAAAAATGACCa      <  7:322976/149‑1 (MQ=255)                                                                                                                                            cTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGAGTAAAAAATGACCa      <  8:304666/149‑1 (MQ=255)                                                                                                                                             tCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGAGTAAAAAATGACCAt     >  2:37070/1‑149 (MQ=255)                                                                                                                                                atgatgGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGAGTAAAAAATGACCATTTc  <  4:152955/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  CCTCTCTGACACTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGAGTAAAAAATGACCATTTC  >  NC_000913/4270377‑4270668

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