Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,810,280	C→A	6.7%	intergenic (+41/‑128)	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,280	0	C	A	6.7%	222.1 / 4.0	69	intergenic (+41/‑128)	eno/ftsB	enolase/cell division protein
Reads supporting (aligned to +/- strand):  ref base C (30/34);  new base A (3/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CTGTGCCGTTCCGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTG  >  NC_002947/1810135‑1810424                                                                                                                                                   |                                                                                                                                                  cTGTGCCGTTCCGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGtt                                                                                                                                                >  2:241338/1‑149 (MQ=255)    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                                                                                                                             >  6:50719/1‑120 (MQ=255)                                                cGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAAAGTTCAGTAC‑GAACGTTCCTAc                                                                                                                             <  5:50719/120‑1 (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                                                                                                                                <  5:164817/105‑1 (MQ=255)                                                            gggCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTcc                                                                                                                                >  6:164817/1‑105 (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)                                                                                                                                                   |                                                                                                                                                  CTGTGCCGTTCCGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCA‑GCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTG  >  NC_002947/1810135‑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: ‑