Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,081,115	A→C	7.7%	intergenic (+33/+82)	efp → / ← ohr	elongation factor P/organic hydroperoxide resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,081,115	0	A	C	7.7%	180.9 / 5.3	64	intergenic (+33/+82)	efp/ohr	elongation factor P/organic hydroperoxide resistance protein
Reads supporting (aligned to +/- strand):  ref base A (28/30);  new base C (2/3);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

CCAAGCTGGCGAACGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGAGTACGGGCACACCTG  >  NC_002947/2080973‑2081262                                                                                                                                                |                                                                                                                                                     ccAAGCTGGCGAACGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCggg                                                                                                                                                >  6:65261/1‑149 (MQ=255) ccAAGCTGGCGAACGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCg                                                                                                                                                  <  8:155190/147‑1 (MQ=255)        ggCGAACGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTtg                                                                                                                                         >  6:195245/1‑149 (MQ=255)         gCGAACGGTACGGAGATCTCGGTTGCCGACTTCATCCAGAGCGTCGTGTGGATCGACATCGTCACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGTCGCAAAAAA‑CCCGGCGCAAGCCGGGtttttt                                                                                                                                          <  4:293012/147‑1 (MQ=255)         gCGAACGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCTTGCGCCGGGtttttt                                                                                                                                          <  5:187619/147‑1 (MQ=255)         gCGAACGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCACGCCGGGtttttt                                                                                                                                          <  5:91871/147‑1 (MQ=255)           gAACGGTACCGAGAGCTCGGTTGCCGACTTCATCCAGATCGACGACTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGc                                                                                                                                      <  5:127013/149‑1 (MQ=255)            aaCGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTTCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGcc                                                                                                                                     >  6:39212/1‑149 (MQ=255)                    cGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCgg                                                                                                                                                 >  8:53371/1‑129 (MQ=255)                    cGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCgg                                                                                                                                                 <  7:53371/129‑1 (MQ=255)                    cGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCGTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGtttt                                                                                                                             >  3:170125/1‑149 (MQ=255)                     gagaTCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTc                                                                                                                            >  6:22878/1‑149 (MQ=255)                      agaTCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTcc                                                                                                                           <  2:40271/149‑1 (MQ=255)                       gaTCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTcc                                                                                                                           <  1:225366/149‑1 (MQ=255)                       gaTCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTcc                                                                                                                           >  3:264028/1‑148 (MQ=255)                          ctcGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCaa                                                                                                                                                      >  8:4244/1‑118 (MQ=255)                          ctcGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAAt                                                                                                                       <  8:159833/149‑1 (MQ=255)                             ggTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAAt                                                                                                                       <  4:98357/146‑1 (MQ=255)                                     aCTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTTCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGt                                                                                                            <  5:39212/149‑1 (MQ=255)                                     aCTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGt                                                                                                            <  7:77819/149‑1 (MQ=255)                                        tCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑Ct                                                                                                        >  4:189041/1‑149 (MQ=255)                                          aTCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGGTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑Ctt                                                                                                       <  3:304117/148‑1 (MQ=255)                                           tCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTtgtg                                                                                                                                       >  2:151877/1‑116 (MQ=255)                                               gATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTc                                                                                                 >  3:104492/1‑149 (MQ=255)                                                     cgaGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAg                                                                                           <  4:264028/149‑1 (MQ=255)                                                            atcgacatcgacACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATcg                                                                                    <  5:65261/149‑1 (MQ=255)                                                                 catcgacaCTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCaa                                                                               <  2:308774/149‑1 (MQ=255)                                                                 catcgacaCTCCCGACAACATCTTCAACCGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCaa                                                                               <  8:291344/149‑1 (MQ=255)                                                                 catcgacaATCCCGACAACAGCTTCAAAGGCCTCTCCATCAAGTAATGCTTCTTGCGAGACTCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGAGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CGTGTGCTCCCACAGGGGATCGCGCaa                                                                               <  3:46195/149‑1 (MQ=255)                                                                     gacaCTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCt                                                                           >  8:276672/1‑149 (MQ=255)                                                                       cacTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCttt                                                                         <  5:22878/149‑1 (MQ=255)                                                                        acTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTg                                                                        >  4:132645/1‑149 (MQ=255)                                                                          tCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTgcg                                                                      >  8:108544/1‑149 (MQ=255)                                                                           cgcgACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAAt                                                                                                                       <  4:302593/100‑1 (MQ=255)                                                                           cgcgACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAAt                                                                                                                       >  3:302593/1‑100 (MQ=255)                                                                             cgACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCt                                                                   <  7:108544/149‑1 (MQ=255)                                                                              gACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCtt                                                                  <  5:98037/149‑1 (MQ=255)                                                                                    aGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACgg                                                            <  7:276672/149‑1 (MQ=255)                                                                                      cTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGt                                                           >  4:296369/1‑148 (MQ=255)                                                                                      cTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCAGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCCTG‑CTTCTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGt                                                           >  2:264014/1‑148 (MQ=255)                                                                                            aGGCCGCTCCAAGTAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGc                                                    <  3:233826/149‑1 (MQ=255)                                                                                              gCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTAAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAg                                                  >  8:186719/1‑149 (MQ=255)                                                                                               ccGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAgg                                                 >  6:112849/1‑149 (MQ=255)                                                                                                cGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGgc                                                <  5:195245/149‑1 (MQ=255)                                                                                                  cTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGgcgc                                              <  3:102838/149‑1 (MQ=255)                                                                                                  cTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGgcg                                               >  7:96042/1‑148 (MQ=255)                                                                                                   tCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGca                                             >  6:212906/1‑149 (MQ=255)                                                                                                   tCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGca                                             >  6:63150/1‑149 (MQ=255)                                                                                                     cAAGAAGCATTACTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGcaca                                           >  3:109648/1‑149 (MQ=255)                                                                                                      aagaagTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACAt                                          >  5:282777/1‑149 (MQ=255)                                                                                                       agaagTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGgc                                                >  6:316733/1‑142 (MQ=255)                                                                                                       agaagTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGgc                                                <  5:316733/142‑1 (MQ=255)                                                                                                              aTGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGTCGCACATCCACAtt                                   >  8:185492/1‑148 (MQ=255)                                                                                                                gCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTAcc                                >  2:135330/1‑149 (MQ=255)                                                                                                                 cttcttGCGATACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTAcc                                <  3:189041/148‑1 (MQ=255)                                                                                                                        ggagaCGCAAAAAA‑CCCGGCGCATGCGGGGTTTTTTGGGCGTGGGTTTTCCCAATCTTCAGCGGGGG‑GTTGTGGTCCAACTGGTTCTCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGtgtg                         <  1:264014/147‑1 (MQ=255)                                                                                                                            gcacAAAAAA‑CCCGGCTTGCGCCGGGTTTTTTGCGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTAc                                                        >  4:126666/4‑113 (MQ=38)                                                                                                                            gcacAAAAAA‑CCCGGCTTGCGCCGGGTTTTTTGCGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTAc                                                        <  3:126666/110‑1 (MQ=38)                                                                                                                             cGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACATGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTg                   >  2:84281/1‑149 (MQ=255)                                                                                                                               cAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGa                 <  1:135330/149‑1 (MQ=255)                                                                                                                               cAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGa                 <  7:186719/149‑1 (MQ=255)                                                                                                                                aaaaaa‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTAACGTGCAGGCGCACATCCACATTACCGCGAGTGGAATTg                   >  8:116978/1‑146 (MQ=255)                                                                                                                                   aaa‑CCCGGCTTGCGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGAGTAc             >  4:182691/1‑149 (MQ=255)                                                                                                                                    aa‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGAGTACg            <  4:109648/149‑1 (MQ=255)                                                                                                                                       cccGGCGCAAGCCGGGTTTTTTGTGCCTGTGTCTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGCTTACGTGCAGGCGCACATCCACATTCCCGCGTGTGCCATACGACTACggg          <  4:170125/149‑1 (MQ=255)                                                                                                                                               aaGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGACTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGAGTACGGGCACACCt   >  1:233550/1‑149 (MQ=255)                                                                                                                                               aaGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGAGTACGGGCACACCTg  >  8:174237/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     CCAAGCTGGCGAACGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAA‑CCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTG‑CTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGAGTACGGGCACACCTG  >  NC_002947/2080973‑2081262

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