Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,289,861:1	+A	8.9%	coding (1128/2349 nt)	PP_1126 →	hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,289,861	1	.	A	8.9%	204.2 / 16.1	56	coding (1128/2349 nt)	PP_1126	hydrolase
Reads supporting (aligned to +/- strand):  ref base . (31/20);  new base A (3/2);  total (34/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-01

TCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACCGCTGGCA  >  NC_002947/1289721‑1290010                                                                                                                                              |                                                                                                                                                      tCCCTGGGCGTGAACGGCGTGGCAGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAACGCATGTTCATGGCCTTGCCCGCTGAGGCTGATGCACTAGCCGCTTTCTACGGAGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCaa                                                                                                                                              <  2:272702/150‑1 (MQ=255)   ccTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGA‑GGCATTGCCAtcgtcgttgtt                                                                                                                                            <  2:1589663/150‑6 (MQ=255)   ccTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGA‑GGCATTGCCAtcgtc                                                                                                                                                  >  2:1117206/1‑144 (MQ=255)   ccTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGA‑GGCATTGCCAtcgtc                                                                                                                                                  <  1:1117206/144‑1 (MQ=255)             aaCGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑tcgtc                                                                                                                                                  >  2:1540599/1‑134 (MQ=255)             aaCGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑tcgtc                                                                                                                                                  <  1:1540599/134‑1 (MQ=255)               cGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCaa                                                                                                                                <  1:675011/150‑1 (MQ=255)                gCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAg                                                                                                                               >  1:348125/1‑150 (MQ=255)                    tGGCGGCAGACGGCGTGTGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAgg                                                                                                                           <  1:491439/150‑1 (MQ=255)                      gcggcAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAgggg                                                                                                                         <  1:46221/150‑1 (MQ=255)                       cggcAGAGGGCGGGCCTACCTATGTCGCCGGGCGCGAAGGCATGTTGCAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAACGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGgc                                                                                                                        <  2:600286/150‑1 (MQ=255)                       cggcAGACGGCGGTCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGgc                                                                                                                        >  1:1215311/1‑150 (MQ=255)                                gcgggcgTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAg                                                                                                                            >  1:630002/1‑137 (MQ=255)                                gcgggcgTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAg                                                                                                                            <  2:630002/137‑1 (MQ=255)                                  gggcgTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGgcgc                                                                                                                      <  1:888779/141‑1 (MQ=255)                                  gggcgTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGgcgc                                                                                                                      >  2:888779/1‑141 (MQ=255)                                           tATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTAcgcc                                                                                                    <  2:1215311/150‑1 (MQ=255)                                                          gAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGt                                                                                     >  2:1590887/1‑150 (MQ=255)                                                                aTGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTAt                                                                               >  2:1145791/1‑150 (MQ=255)                                                                    tCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCACTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCa                                                                           >  2:1186715/1‑150 (MQ=255)                                                                      aaGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCaaa                                                                         >  2:1431663/1‑150 (MQ=255)                                                                       aGGTCTTGCCCGCTGAGGGCGATGCACTGGCCGATATGCCCGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAAc                                                                        <  1:499004/150‑1 (MQ=255)                                                                        ggCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAAcc                                                                       <  1:1145791/150‑1 (MQ=255)                                                                         gCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAAccc                                                                      >  1:1522507/1‑150 (MQ=255)                                                                           cTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTAAGCCGCCAACCTGCAGCGTACCTATGGCCAATCCTg                                                                    >  2:92850/1‑150 (MQ=255)                                                                                       gggCGATGCACTGACCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGt                                                        >  2:1573174/1‑150 (MQ=255)                                                                                          cGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCCGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGa                                                     >  1:192926/1‑150 (MQ=255)                                                                                           gATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAg                                                    >  1:686570/1‑150 (MQ=255)                                                                                            aTGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCATGCGCGGCTGCTGGTTCAACTCGGACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAgg                                                   <  1:1019734/150‑1 (MQ=255)                                                                                                  tGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGgccagc                                             >  1:173817/1‑150 (MQ=255)                                                                                                  tGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGgccagc                                             >  1:173956/1‑150 (MQ=255)                                                                                                      cGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCcagccagc                                         >  1:1398811/1‑150 (MQ=255)                                                                                                      cGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCcagccagc                                         <  1:1431663/150‑1 (MQ=255)                                                                                                      cGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCcagccagc                                         >  2:1039001/1‑150 (MQ=255)                                                                                                         tttCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGa                                      >  2:499919/1‑150 (MQ=255)                                                                                                               tGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTgcgc                                >  2:202058/1‑150 (MQ=255)                                                                                                               tGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGACCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTgcgc                                >  1:377503/1‑150 (MQ=255)                                                                                                                   gAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTc                                                                                                                >  2:980831/1‑66 (MQ=255)                                                                                                                   gAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTc                                                                                                                <  1:980831/66‑1 (MQ=255)                                                                                                                   gAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGcgccg                            <  1:202058/150‑1 (MQ=255)                                                                                                                    aaaGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGcgccgc                           >  1:1048133/1‑150 (MQ=255)                                                                                                                     aaGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCgccgcc                          >  1:787313/1‑150 (MQ=255)                                                                                                                       gCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGccgccgc                        >  1:865351/1‑149 (MQ=255)                                                                                                                         ccttgcgacga‑GGCATTGCCATCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTg                      >  1:183329/6‑150 (MQ=255)                                                                                                                         aacaacGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTg                      >  2:1211681/1‑150 (MQ=255)                                                                                                                           ttgcgacga‑GGCATTGCCATCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCt                    >  2:1492516/4‑150 (MQ=255)                                                                                                                           caacGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCt                    >  2:17854/1‑150 (MQ=255)                                                                                                                            aacGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGc                                                                                          >  2:1322352/1‑79 (MQ=255)                                                                                                                            aacGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGc                                                                                          <  1:1322352/79‑1 (MQ=255)                                                                                                                            aacGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTg                   <  2:1398811/150‑1 (MQ=255)                                                                                                                             acgacgATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTgg                  >  2:1266019/1‑150 (MQ=255)                                                                                                                               gacgaTGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGa                                                     >  1:596602/1‑113 (MQ=255)                                                                                                                               gacgaTGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGa                                                     <  2:596602/113‑1 (MQ=255)                                                                                                                                  gaTGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAAACCTACCCGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGcc             >  1:1601484/1‑150 (MQ=255)                                                                                                                                    tGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGa           <  1:17854/150‑1 (MQ=255)                                                                                                                                       cAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGAcc         <  2:1055072/149‑1 (MQ=255)                                                                                                                                       cAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGAcc         >  1:1055072/1‑149 (MQ=255)                                                                                                                                          tGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACCGCTg     <  2:787313/150‑1 (MQ=255)                                                                                                                                             c‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACCGc       <  2:1525853/145‑1 (MQ=255)                                                                                                                                             c‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACCGc       >  1:1525853/1‑145 (MQ=255)                                                                                                                                             c‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACCGCTGGCa  >  2:450180/1‑150 (MQ=255)                                                                                                                                              |                                                                                                                                                      TCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCC‑TCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACCGCTGGCA  >  NC_002947/1289721‑1290010

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