Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	902,152	T→A	6.7%	F181I (TTC→ATC)	trxB →	thioredoxin reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	902,152	0	T	A	6.7%	183.7 / 3.1	60	F181I (TTC→ATC)	trxB	thioredoxin reductase
Reads supporting (aligned to +/- strand):  ref base T (28/28);  new base A (2/2);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GGTGTTTCGGCCTGTGCCACCTGTGACGGTTTCTTCTACCGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGACAACATGGGTGTAACCGGTGCGCGCCTGAAGAACAACGACGGCA  >  NC_002947/902005‑902296                                                                                                                                                    |                                                                                                                                                 ggTGTTTCGGCCTGTGCCACCTGTGACGGTTTCTTCTACCGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCtt                                                                                                                                                 <  1:372759/149‑1 (MQ=255)  gtgtTTCGGCCTGTGCCACCTGTGACGGTTTCTTCTACCGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTc                                                                                                                                                <  7:351180/149‑1 (MQ=255)  gtgtTTCGGCCTGTGCCACCTGTGACGGTTTCTTCTACCGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTc                                                                                                                                                <  2:73675/149‑1 (MQ=255)           ccTGTGCCACCTGTGACGGTTTCTTCTACCGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGcc                                                                                                                                          >  3:393591/1‑146 (MQ=255)           ccTGTGCCACCTGTGACGGTTTCTTCTACCGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGcc                                                                                                                                          >  2:181029/1‑146 (MQ=255)           ccTGTGCCACCTGTGACGGTTTCTTCTACCGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCgag                                                                                                                                       <  4:393591/149‑1 (MQ=255)           ccTGTGCCACCTGTGACGGTTTCTTCTACCGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCATCCGCGcc                                                                                                                                          >  4:202579/1‑146 (MQ=255)           ccTGTGCCACCTGTGACGGTTTCTTCTACCGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCATCCGCGcc                                                                                                                                          >  4:196430/1‑146 (MQ=255)                  cACCTGTGACGGTTTCTTCTACCGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATcc                                                                                                                                <  5:167862/149‑1 (MQ=255)                        tgACGGTTTCTTCTACCGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTc                                                                                                                           <  2:274280/148‑1 (MQ=255)                          aCGGTTTCTTCTACCGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGa                                                                                                                         >  3:159229/1‑148 (MQ=255)                          aCGGTTTCTTCTACCGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGAc                                                                                                                        <  2:386076/149‑1 (MQ=255)                            ggTTTCTTCTACCGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACaa                                                                                                                      >  5:80078/1‑149 (MQ=255)                                  ttctACCGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGaa                                                                                                                >  6:204369/1‑149 (MQ=255)                                  ttctACCGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGaa                                                                                                                >  2:266931/1‑149 (MQ=255)                                    ctACCGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGaa                                                                                                                <  5:204369/147‑1 (MQ=255)                                       ccGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGccc                                                                                                           <  2:361303/149‑1 (MQ=255)                                                gCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGcc                                                                                                   >  2:391417/1‑148 (MQ=255)                                                   ggTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAg                                                                                               <  1:266931/149‑1 (MQ=255)                                                         cGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGc                                                                                             >  8:93277/1‑145 (MQ=255)                                                           tGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGc                                                                                             >  1:270083/1‑143 (MQ=255)                                                             gTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATc                                                                                       >  6:342382/1‑147 (MQ=255)                                                                    gtggtAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGagctcaag                                                                              >  7:372994/1‑149 (MQ=255)                                                                          acacTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAAc                                                                        >  8:249682/1‑149 (MQ=255)                                                                            acTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGc                                                                      >  6:186044/1‑149 (MQ=255)                                                                                ccGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCAcc                                                                  >  8:171736/1‑149 (MQ=255)                                                                                  gTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCt                                                                <  6:80078/149‑1 (MQ=255)                                                                                    ggAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTg                                                               >  1:167031/1‑148 (MQ=255)                                                                                           gCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGt                                                       >  3:390544/1‑149 (MQ=255)                                                                                               tGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTg                                                   <  1:391417/149‑1 (MQ=255)                                                                                               tGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTg                                                   <  7:171736/149‑1 (MQ=255)                                                                                                 tACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTTCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCTAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTgg                                                  <  1:367228/148‑1 (MQ=255)                                                                                                 tACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTggg                                                 <  2:270083/149‑1 (MQ=255)                                                                                                 tACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTgg                                                  <  7:93277/148‑1 (MQ=255)                                                                                                 tACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCATCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTgg                                                  <  3:196430/148‑1 (MQ=255)                                                                                                 tACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCATCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTgg                                                  <  3:202579/148‑1 (MQ=255)                                                                                                   ccTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGt                                                <  4:390544/148‑1 (MQ=255)                                                                                                     tGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGa                                              >  3:260073/1‑148 (MQ=255)                                                                                                         cAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAgctcaagctcaa                                                                         <  8:20362/117‑1 (MQ=255)                                                                                                         cAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAgctcaagctcaa                                                                         >  7:20362/1‑117 (MQ=255)                                                                                                         cAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGacaaca                                         >  3:317719/1‑149 (MQ=255)                                                                                                          aaCATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGACAACAt                                        >  6:5756/1‑149 (MQ=255)                                                                                                           aCATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGACAACATg                                       <  4:317719/149‑1 (MQ=255)                                                                                                                 ccAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTgg                                                  >  8:230895/1‑132 (MQ=255)                                                                                                                 ccAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGACAACATGGgtgt                                  >  8:34237/1‑148 (MQ=255)                                                                                                                 ccAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGAAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTgg                                                  <  7:230895/132‑1 (MQ=255)                                                                                                                  cAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGACAACATGGGTGTaa                                >  3:81635/1‑149 (MQ=255)                                                                                                                         gTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGACAACATGGGTGTAACCGGTgc                         <  8:372994/149‑1 (MQ=255)                                                                                                                             cccTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAATATCCTTGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGACAACATGGGTGTAACCGGTgcgcgc                     <  1:186664/149‑1 (MQ=255)                                                                                                                               cTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGACAACATGGGTGTAACCGGTGCGCGCCt                   <  5:342382/149‑1 (MQ=255)                                                                                                                                     cacCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAg                                                                                               <  7:308555/67‑1 (MQ=255)                                                                                                                                     cacCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAg                                                                                               >  8:308555/1‑67 (MQ=255)                                                                                                                                       ccGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGACAACATGGGTGTAACCGGTGCGCGCCTGAAGaaca           <  7:249682/149‑1 (MQ=255)                                                                                                                                        cGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGac                                             <  4:88895/114‑1 (MQ=255)                                                                                                                                        cGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGac                                             >  3:88895/1‑114 (MQ=255)                                                                                                                                        cGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGACAACATGGGTGTAACCGGTGCGCGCCTGAAGaacaa          >  7:383984/1‑149 (MQ=255)                                                                                                                                         gTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGACAACATGGGTGTAACCGGTGCGCGCCTGAAGaacaac         >  2:179013/1‑149 (MQ=255)                                                                                                                                             cgAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGACAACATGGGTGTAACCGGTGCGCGCCTgaaga              <  5:425949/140‑1 (MQ=255)                                                                                                                                             cgAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGACAACATGGGTGTAACCGGTGCGCGCCTgaaga              >  6:425949/1‑140 (MQ=255)                                                                                                                                             cgAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGACAACATGGGTGTAACCGGTGCGCGCCTGAAGAACAacgacg     <  7:34237/149‑1 (MQ=255)                                                                                                                                                aaCCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGACAACATGGGTGTAAGCGGTGCGGGCCTGAAGAACAACGACGGCa  <  4:159229/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 GGTGTTTCGGCCTGTGCCACCTGTGACGGTTTCTTCTACCGCAACAAGCCGGTCGCCGTGGTCGGTGGTGGTAACACTGCCGTGGAAGAGGCACTGTACCTGGCCAACATCGCCAGCAAAGTGACCCTGGTACACCGTCGCGAAACCTTCCGCGCCGAGAAGATCCTGGTCGACAAGCTGAACGCCCGTGTGGCCGAAGGCAAGATCGAGCTCAAGCTCAACGCCACCCTGGACGAAGTGCTGGGTGACAACATGGGTGTAACCGGTGCGCGCCTGAAGAACAACGACGGCA  >  NC_002947/902005‑902296

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: ‑