Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,921,550	T→A	11.4%	intergenic (+25/+36)	PP_1720 → / ← PP_1721	alcohol dehydrogenase/phosphoserine phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,921,550	0	T	A	11.4%	163.9 / 6.5	68	intergenic (+25/+36)	PP_1720/PP_1721	alcohol dehydrogenase/phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base T (32/28);  new base A (3/5);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAA  >  NC_002947/1921405‑1921690                                                                                                                                                  |                                                                                                                                              gTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTAcc                                                                                                                                            >  3:212686/1‑149 (MQ=255)       gccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAACCCGCGaa                                                                                                                                      <  1:363204/149‑1 (MQ=255)       gccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCCGGCCCCTTCGCGGGTAAACCcgcg                                                                                                                                        >  3:20175/1‑147 (MQ=255)       gccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCCGGCCCCTTCGCGGGTAAACCCGCGaa                                                                                                                                      <  4:20175/149‑1 (MQ=255)          gcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAggg                                                                                                                                   >  2:281575/1‑149 (MQ=255)          gcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAACCCGCGAAggg                                                                                                                                   <  3:263103/149‑1 (MQ=255)               tGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑gg                                                                                                                             <  7:192540/149‑1 (MQ=255)                           aTCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCCAAGGGGCC‑GGCACAGCCCATa                                                                                                                  >  8:127409/1‑148 (MQ=255)                              cGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATaaaaa                                                                                                              <  5:144449/149‑1 (MQ=255)                               ggTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAgg                                                                                                                                    <  5:212165/127‑1 (MQ=255)                               ggTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAgg                                                                                                                                    >  6:212165/1‑127 (MQ=255)                                gTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTAGCGGGTTTACCCGCGAAGGGGCC‑GGCCCAGCCCATAAAAAcc                                                                                                            >  1:182577/1‑149 (MQ=255)                                  ggAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCtt                                                                                                          <  4:210441/149‑1 (MQ=255)                                   gtttttGCCCTGGGGGGGTATGCCGAGGTGGATACCAACCTGGTGTGGGGAAAGTTGGTGATGGTTATCTACCTTAGCCTGGTTTATGGGGGGCTGGCCCTTTCGCGGTTTTACCCGCTAAGGGGCC‑GGCACAGCCCATAAAAACCTTc                                                                                                         <  6:7992/146‑1 (MQ=255)                                   gAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTTATGGACCCTCGCCTGGTGTACGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCTAAGGGGCC‑GGCACAGCCCATAAAAACCTTc                                                                                                         <  4:100548/149‑1 (MQ=255)                                               gcggcgTATGCCGAGCTGGAGACCAACGAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCa                                                                                             <  4:71330/149‑1 (MQ=255)                                               gcggcgTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCa                                                                                             <  3:248257/149‑1 (MQ=255)                                               gcggcgTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCa                                                                                             <  5:339246/149‑1 (MQ=255)                                                       tGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCgg                                                                                      <  1:43925/148‑1 (MQ=255)                                                                 gagaCAAACCTGGTGGCGGGCAAGCTGGTTATGGTGATCGCCCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCAtt                                                                             <  8:154093/147‑1 (MQ=255)                                                                      caaccaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTcagctccagc                                                                                                <  5:331756/123‑1 (MQ=255)                                                                      caaccaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTcagctccagc                                                                                                >  6:331756/1‑123 (MQ=255)                                                                      caaccaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATATCg                                                                      >  3:38927/1‑149 (MQ=255)                                                                       aaccaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCgg                                                                     >  1:208876/1‑149 (MQ=255)                                                                           aGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCAATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATg                                                                 >  2:64105/1‑149 (MQ=255)                                                                               gtCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTcg                                                             >  6:182706/1‑149 (MQ=255)                                                                                 cGGGCAAGGGGGGGATGGGGAGCGCCCCTCGCCGGGTGTACGCGGGGCGGCCCCCTTCGCGGGTTACCCCGCTAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTcgcg                                                           <  2:210808/149‑1 (MQ=255)                                                                                   ggCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCgcgc                                                          >  3:293484/1‑148 (MQ=255)                                                                                    gCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCaa                                                        <  7:96529/149‑1 (MQ=255)                                                                                      aaGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑CGCACAGCCCATAAAAACCTTCACCACCACCTCAGAGTGGGCCCCCCAATTAAATCGGCATTCTCCCGCAaca                                                      >  2:198645/1‑149 (MQ=255)                                                                                         cTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTg                                                   >  1:50559/1‑149 (MQ=255)                                                                                          tggtgatggtgaTCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGc                                                  >  5:115594/1‑149 (MQ=255)                                                                                          tggtgatggtgaTCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCAATTCTCGGGTTTACAAGCGAATGGGCC‑GGCACAGACCATAAAAACCTTCAGCTCCAGGTCAGAGTCGGCCCACCATTTATCTAGGCATTCTCGCACAACACTGc                                                  >  5:349666/1‑149 (MQ=255)                                                                                            gtgatggtgatCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTACAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGt                                                >  4:271966/1‑149 (MQ=255)                                                                                                tggtgatCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCgg                                            <  5:182706/149‑1 (MQ=255)                                                                                                 ggtgatCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCggg                                           >  7:354568/1‑149 (MQ=255)                                                                                                  gtgatCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTAGCGGGTTTAACCGCGAAGGGGCC‑GGCACAGCCCATAAAAAACTTAAGCAACAGCTCAGAGTAGGCCACCCATTTATATCGGCAAGCTCGCGCAACACATCGTCCCGGt                                          >  4:366241/1‑149 (MQ=255)                                                                                                       cGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTcac                                      >  1:109697/1‑148 (MQ=255)                                                                                                        gACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAACAACCTTCAGCTCCCGCTCAGAGTCGGCCACCCATTGATCTCGGCATCCTCGCCCCACACTGCGTCCGGGTTcacca                                    >  8:96609/1‑149 (MQ=255)                                                                                                         aCCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTcacca                                    >  7:5209/1‑148 (MQ=255)                                                                                                           ccTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCCTTGATCTCGGCATTCTCGCGCAACACTGCGTCCGGGTTCACCACAt                                 >  7:89656/1‑149 (MQ=255)                                                                                                           ccTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACAt                                 <  1:281575/149‑1 (MQ=255)                                                                                                           ccTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTAAGCTCCAGCTCAGCGTCGGCCACCCATTGATCTCGGCATGCTCGCGCACCCCTGCGTCCGGGTTCACCCCAt                                 >  4:73961/1‑149 (MQ=255)                                                                                                            cTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATg                                <  6:115594/149‑1 (MQ=255)                                                                                                            cTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATg                                >  6:280864/1‑149 (MQ=255)                                                                                                             tAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGc                               >  6:1487/1‑149 (MQ=255)                                                                                                               gCCGGGTGAACGCGGGGGGGCCCCCTTCGCGGGTTTCCCCGCTACGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGc                               <  2:182577/147‑1 (MQ=255)                                                                                                                ccTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTca                                       >  1:256307/1‑138 (MQ=255)                                                                                                                ccTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTca                                       <  2:256307/138‑1 (MQ=255)                                                                                                                ccTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCggg                            <  4:293484/149‑1 (MQ=255)                                                                                                                ccTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCggg                            <  2:208876/149‑1 (MQ=255)                                                                                                                ccTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCggg                            <  2:109697/149‑1 (MQ=255)                                                                                                                    gtgtAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGa                         <  5:1487/148‑1 (MQ=255)                                                                                                                        aaGCGGGGCTGCCCCCTTCGCGGGTTAACCCGCGATGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCa                     <  7:127409/148‑1 (MQ=255)                                                                                                                         aGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCAct                   <  1:1419/149‑1 (MQ=255)                                                                                                                         aGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCAct                   >  2:1419/1‑149 (MQ=255)                                                                                                                         aGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCAc                    >  1:267489/1‑148 (MQ=255)                                                                                                                         aGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCAc                    <  2:267489/148‑1 (MQ=255)                                                                                                                          gCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCActc                  <  5:280864/149‑1 (MQ=255)                                                                                                                              ggCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCAATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTcaa              <  1:64105/149‑1 (MQ=255)                                                                                                                                  ggCCCCTTCGCGGGTAAACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTcaacaa           >  1:233776/1‑148 (MQ=255)                                                                                                                                   gCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGtt                                         >  1:358376/1‑117 (MQ=255)                                                                                                                                     cccTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCaacaacaac       <  1:49255/149‑1 (MQ=255)                                                                                                                                      ccTTCGCGGGTAAACCCGCGAAGGGGCCGGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCaacaacaac       >  1:27977/1‑149 (MQ=255)                                                                                                                                       cTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACgg     >  2:7602/1‑149 (MQ=255)                                                                                                                                       cTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTAAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACAACATGCGGGTGCTCCACTCTCACCAACAACgg     >  2:150975/1‑149 (MQ=255)                                                                                                                                        ttCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACgg     >  6:244833/1‑148 (MQ=255)                                                                                                                                           gcgGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCaa  <  8:5209/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                              GTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCC‑GGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAA  >  NC_002947/1921405‑1921690

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