Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,243,742	T→A	12.9%	intergenic (+31/+68)	tsaA → / ← rnt	peroxiredoxin/ribonuclease T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,243,742	0	T	A	12.9%	153.9 / 14.2	62	intergenic (+31/+68)	tsaA/rnt	peroxiredoxin/ribonuclease T
Reads supporting (aligned to +/- strand):  ref base T (28/26);  new base A (4/4);  total (32/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 = 6.79e-01

TTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCGGT  >  NC_002947/1243595‑1243890                                                                                                                                                    |                                                                                                                                                      ttCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTg                                                                                                                                                      >  5:169007/1‑149 (MQ=255) ttCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTg                                                                                                                                                      <  2:177222/149‑1 (MQ=255)            gCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGtttttt                                                                                                                                           <  2:220792/149‑1 (MQ=255)            gCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCACGAAGGCCGGtttttt                                                                                                                                           <  5:17330/149‑1 (MQ=255)             cACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGttttttt                                                                                                                                          <  6:169007/149‑1 (MQ=255)             cACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCACGAAGGCCGGttttt                                                                                                                                            <  7:63668/147‑1 (MQ=255)             cACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCACGAAGGCCGGttttt                                                                                                                                            >  8:63668/1‑147 (MQ=255)                 gCGAAGTGTGCCCTGCCAACTGG‑AAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGcc                                                                                                                                     >  5:204233/1‑149 (MQ=255)                     aGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTc                                                                                                                                  <  6:58785/149‑1 (MQ=255)                      gtgtGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGATTTTTTATGCCtaca                                                                                                                                 >  4:179081/1‑146 (MQ=255)                              tGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCgg                                                                                                                         >  6:295312/1‑149 (MQ=255)                                ccAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGttttttt                                                                                                                                          <  2:38959/130‑1 (MQ=255)                                ccAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGttttttt                                                                                                                                          >  1:38959/1‑130 (MQ=255)                                ccAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCgg                                                                                                                         <  6:129503/147‑1 (MQ=255)                                  aaCTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑ctc                                                                                                                    >  3:160892/1‑149 (MQ=255)                                   aCTGG‑AAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCtt                                                                                                                  <  6:204233/149‑1 (MQ=255)                                     tGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTc                                                                                                                 <  3:17584/149‑1 (MQ=255)                                       gAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt                                                                                                                                        <  7:51981/125‑1 (MQ=255)                                       gAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt                                                                                                                                        >  8:51981/1‑125 (MQ=255)                                             aGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCACGAAGGCCGGttttttc                                                                                                                                          >  7:40925/1‑116 (MQ=38)                                             aGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCACGAAGGCCGGttttttc                                                                                                                                          <  8:40925/117‑2 (MQ=38)                                                 gACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGccc                                                                                                     >  7:129819/1‑149 (MQ=255)                                                    aaGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCcgc                                                                                                   >  8:275130/1‑148 (MQ=255)                                                      gggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGaa                                                                                                >  4:136800/1‑149 (MQ=255)                                                        gCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATg                                                                                              >  1:281929/1‑149 (MQ=255)                                                         cATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGc                                                                                             >  5:199712/1‑149 (MQ=255)                                                              aCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCttt                                                                                         <  2:281929/148‑1 (MQ=255)                                                              aCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTa                                                                                        <  5:295312/149‑1 (MQ=255)                                                                 cTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCa                                                                                     <  4:160892/149‑1 (MQ=255)                                                                  ttCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCt                                                                                           >  8:65657/1‑142 (MQ=255)                                                                  ttCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCt                                                                                           <  7:65657/142‑1 (MQ=255)                                                                    cGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGa                                                                                   <  1:115691/148‑1 (MQ=255)                                                                            ggCTTTGCTGTCTACCTGAGCGAGAAGTCTGGCAAGTTGTATTTGCTGCGTGGTTTAAAAACCCGGCGTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑GTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGtt                                                                          <  3:179081/149‑1 (MQ=255)                                                                                         aCCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGACATCGCGCCAg                                                             >  6:137649/1‑149 (MQ=255)                                                                                          ccTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGTAATCGCGCCAgc                                                            >  1:266781/1‑149 (MQ=255)                                                                                          ccTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgc                                                            >  4:24172/1‑149 (MQ=255)                                                                                           cTGAGCGAGAACGCTGGCAATCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTa                                                                                                                                         <  1:82287/72‑1 (MQ=255)                                                                                           cTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTa                                                                                                                                         >  2:82287/1‑72 (MQ=255)                                                                                           cTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgcc                                                           >  8:163113/1‑149 (MQ=255)                                                                                            tGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTa                                                                                        <  8:225013/119‑1 (MQ=255)                                                                                            tGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTa                                                                                        >  7:225013/1‑119 (MQ=255)                                                                                              aGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgccgc                                                         <  2:293044/148‑1 (MQ=255)                                                                                               gCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgccgcc                                                        >  3:299269/1‑148 (MQ=255)                                                                                                           gCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAg                                           <  2:218078/149‑1 (MQ=255)                                                                                                            cAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGc                                          <  7:275130/149‑1 (MQ=255)                                                                                                             aaGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCg                                         >  4:251214/1‑149 (MQ=255)                                                                                                             aaGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCg                                         <  2:266781/149‑1 (MQ=255)                                                                                                             aaGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCg                                         >  3:84429/1‑149 (MQ=255)                                                                                                              aGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCgg                                        >  7:236637/1‑149 (MQ=255)                                                                                                               gCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGc                                          >  1:34770/1‑146 (MQ=255)                                                                                                               gCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGc                                          <  2:34770/146‑1 (MQ=255)                                                                                                               gCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGt                                       >  1:27315/1‑149 (MQ=255)                                                                                                                 tGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTg                                     >  1:120420/1‑149 (MQ=255)                                                                                                                             gTGCATTAAAAAACCGGCCTTCGTGGCGGGTTTTTTTATGCCTTCCTGTGCGATCTCTCTTCGAATGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCa                          <  1:260175/149‑1 (MQ=255)                                                                                                                                cATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa                       >  8:71697/1‑148 (MQ=38)                                                                                                                                 ataaaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttcttt                                               >  4:280584/4‑123 (MQ=38)                                                                                                                                 ataaaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttcttt                                               <  3:280584/120‑1 (MQ=38)                                                                                                                                    aaaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCATTAGCAGATCAATAGTTGAAATAGCGCCAGCCGCACATTTCATTCCAGCGGGTGACGATGCCGCAGAACAGCt                  >  2:162222/1‑149 (MQ=255)                                                                                                                                         aCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCg             <  6:199712/149‑1 (MQ=255)                                                                                                                                               cTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTtc       >  1:22104/1‑149 (MQ=255)                                                                                                                                                  cGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTtctc     <  2:120420/148‑1 (MQ=255)                                                                                                                                                  cGTGGCCGGTTTTTTTATGCCTTCCTGGGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTtctc     <  2:64889/148‑1 (MQ=255)                                                                                                                                                    tGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCGgt  <  4:84429/149‑1 (MQ=255)                                                                                                                                                    tGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCGgt  <  3:24172/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      TTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGC‑CTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCGGT  >  NC_002947/1243595‑1243890

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