Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,810,284:1	+A	6.2%	intergenic (+45/‑124)	eno → / → ftsB	enolase/cell division protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,810,284	1	.	A	6.2%	255.5 / 11.7	64	intergenic (+45/‑124)	eno/ftsB	enolase/cell division protein
Reads supporting (aligned to +/- strand):  ref base . (28/32);  new base A (2/2);  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.49e-01

TGCCGTTCCGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTG  >  NC_002947/1810138‑1810424                                                                                                                                                     |                                                                                                                                             tgCCGTTCCGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTTGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑g                                                                                                                                             <  7:61202/148‑1 (MQ=255)  gCCGTTCCGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTa                                                                                                                                           >  3:56981/1‑149 (MQ=255)   ccGTTCCGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTAc                                                                                                                                          >  3:222575/1‑149 (MQ=255)   ccGTTCCGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTAc                                                                                                                                          >  4:189860/1‑149 (MQ=255)   ccGTTCCGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTAc                                                                                                                                          >  7:132866/1‑149 (MQ=255)      ttCCGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAAAGTTcagtacg                                                                                                                                         >  6:285503/1‑142 (MQ=255)        ccGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGaa                                                                                                                                      <  2:45347/148‑1 (MQ=255)        ccGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGaa                                                                                                                                      <  2:297646/148‑1 (MQ=255)             cgcgTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTcc                                                                                                                                >  3:291494/1‑149 (MQ=255)                 tttCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACg                                                                                                                            >  1:217342/1‑149 (MQ=255)                 tttCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAAAGTTCAGTAC‑GAACGTTCCTAca                                                                                                                            <  5:285503/149‑2 (MQ=255)                       aGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCCGTTTACCCCCGTCGTGCCGAGTTTCGCGCCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGATC‑GTACTGAACTTTCCTACGTACTaa                                                                                                                      >  3:78698/1‑149 (MQ=255)                         tACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACg                                                                                                                    >  5:33392/1‑149 (MQ=255)                          aCAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACgg                                                                                                                   >  3:83668/1‑149 (MQ=255)                               cAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCACGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGttt                                                                                                               <  1:241338/149‑1 (MQ=255)                                   tgctgcGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGATTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTc                                                                                                          >  6:282158/1‑149 (MQ=255)                                          cATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACg                                                                                                                            <  6:308160/124‑1 (MQ=255)                                          cATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACg                                                                                                                            >  5:308160/1‑124 (MQ=255)                                           aTCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGc                                                                                                  >  6:134088/1‑149 (MQ=255)                                             cGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAAAGTTCAGTAC‑GAACGTTCCTAc                                                                                                                             <  5:50719/120‑1 (MQ=255)                                             cGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAAAGTTCAGTAC‑GAACGTTCCTAc                                                                                                                             >  6:50719/1‑120 (MQ=255)                                              gAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTg                                                                                               <  8:132866/149‑1 (MQ=255)                                                ccagCCACTCGGCGCCAAACCCGTCTACCGCGGTCGCGCCGAGTTCCGCGGCTAAGCAAGAGATGGTAAAAACGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGc                                                                                             <  4:230437/147‑1 (MQ=255)                                                  agCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCt                                                                                           >  7:60556/1‑149 (MQ=255)                                                   gCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACgg                                                                                                                   >  2:77524/1‑124 (MQ=255)                                                   gCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACgg                                                                                                                   <  1:77524/124‑1 (MQ=255)                                                   gCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTc                                                                                          <  1:221685/149‑1 (MQ=255)                                                   gCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTc                                                                                          <  2:107426/149‑1 (MQ=255)                                                    cAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTATGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTAc                                                                                                                             <  1:283550/113‑1 (MQ=255)                                                    cAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTAc                                                                                                                             >  2:283550/1‑113 (MQ=255)                                                    cAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCg                                                                                         >  6:66607/1‑149 (MQ=255)                                                     aaCTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCgg                                                                                        >  5:56186/1‑149 (MQ=255)                                                       cTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGcc                                                                                      <  3:56620/149‑1 (MQ=255)                                                         gggCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTcc                                                                                                                                >  6:164817/1‑105 (MQ=255)                                                         gggCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTcc                                                                                                                                <  5:164817/105‑1 (MQ=255)                                                          ggCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAgg                                                                                   <  6:56186/149‑1 (MQ=255)                                                           gcgcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCACGc                                                                                  >  4:236863/1‑149 (MQ=255)                                                            cgcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCt                                                                                 <  5:257830/149‑1 (MQ=255)                                                            cgcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGCt                                                                                 <  4:291494/149‑1 (MQ=255)                                                             gcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCtt                                                                                <  2:217342/149‑1 (MQ=255)                                                                 aaGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTg                                                                            <  3:81836/149‑1 (MQ=255)                                                                  aGCGGTTTACCGCGGTCGTGCCGAGTTTCTCTGCTAAGCAAGTGTTGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGTAGGCTGGCGTCGGCCAGGCTTCGt                                                                             <  4:78698/147‑1 (MQ=255)                                                                  aGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGc                                                                           <  6:33392/149‑1 (MQ=255)                                                                       tttACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGttt                                                                                                               >  1:21479/1‑108 (MQ=255)                                                                       tttACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGttt                                                                                                               <  2:21479/108‑1 (MQ=255)                                                                         tACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCg                                                                     <  3:189860/148‑1 (MQ=255)                                                                                   gTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGaa                                                           <  4:83668/148‑1 (MQ=255)                                                                                    tGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGt                                                         <  5:66607/149‑1 (MQ=255)                                                                                      ccGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTaa                                                       >  7:99720/1‑149 (MQ=255)                                                                                      ccGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCCCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTa                                                        <  5:134088/148‑1 (MQ=255)                                                                                       cGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTaa                                                       >  3:220791/1‑148 (MQ=255)                                                                                             tCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTAc                                                >  2:194122/1‑149 (MQ=255)                                                                                                cgGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTAcgcg                                             <  4:152821/149‑1 (MQ=255)                                                                                                        cAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑g                                                                                                                                             >  1:248943/1‑45 (MQ=255)                                                                                                        cAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑g                                                                                                                                             <  2:248943/45‑1 (MQ=255)                                                                                                              aTGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAAc                               >  5:274507/1‑149 (MQ=255)                                                                                                                  tAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTgg                           <  1:194122/149‑1 (MQ=255)                                                                                                                    aaaaGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGa                          <  4:220791/148‑1 (MQ=255)                                                                                                                      aaGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCCGCCAGGCTTCGTGCTATTCGAGGACCCTAAGTAACAATTACGCGGCAGACTTTTTAACCGGGAGaa                       >  2:201021/1‑149 (MQ=255)                                                                                                                           agca‑gcCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGa                  <  8:60556/149‑1 (MQ=255)                                                                                                                           agca‑gcCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGa                  <  8:99720/149‑1 (MQ=255)                                                                                                                              a‑gcCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATgc               <  6:274507/149‑1 (MQ=255)                                                                                                                                   ggAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAg           >  8:88996/1‑149 (MQ=255)                                                                                                                                      gTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTcc        <  4:222575/149‑1 (MQ=255)                                                                                                                                          tAGGAAAGTTCAGTAC‑GAACGTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTAt    >  5:277532/1‑149 (MQ=255)                                                                                                                                            ggAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTAc                                                >  6:234296/1‑103 (MQ=255)                                                                                                                                            ggAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTAc                                                <  5:234296/103‑1 (MQ=255)                                                                                                                                            ggAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTg  >  6:76988/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                             TGCCGTTCCGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTG  >  NC_002947/1810138‑1810424

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