Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	243,632	T→C	6.1%	I30I (ATT→ATC)	gmhA →	D‑sedoheptulose 7‑phosphate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	243,632	0	T	C	6.1%	206.8 / 3.8	65	I30I (ATT→ATC)	gmhA	D‑sedoheptulose 7‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base T (27/34);  new base C (2/2);  total (29/36)
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.25e-01

GCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGA  >  NC_000913/243488‑243780                                                                                                                                                 |                                                                                                                                                     gCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAgc                                                                                                                                                  <  4:84580/149‑1 (MQ=255) gCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAgc                                                                                                                                                  <  4:17400/149‑1 (MQ=255) gCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAgc                                                                                                                                                  <  2:451239/149‑1 (MQ=255)      aTCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGcgcgcg                                                                                                                                             <  1:538125/149‑1 (MQ=255)       tCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCgcgg                                                                                                                                            >  4:567442/1‑149 (MQ=255)        ccGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCgcggc                                                                                                                                           >  3:284938/1‑149 (MQ=255)          ggTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGcggcgg                                                                                                                                         <  1:153830/149‑1 (MQ=255)             acacTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTcc                                                                                                                                      <  4:265523/149‑1 (MQ=255)                 tGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGtt                                                                                                                                  >  3:524249/1‑149 (MQ=255)                 tGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGtt                                                                                                                                  >  1:380704/1‑149 (MQ=255)                 tGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGtt                                                                                                                                  >  1:114782/1‑149 (MQ=255)                 tGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGtt                                                                                                                                  <  2:154867/149‑1 (MQ=255)                  gCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTa                                                                                                                                 <  2:67849/149‑1 (MQ=255)                       ttGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTGACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGa                                                                                                                            >  1:375599/1‑149 (MQ=255)                       ttGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGa                                                                                                                            <  4:223283/149‑1 (MQ=255)                            tATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCt                                                                                                                       >  2:565710/1‑149 (MQ=255)                              ttACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCttt                                                                                                                     <  2:468477/149‑1 (MQ=255)                                 cATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTaaa                                                                                                                  <  3:264033/149‑1 (MQ=255)                                   tttATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGc                                                                                                                <  3:183425/149‑1 (MQ=255)                                   tttATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATCCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGc                                                                                                                >  2:274632/1‑149 (MQ=255)                                   tttATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATCCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGc                                                                                                                >  2:96384/1‑149 (MQ=255)                                     tATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCg                                                                                                              >  3:316215/1‑149 (MQ=255)                                      aTGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCgg                                                                                                             <  3:97231/149‑1 (MQ=255)                                       tGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATCCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGt                                                                                                            <  1:96384/149‑1 (MQ=255)                                       tGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATCCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGt                                                                                                            <  1:274632/149‑1 (MQ=255)                                        gCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTg                                                                                                           <  1:559666/149‑1 (MQ=255)                                            aaGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCaa                                                                                                       >  2:452374/1‑149 (MQ=255)                                              ggATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGa                                                                                                                            <  4:399628/126‑1 (MQ=255)                                              ggATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGa                                                                                                                            >  3:399628/1‑126 (MQ=255)                                              ggATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAg                                                                                                     >  2:397702/1‑149 (MQ=255)                                               gATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGTGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGt                                                                                                    >  4:83721/1‑149 (MQ=255)                                                    ccTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCttt                                                                                               >  2:184122/1‑149 (MQ=255)                                                       cATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCt                                                                                            <  2:309790/149‑1 (MQ=255)                                                         tGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGc                                                                                          <  4:494512/149‑1 (MQ=255)                                                         tGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGc                                                                                          <  2:229035/149‑1 (MQ=255)                                                          gTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCg                                                                                         <  3:35677/149‑1 (MQ=255)                                                              cAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCaa                                                                                     <  3:567442/149‑1 (MQ=255)                                                                  aTCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAAcggc                                                                                 >  1:486770/1‑149 (MQ=255)                                                                    cTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCCATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAAcggcgg                                                                               <  3:348209/149‑1 (MQ=255)                                                                     ttATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGt                                                                              >  1:325337/1‑149 (MQ=255)                                                                     ttATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGt                                                                              >  4:451908/1‑149 (MQ=255)                                                                       aTTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGACGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTc                                                                            <  1:529503/149‑1 (MQ=255)                                                                         tCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTccc                                                                          <  3:282986/149‑1 (MQ=255)                                                                          cGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCa                                                                         <  2:380704/149‑1 (MQ=255)                                                                                                    aaCGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAag                                               >  2:555321/1‑149 (MQ=255)                                                                                                         tGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTg                                          <  3:323013/149‑1 (MQ=255)                                                                                                            cTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGAcc                                       >  4:287317/1‑149 (MQ=255)                                                                                                            cTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGAcc                                       <  2:486770/149‑1 (MQ=255)                                                                                                              aaCTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCgg                                     >  3:259392/1‑149 (MQ=255)                                                                                                                cTTTTTAAAAGATAACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTc                                   <  4:358939/149‑1 (MQ=255)                                                                                                                 tttttAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTc                                                                                                                                       >  2:523845/1‑48 (MQ=255)                                                                                                                 tttttAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTc                                                                                                                                       <  1:523845/48‑1 (MQ=255)                                                                                                                 tttttAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGAGCGGTCg                                  >  2:217580/1‑149 (MQ=255)                                                                                                                 tttttAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTCg                                  >  4:175878/1‑149 (MQ=255)                                                                                                                 tttttAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTCg                                  >  3:202863/1‑149 (MQ=255)                                                                                                                     tAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTCGCTAc                              >  1:246634/1‑149 (MQ=255)                                                                                                                       aaaGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTCGCTACCg                            <  1:268994/149‑1 (MQ=255)                                                                                                                       aaaGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTCGCTACCg                            <  3:18235/149‑1 (MQ=255)                                                                                                                       aaaGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTCGCTACCg                            <  4:316215/149‑1 (MQ=255)                                                                                                                          gATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTCGCTACCGTGa                         <  4:264488/149‑1 (MQ=255)                                                                                                                                  cAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTCGCTACCGTGAAAACCGTc                 >  1:149785/1‑149 (MQ=255)                                                                                                                                   aaTATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTCGCTACCGTGAAAACCGTcc                <  4:284938/149‑1 (MQ=255)                                                                                                                                      atTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTCGCTACCGTGAAAACCGTCCggg             <  4:524249/149‑1 (MQ=255)                                                                                                                                       ttCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGc            >  2:498182/1‑149 (MQ=255)                                                                                                                                            gCCATTCAGCGCGCGGCGGTCCTGTCAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTAccc       >  2:81886/1‑149 (MQ=255)                                                                                                                                                 tCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGa  <  2:375599/149‑1 (MQ=255)                                                                                                                                                 tCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGa  <  1:503456/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     GCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGA  >  NC_000913/243488‑243780

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