Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	817,870	C→T	7.2%	A100T (GCG→ACG)	PP_0703 ←	transmembrane sensor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	817,870	0	C	T	7.2%	160.7 / 2.4	55	A100T (GCG→ACG)	PP_0703	transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base C (25/26);  new base T (2/2);  total (27/28)
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.94e-01

CGTTGCTCCGCGCTGTAGTGCAAGTCGATGGCACTGGCGGTGTTGAGCTGCAGACGGCTGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGTGCGCCGGCTCAGGTCCTTGGGCAGCCCGGCCAATTCATCACCCAGATTACTGA  >  NC_002947/817725‑818016                                                                                                                                                  |                                                                                                                                                   cGTTGCTCCGCGCTGTAGTGCAAGTCGATGGCACTGGCGTTGTTGAGCTGCAGACGGCTGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGgc                                                                                                                                                 <  6:376040/149‑1 (MQ=255)         cgcgCTGTAGTGCAAGTCGATGGCACTGGCGGTGTTGAGCTGCAGACGGCTGCCGTCACTTCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGTGgcccagcc                                                                                                                                           >  1:67549/1‑147 (MQ=255)         cgcgCTGTAGTGCAAGTCGATGGCACTGGCGGTGTTGAGCTGCAGACGGCTGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGTGgcccagcc                                                                                                                                           >  1:72885/1‑147 (MQ=255)         cgcgCTGTAGTGCAAGTCGATGGCACTGGCGGTGTTGAGCTGCAGACGGCTGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGAGGCAGGCCTAGCGGTTCGCGCGCGgcccagcc                                                                                                                                           >  7:100271/1‑147 (MQ=255)                       aGTCGATGGCACTGGCGGTGTTGAGCTGCAGACGGCTGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGcccagcccag                                                                                                                                        <  6:243558/136‑1 (MQ=255)                       aGTCGATGGCACTGGCGGTGTTGAGCTGCAGACGGCTGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGcccagcccag                                                                                                                                        >  5:243558/1‑136 (MQ=255)                           gATGGCACTGGCGGTGTTGAGCTGCAGACGGCTGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGcc                                                                                                                        >  3:211501/1‑148 (MQ=255)                             tGGCACTGGCGGTGTTGAGCTGCAGACGGCTGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCaa                                                                                                                      >  5:18459/1‑148 (MQ=255)                               gCACTGGCGGTGTTCAGCTGCAGACGGATGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACa                                                                                                                    <  1:244808/148‑1 (MQ=255)                                cACTGGCGGTGTTGAGCTGCAGACGGCTGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGc                                                                                                                  <  5:241871/149‑1 (MQ=255)                                     gCGGTGTTGAGCTGCAGACGTCTGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCa                                                                                                              <  7:324272/148‑1 (MQ=255)                                     gCGGTGTTGAGCTGCAGACGGCTGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCa                                                                                                              <  8:8934/148‑1 (MQ=255)                                     gCGGTGTTGAGCTGCAGACGGCTGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCa                                                                                                              <  1:238507/148‑1 (MQ=255)                                       ggTGTTGAGCTGCAGACGGCTGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCa                                                                                                           <  8:163141/149‑1 (MQ=255)                                        gtgtTGAGCTGCAGACGGCTGCCGTCACTTCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGTGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCaa                                                                                                          <  2:67549/149‑1 (MQ=255)                                        gtgtTGAGCTGCAGACGGCTGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGTGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCaa                                                                                                          <  2:72885/149‑1 (MQ=255)                                         tgtTGAGCTGCAGACGGCTGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAAc                                                                                                         >  4:305404/1‑149 (MQ=255)                                                   cAGACGGCTGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAggtggt                                                                                               >  8:34182/1‑149 (MQ=255)                                                        ggCTGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTg                                                                                          <  7:34182/149‑1 (MQ=255)                                                              cGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCggctg                                                                                    >  4:51706/1‑149 (MQ=255)                                                               gTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTcc                                                                                             >  1:267496/1‑139 (MQ=255)                                                               gTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTcc                                                                                             <  2:267496/139‑1 (MQ=255)                                                                  aCTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTc                                                                                              <  4:346110/135‑1 (MQ=255)                                                                  aCTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTc                                                                                              >  3:346110/1‑135 (MQ=255)                                                                  aCTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCtt                                                                                >  8:4872/1‑149 (MQ=255)                                                                      ccTTGCCATTGGCGCCGTTCGCGTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCtttt                                                                              <  1:240341/147‑1 (MQ=255)                                                                      ccTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCtttt                                                                              <  7:4872/147‑1 (MQ=255)                                                                      ccTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTgg                                                                            <  5:143750/149‑1 (MQ=255)                                                                      ccTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTg                                                                             <  4:211501/148‑1 (MQ=255)                                                                      ccTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTg                                                                             <  6:221323/148‑1 (MQ=255)                                                                      ccTTGACATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTgg                                                                            >  5:221323/1‑149 (MQ=255)                                                                        ttGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGt                                                                          >  4:341714/1‑149 (MQ=255)                                                                        ttGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGt                                                                          >  7:347447/1‑149 (MQ=255)                                                                             aTTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTgg                                                                     >  4:326781/1‑149 (MQ=255)                                                                             aTTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTgg                                                                     >  6:169912/1‑149 (MQ=255)                                                                                   gcCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCtttt                                                                              >  1:237731/1‑134 (MQ=255)                                                                                   gcCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCtttt                                                                              <  2:237731/134‑1 (MQ=255)                                                                                           ccTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGt                                                       <  8:347447/149‑1 (MQ=255)                                                                                            cTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGTg                                                      >  2:4989/1‑149 (MQ=255)                                                                                                 gCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACgg                                                                                                                          >  7:164752/1‑76 (MQ=255)                                                                                                 gCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACgg                                                                                                                          <  8:164752/76‑1 (MQ=255)                                                                                                 gCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGTGCGCCg                                                 <  6:18459/149‑1 (MQ=255)                                                                                                          gtgtCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCa                                                           >  7:162339/1‑130 (MQ=255)                                                                                                          gtgtCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCa                                                           <  8:162339/130‑1 (MQ=255)                                                                                                            gtCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGTGCGCCGGCTCAGGTCCt                                      >  3:170521/1‑149 (MQ=255)                                                                                                               ggCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGTGCGCCGGCTCAGGTCCTTgg                                   >  8:334951/1‑149 (MQ=255)                                                                                                                gCCAGCAGCTGTGGCAGGCCTAGCGGTTAGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCAAACAGCACCAGCAACTTCAGGTGGTCCCTGAGGCTGGCTATTGGGTGCTGGCTAGCGACCAGGGTGCGACGGCTCAGGTCCTTggg                                  >  6:230624/1‑149 (MQ=255)                                                                                                                         tgtgGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGTGCGCCGGCTCAGGTCCTTGGGCAGCCCGGc                         >  7:149438/1‑149 (MQ=255)                                                                                                                           tgGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCATCTTCAGGTGGTCCCCGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGTTGCGCCGGCTCAGGTCCTTGGGCAGCCCGGcc                        <  6:295760/148‑1 (MQ=255)                                                                                                                           tgGCAGGCCTAGCGGGTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGGCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGTGCGCCGGCTCAGGTCCTTGGGCAGCCCGGcc                        <  8:249668/148‑1 (MQ=255)                                                                                                                             gCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGTGCGCCGGCTCAGGTCCTTGGGCAGCCCGGCCaa                      <  3:341714/148‑1 (MQ=255)                                                                                                                             gCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGTGCGCCGGCTCAGGTCCTTGGGCAGCCCGGCCAAt                     >  2:210420/1‑149 (MQ=255)                                                                                                                                      gCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGTGCGCCGGCTCAGGTCCTTGGGCAGCCCGGCCAATTCATCAcc             >  2:69514/1‑148 (MQ=255)                                                                                                                                      gCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGTGCGCCGGCTCAGGTCCTTGGGCAGCCCGGCCAATTCATCAc              <  3:326781/147‑1 (MQ=255)                                                                                                                                          ttCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGTGCGCCGGCTCAGGTCCTTGGGCAGCCCGGCCAATTCATCACCCAGAt        <  3:51706/149‑1 (MQ=255)                                                                                                                                              cgcgcgGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGTGCGCCGGCTCAGGTCCTTGGGCAGCCCGGCCAATTCATCACCCAGATTACt    >  7:217671/1‑149 (MQ=255)                                                                                                                                               gcgcgGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGTGCGCCGGCTCAGGTCCtt                                     <  8:29633/115‑1 (MQ=255)                                                                                                                                               gcgcgGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGTGCGCCGGCTCAGGTCCtt                                     >  7:29633/1‑115 (MQ=255)                                                                                                                                                cgcgGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGTGCGCCGGCTCAGGTCCTTGGGCAGCCCGGCCAATTCATCACCCAGATTACTGa  >  3:300612/1‑149 (MQ=255)                                                                                                                                                cgcgGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGTGCGCCGGCTCAGGTCCTTGGGCAGCCCGGCCAATTCATCACCCAGATTACTGa  >  8:344000/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CGTTGCTCCGCGCTGTAGTGCAAGTCGATGGCACTGGCGGTGTTGAGCTGCAGACGGCTGCCGTCACTGCCTTGCCATTGGCGCCGTTCGCCTGTGGCAGTGCTGGTGTCGGCCAGCAGCTGTGGCAGGCCTAGCGGTTCGCGCGCGGCCCAGCCCAGGCTGCCGCCCACGGCCAACAGCCCCAGCAACTTCAGGTGGTCCCTGCGGCTGGCTTTTGGGTGCTGGCTAGCGTCCAGGGTGCGCCGGCTCAGGTCCTTGGGCAGCCCGGCCAATTCATCACCCAGATTACTGA  >  NC_002947/817725‑818016

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