Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,569,304	A→G	6.6%	N412D (AAC→GAC)	malQ →	4‑alpha‑glucanotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,569,304	0	A	G	6.6%	174.4 / 2.3	61	N412D (AAC→GAC)	malQ	4‑alpha‑glucanotransferase
Reads supporting (aligned to +/- strand):  ref base A (27/29);  new base G (2/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

CACGAAGCGGCCCGCAGCAACGGCATGGCCGTGGGCCTGATCGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGCCCACGCTGGCGGCCTGCGCATCGACCATGTGATGGGCCTGCGCCGG  >  NC_002947/4569160‑4569450                                                                                                                                                 |                                                                                                                                                      cACGAAGCGGCCCGCAGCAACGGCATGGCCGTGGGCCTGATCGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCa                                                                                                                                                   >  7:9098/1‑149 (MQ=255) cACGAAGCGGCCCGCAGCAACGGCATGGCCGTGGGCCTGATCGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCa                                                                                                                                                   >  7:9942/1‑149 (MQ=255)   cGAAGCGGCCCGCAGCAACGGCATGGCCGTGGGCCTGATCGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCGACCACg                                                                                                                                                 >  3:127659/1‑149 (MQ=255)   cGAAGCGGCCCGCAGCAACGGCATGGCCGTGGGCCTGATCGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCGACCACg                                                                                                                                                 >  3:132687/1‑149 (MQ=255)           cccGCAGCAACGGCATGGCCGTGGGCCTGATCGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAg                                                                                                                                         <  6:4774/149‑1 (MQ=255)           cccGCAGCAACGGCATGGCCGTGGGCCTGATCGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAg                                                                                                                                         <  3:292648/149‑1 (MQ=255)            ccGCAGCAACGGCATGGCCGTGGGCCTGATCGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAg                                                                                                                                         >  2:341834/1‑148 (MQ=255)                 gcaACGGCATGGCCGTGGGCCTGATCGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTgg                                                                                                                                   >  5:126563/1‑149 (MQ=255)                  caACGGCATGGCCGTGGGCCTGATCGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTggg                                                                                                                                  >  1:331621/1‑149 (MQ=255)                    aCGGCATGGCCGTGGGCCTGATCGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACTACGTAGGGCAGGATTGGGGc                                                                                                                                <  7:142663/149‑1 (MQ=255)                          tGGCCGTGGGCCTGATCGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTgc                                                                                                                       <  1:262184/149‑1 (MQ=255)                           ggCCGTGGGCCTGATCGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTgc                                                                                                                       <  4:250440/148‑1 (MQ=255)                            gCCGTGGGCCTGATCGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTgcgc                                                                                                                     <  6:126563/149‑1 (MQ=255)                                  ggCCTGATCGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTtctc                                                                                                               >  5:312052/1‑149 (MQ=255)                                   gCCTGATCGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCg                                                                                                              <  6:312052/149‑1 (MQ=255)                                   gCCTGATCGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCg                                                                                                              <  2:95607/149‑1 (MQ=255)                                          cGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGt                                                                                                                                                <  8:203343/111‑1 (MQ=255)                                          cGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGt                                                                                                                                                >  7:203343/1‑111 (MQ=255)                                            ccGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGc                                                                                                     <  2:331621/149‑1 (MQ=255)                                            ccGACCTGGCGGTGGGCGCCGATGATGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCCACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGc                                                                                                     >  7:213633/1‑149 (MQ=255)                                               aCCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTgc                                                                                                                       >  3:284386/1‑128 (MQ=255)                                               aCCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTgc                                                                                                                       <  4:284386/128‑1 (MQ=255)                                                  tGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGCCGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAg                                                                                               >  2:301056/1‑149 (MQ=255)                                                      ggTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCGACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCa                                                                                           <  4:127659/149‑1 (MQ=255)                                                      ggTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCGACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCa                                                                                           <  4:132687/149‑1 (MQ=255)                                                          ggCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATgg                                                                                       <  1:341834/149‑1 (MQ=255)                                                             gcCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGc                                                                                                     >  3:313171/1‑132 (MQ=255)                                                             gcCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGc                                                                                                     <  4:313171/132‑1 (MQ=255)                                                               cGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATc                                                                                  <  8:9098/149‑1 (MQ=255)                                                                    gTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGcc                                                                             <  5:338209/149‑1 (MQ=255)                                                                        cGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGCCCGGCCTGAGCGAGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGCGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCtac                                                                          >  3:249414/1‑146 (MQ=255)                                                                          gCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATc                                                                       <  1:161432/149‑1 (MQ=255)                                                                          gCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATc                                                                       >  4:222421/1‑149 (MQ=255)                                                                               cAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGAt                                                                  >  6:277319/1‑149 (MQ=255)                                                                                  gCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGCCTATCGCGCCTTCATGGAGAtgct                                                               >  1:47464/1‑149 (MQ=255)                                                                                      ggAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACAACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTCCgg                                                           >  2:285850/1‑149 (MQ=255)                                                                                       gAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑t‑‑tt                                                                                                                         <  5:253019/86‑1 (MQ=255)                                                                                       gAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑t‑‑tt                                                                                                                         >  6:253019/1‑86 (MQ=255)                                                                                       gAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCt                                                                                     <  8:137852/122‑1 (MQ=255)                                                                                       gAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCt                                                                                     >  7:137852/1‑122 (MQ=255)                                                                                       gAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCggg                                                          <  7:261540/149‑1 (MQ=255)                                                                                             ccAGGACGAACTGCTTGCCTGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGTGCA‑T‑‑TTGCGCCTTCTCGCCCGAAGGCCTCAAGCGCAATGGCTATCTCGCCTTCATCGAGATGCTGCGGGCCAAcc                                                    <  1:285850/149‑1 (MQ=255)                                                                                                  aCGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGcc                                               >  1:145555/1‑149 (MQ=255)                                                                                                    gAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGCCCa                                             >  7:71114/1‑149 (MQ=255)                                                                                                    gAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGAAGGGCAGGATTGGGGCAGTAGTTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGc                                                >  1:279729/1‑149 (MQ=255)                                                                                                          ctgGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGCCCACGCTgg                                       >  5:45337/1‑149 (MQ=255)                                                                                                              acgtccTGAGAGTGGGTGCGCCGCCCGACTTCCTCAAACTGGTAGGTCAGGATTGGGGGA‑T‑‑TTGCGCCTTGTGGAGCGTCGGCGTCAAGCGCAATGGCTATCGCGCGTTCATCGAGATGCTGCGGGCGAAGCTAGCCCACGCTggtggc                                   <  2:47464/145‑1 (MQ=255)                                                                                                               cggccTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGCCCACGCTGGCGGcc                                  >  8:208285/1‑149 (MQ=255)                                                                                                                ggccTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGCCCACGCTGGCGGCCt                                 <  2:145555/149‑1 (MQ=255)                                                                                                                  ccTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGCCCACGCTGGCGGCCTgc                               <  8:9942/149‑1 (MQ=255)                                                                                                                        cGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGCCCACGCTGGCGGCCTGCGCATc                          <  8:213633/148‑1 (MQ=255)                                                                                                                        cGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGCCCACGCTGGCGGCCTGCGCATCg                         <  3:222421/149‑1 (MQ=255)                                                                                                                          tGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGCCCACGCTGGCGGCCTGCGCATCGAc                       <  6:45337/149‑1 (MQ=255)                                                                                                                            ggTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGc                                                <  7:66399/122‑1 (MQ=255)                                                                                                                            ggTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGc                                                >  8:66399/1‑122 (MQ=255)                                                                                                                                  ccgccCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGCCCACGCTGGCGGCCTGCGCATCGACCATGTGAt               <  8:82229/149‑1 (MQ=255)                                                                                                                                    gccCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGCCCACGCTGGCGGCCTGCGCATCGACCATGTGATgg             <  5:277319/149‑1 (MQ=255)                                                                                                                                           aTCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGCCCACGCTGGCGGCCTGCGCATCGACCATGTGATGGGCCTgcg      >  5:147501/1‑149 (MQ=255)                                                                                                                                             ccTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGCCCACGCTGGCGGCCTGCGCATCGACCATGTGATGGGCCTGCGcc    >  7:94848/1‑149 (MQ=255)                                                                                                                                              cTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGCCCACGCTGGCGGCCTGCGCATCGACCATGTGATGGGCCTGCGCCg   >  6:96173/1‑149 (MQ=255)                                                                                                                                               tCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGCCCACGCTGGCGGCCTGCGCATCGACCATGTGATGGGCCTGCGCCgg  >  3:125082/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                      CACGAAGCGGCCCGCAGCAACGGCATGGCCGTGGGCCTGATCGCCGACCTGGCGGTGGGCGCCGATGGTGCCGGCAGCCAGGCCTGGAGCCGCCAGGACGAACTGCTGGCCGGCCTGAGCGTGGGTGCGCCGCCCGACATCCTCAACCACGTAGGGCAGGATTGGGGCA‑T‑‑TTGCGCCTTCTCGCCCGACGGCCTCAAGCGCAATGGCTATCGCGCCTTCATCGAGATGCTGCGGGCCAACCTCGCCCACGCTGGCGGCCTGCGCATCGACCATGTGATGGGCCTGCGCCGG  >  NC_002947/4569160‑4569450

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