Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,653,007	A→T	6.8%	intergenic (‑108/+53)	PP_4961 ← / ← PP_4962	lipoprotein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,653,007	0	A	T	6.8%	151.0 / 2.9	59	intergenic (‑108/+53)	PP_4961/PP_4962	lipoprotein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/21);  new base T (2/2);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

CCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCTTGCCGATGACCGGCACAGTCTCGGC  >  NC_002947/5652869‑5653155                                                                                                                                           |                                                                                                                                                      ttttggCCGTGACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCctt                                                                                                                                             >  4:265091/5‑149 (MQ=255)  cAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCcttc                                                                                                                                            >  2:45683/1‑149 (MQ=255)     ggCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCTGTGCGGcccttcgcgg                                                                                                                                         <  2:48849/149‑8 (MQ=38)           cACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGAc                                                                                                                                    >  2:80428/1‑148 (MQ=255)                 cAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAAt                                                                                                                             >  6:67487/1‑149 (MQ=255)                 cAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCGTCTTCAAGACCACaaa                                                                                                                              >  3:229373/1‑148 (MQ=255)                   aaaGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCa                                                                                                                           <  1:80428/149‑1 (MQ=255)                   aaaGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCa                                                                                                                           <  5:67487/149‑1 (MQ=255)                               ttcaGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTATCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCttcttc                                                                                                                                         >  5:237278/1‑123 (MQ=255)                                  aGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCa                                                                                                                           >  2:15257/1‑134 (MQ=255)                                  aGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCa                                                                                                                           <  1:15257/134‑1 (MQ=255)                                  aGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCACCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAAc                                                                                                                         >  3:108093/1‑136 (MQ=255)                                     ttCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTAt                                                                                                          >  7:181341/1‑148 (MQ=255)                                      tCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTa                                                                                                        >  2:210809/1‑149 (MQ=255)                                        cTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGTGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAc‑c                                                                                                     >  4:171618/1‑149 (MQ=255)                                         tCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAc‑cc                                                                                                    >  3:192139/1‑149 (MQ=255)                                                ggggCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGTAGCGGGCTAGCCCGCTAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCaa                                                                                              <  1:45683/148‑1 (MQ=255)                                                 gggCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAgg                                                                                            <  8:159144/149‑1 (MQ=255)                                                     tgctgGGCCGCCCGTCTCGGGCTGGCCCGCTCCCACCGGTACAGTGCTGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCaa                                                                                        <  4:229373/149‑1 (MQ=38)                                                      gctgTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCt                                                                                                            >  5:77047/1‑129 (MQ=255)                                                      gctgTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCt                                                                                                            <  6:77047/129‑1 (MQ=255)                                                         gtgCAGCCCTTCTCGGGCTCGCCCGCTCCCCCAGGTACAGTGCAGTCTCGGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGAGCACCACTCAACGTTCAGGTTGGCTATTAC‑CCGCTCTAGGGGAACCCg                                                                                    >  5:278295/1‑149 (MQ=38)                                                               cccTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCt                                                                              >  1:47424/1‑149 (MQ=255)                                                                 cTTCTCGGGCTCGCCCGCACCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTg                                                                            <  1:168411/149‑1 (MQ=255)                                                                        ggCTCGCCCGCTCCCACCGGAACAGTGCAGTCTCTGGGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCtggt                                                                     <  3:66077/149‑1 (MQ=255)                                                                         gcccccccGCTCCCACCGGTCCAGGGCCGTCTCTGGGGGCGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCtggtg                                                                    <  4:192139/144‑1 (MQ=255)                                                                         gCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCCCAGCGGCCCCGTCTTGAAGAACACAAATAAACCTTAAGGTTGGCAATTTC‑CCGATCAAGGGCAACCCGGCGCCGTGCGCggggg                                                                    >  6:52366/1‑147 (MQ=37)                                                                           tCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCAGCAATCAACGTTCAGGTTGCCTATTAC‑CCGGTCCAGGGCAACCCGGCGCCTTGCGCGggtggt                                                                  >  3:36894/1‑149 (MQ=39)                                                                           tCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTggtggt                                                                  >  6:93097/1‑149 (MQ=255)                                                                               ccGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtc                                                              >  8:57822/1‑149 (MQ=255)                                                                                   tCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtcgtcg                                                          >  1:117629/1‑149 (MQ=255)                                                                                     ccccaGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTa                                                        <  3:265091/146‑1 (MQ=255)                                                                                     ccACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCTCTGGTGGTCTTCGTCGTa                                                        >  1:208372/1‑149 (MQ=255)                                                                                        caGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCtggtg                                                                    >  1:161246/1‑134 (MQ=255)                                                                                            tACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACACCGGCCCCTTCTTCAAGAACACAAAACAACCTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCCGCGCCATTCGCTGGTGGTCTTGGTGGTACattta                                                  >  6:208788/1‑147 (MQ=255)                                                                                             aCAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTc                                                  >  2:35134/1‑147 (MQ=255)                                                                                               aGTGCAGTCTCTGTGGGAGCAGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACTCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAg                                               <  1:83830/149‑1 (MQ=255)                                                                                                 tGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGcc                                            >  3:254193/1‑149 (MQ=255)                                                                                                    aGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTc                                         >  1:309405/1‑149 (MQ=255)                                                                                                         ctGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGc                                    >  3:32816/1‑149 (MQ=255)                                                                                                          tgtgGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCgg                                      <  2:47424/146‑1 (MQ=255)                                                                                                          tgtgGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGcc                                   >  2:124956/1‑149 (MQ=255)                                                                                                            tgGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGccc                                  >  8:233950/1‑148 (MQ=255)                                                                                                            tgGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGccc                                  <  7:57822/148‑1 (MQ=255)                                                                                                            tgGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGccc                                  <  7:233950/148‑1 (MQ=255)                                                                                                            tgGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCa                                 >  2:278302/1‑149 (MQ=255)                                                                                                             gggAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCAc                                <  8:181341/149‑1 (MQ=255)                                                                                                             gggAGCGGCCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCAc                                <  5:208788/149‑1 (MQ=255)                                                                                                                       gAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCTTGCCg                      <  1:124956/149‑1 (MQ=255)                                                                                                                          cccGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCTTGCCGAt                    >  1:309095/1‑148 (MQ=255)                                                                                                                                 ggggCCGCTGTGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCt                           >  4:83911/2‑134 (MQ=255)                                                                                                                                 ggggCCGCTGTGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCt                           <  3:83911/133‑1 (MQ=255)                                                                                                                                  gggCCGCTGTGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCTTGCCGATGACCGGca            >  4:261694/1‑148 (MQ=255)                                                                                                                                   ggCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCTTGCCGATGACCGGcaca          <  5:93097/149‑1 (MQ=255)                                                                                                                                   ggCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCTTGCCGATGACCGGca            >  2:16739/1‑147 (MQ=255)                                                                                                                                     ccGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCTTGCCGa                     >  2:51810/1‑136 (MQ=255)                                                                                                                                     ccGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCTTGCCGa                     <  1:51810/136‑1 (MQ=255)                                                                                                                                      cGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTc                                                                 <  6:54788/91‑1 (MQ=255)                                                                                                                                      cGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTc                                                                 >  5:54788/1‑91 (MQ=255)                                                                                                                                         acaGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCTTGCCGATGACCGGCACAGTCTCg    >  4:39426/1‑149 (MQ=255)                                                                                                                                           aGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCTTGCCGATGACCGGCACAGTCTCGGc  <  5:52366/149‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                      CCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAC‑CCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCTTGCCGATGACCGGCACAGTCTCGGC  >  NC_002947/5652869‑5653155

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