Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,074,034	T→C	6.2%	V99A (GTG→GCG)	yihV →	6‑deoxy‑6‑sulfofructose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,074,034	0	T	C	6.2%	186.7 / 3.0	64	V99A (GTG→GCG)	yihV	6‑deoxy‑6‑sulfofructose kinase
Reads supporting (aligned to +/- strand):  ref base T (28/32);  new base C (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

AAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCAC  >  NC_000913/4073885‑4074181                                                                                                                                                      |                                                                                                                                                    aaGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCAtggt                                                                                                                                                     <  2:1085614/150‑1 (MQ=255)   ggCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCAtggtgg                                                                                                                                                   >  1:114361/1‑150 (MQ=255)   ggCGGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGGGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCAtggtgg                                                                                                                                                   <  1:1519899/150‑1 (MQ=255)     cTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGAt                                                                                                                                                 <  2:1510597/150‑1 (MQ=255)       gggggCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCACCATGATGGCGGAttg                                                                                                                                               <  2:6666/150‑3 (MQ=255)        ggggCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCACCATGATGGCGGAttgc                                                                                                                                              <  1:356958/150‑4 (MQ=255)          ggCGCAGGTCGATTTTATTGGTCGCATAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATAc                                                                                                                                               >  1:974310/1‑147 (MQ=255)          ggCGCAGGTCGATTTTATTGGTCGCATAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATAc                                                                                                                                               <  2:974310/147‑1 (MQ=255)            cgcAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAg                                                                                                                                          >  1:976644/1‑150 (MQ=255)               aGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACGAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCg                                                                                                                                       <  2:1393917/150‑1 (MQ=255)                   cGATTTTATTGGTCGCGTAGGTGATGACGAGACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACATCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCg                                                                                                                                   <  2:1383777/150‑1 (MQ=255)                   cGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCg                                                                                                                                   <  2:550120/150‑1 (MQ=255)                   cGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCg                                                                                                                                   <  1:495070/150‑1 (MQ=255)                    gATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAAACAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCgg                                                                                                                                  >  2:163290/1‑150 (MQ=255)                      ttttATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGAt                                                                                                                                <  2:494454/150‑1 (MQ=255)                         tATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAAt                                                                                                                             <  2:1590188/150‑1 (MQ=255)                             ggTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCAtt                                                                                                                         <  2:1503396/150‑1 (MQ=255)                              gTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCCGGCAGAGCTAGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTa                                                                                                                        <  2:1587089/150‑1 (MQ=255)                                      ggTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTAcccc                                                                                                                >  2:627092/1‑150 (MQ=255)                                          atgaCGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGcc                                                                                                            <  1:644849/150‑1 (MQ=255)                                            gacgacACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCg                                                                                                          >  1:1573813/1‑150 (MQ=255)                                                 cacGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATct                                                                                                     <  2:1623952/150‑1 (MQ=255)                                                 cacGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATct                                                                                                     >  1:1138932/1‑150 (MQ=255)                                                      gCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCtgctgc                                                                                                >  2:1254970/1‑150 (MQ=255)                                                              ctgctgGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCa                                                                                        <  2:1414915/150‑1 (MQ=255)                                                                   ggCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTg                                                                                   >  2:1133256/1‑150 (MQ=255)                                                                    gCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTgg                                                                                  >  2:78682/1‑150 (MQ=255)                                                                        agCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGCCGCAGAGTGGTTgg                                                                              >  2:1433631/1‑150 (MQ=255)                                                                          cTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTggag                                                                            <  1:929244/150‑1 (MQ=255)                                                                                  cTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGa                                                                    <  2:347239/150‑1 (MQ=255)                                                                                   tGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGAt                                                                   <  1:424010/150‑1 (MQ=255)                                                                                    ggggCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGAtt                                                                  >  1:101473/1‑150 (MQ=255)                                                                                        cGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTtctc                                                              >  1:1034949/1‑150 (MQ=255)                                                                                            aaCACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTa                                                                                                                        >  1:1444264/1‑88 (MQ=255)                                                                                            aaCACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTa                                                                                                                        <  2:1444264/88‑1 (MQ=255)                                                                                            aaCACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCCGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAg                                                          >  2:968227/1‑150 (MQ=255)                                                                                            aaCACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAg                                                          >  2:412076/1‑150 (MQ=255)                                                                                             acacCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGt                                                         <  2:1371986/150‑1 (MQ=255)                                                                                              cacCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTg                                                        >  2:1409089/1‑150 (MQ=255)                                                                                              cacCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTg                                                        <  2:475903/150‑1 (MQ=255)                                                                                                     tACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGAtgtt                                                 >  1:947857/1‑150 (MQ=255)                                                                                                      acacCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGAtgttg                                                >  2:433318/1‑150 (MQ=255)                                                                                                          cAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCt                                            <  2:1189445/150‑1 (MQ=255)                                                                                                          cAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCt                                            <  2:964659/150‑1 (MQ=255)                                                                                                            aaCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTgg                                          >  2:649486/1‑150 (MQ=255)                                                                                                            aaCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTgg                                          <  2:520006/150‑1 (MQ=255)                                                                                                             aCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGc                                         >  2:1436268/1‑150 (MQ=255)                                                                                                             aCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGc                                         <  2:1316119/150‑1 (MQ=255)                                                                                                             aCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGAGACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGc                                         <  1:751956/150‑1 (MQ=255)                                                                                                              cGGAATAACCAGGCGAAATCTTCGCAATCCGCCATCATGCTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCa                                        >  2:1193195/1‑150 (MQ=255)                                                                                                                       cAGGCGAAATCTTCGCAACCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGc                               >  1:561343/1‑150 (MQ=255)                                                                                                                         ggCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTg                             >  1:811918/1‑150 (MQ=255)                                                                                                                         ggCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTg                             <  1:1160661/150‑1 (MQ=255)                                                                                                                          gCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTgg                            >  2:613634/1‑150 (MQ=255)                                                                                                                          gCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTgg                            <  1:1266654/150‑1 (MQ=255)                                                                                                                           cGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTgcc                           >  2:946365/1‑148 (MQ=255)                                                                                                                                 cTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCacgacg                     <  1:1361931/150‑1 (MQ=255)                                                                                                                                     gCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCa                                        <  1:1487069/127‑1 (MQ=255)                                                                                                                                     gCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCa                                        >  2:1487069/1‑127 (MQ=255)                                                                                                                                      gtaTCCACCATGATGGCGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCt                >  1:1496030/3‑150 (MQ=255)                                                                                                                                      gtaTCCACCATGATGGCGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCt                >  2:1004132/3‑150 (MQ=255)                                                                                                                                       aaTCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTa               <  2:1138932/150‑1 (MQ=255)                                                                                                                                               atcatGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGcc       <  2:1573813/150‑1 (MQ=255)                                                                                                                                               atcatGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGcc       <  2:947857/150‑1 (MQ=255)                                                                                                                                               atcatGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGcc       <  2:1161751/150‑1 (MQ=255)                                                                                                                                                  atGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCtta    >  1:908823/1‑149 (MQ=255)                                                                                                                                                  atGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTc    >  1:1062363/1‑150 (MQ=255)                                                                                                                                                    ggtggATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCAc  <  2:1262490/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    AAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCAC  >  NC_000913/4073885‑4074181

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