Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,415,032	T→C	5.7%	intergenic (+141/+27)	PP_5644 → / ← PP_3899	membrane protein/hicB family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,415,032	0	T	C	5.7%	208.4 / 2.6	70	intergenic (+141/+27)	PP_5644/PP_3899	membrane protein/hicB family protein
Reads supporting (aligned to +/- strand):  ref base T (33/33);  new base C (2/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GTAGATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACATT  >  NC_002947/4414888‑4415178                                                                                                                                                 |                                                                                                                                                     gTAGATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCTGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑a                                                                                                                                                <  8:129530/149‑1 (MQ=255)      tGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCg                                                                                                                                           >  6:33427/1‑149 (MQ=255)       ggCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCgg                                                                                                                                          >  1:321142/1‑149 (MQ=255)        gCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCACTG‑C‑AATGCggg                                                                                                                                         <  5:122769/149‑1 (MQ=255)        gCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTATCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCACTG‑C‑AATGCggg                                                                                                                                         <  6:154110/149‑1 (MQ=255)                gCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCtt                                                                                                                                 >  4:171967/1‑149 (MQ=255)                gCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCtt                                                                                                                                 >  6:36260/1‑149 (MQ=255)                gCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCAAGCAATCACGCTTTTTTCACGACCTCCCCTGCACAGTGAAGGATCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCtt                                                                                                                                 >  4:267860/1‑149 (MQ=255)                 cACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCttt                                                                                                                                <  5:65663/149‑1 (MQ=255)                   ctcctcTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTgg                                                                                                                              <  7:39211/149‑1 (MQ=255)                         ttGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGggcggc                                                                                                                         >  8:171964/1‑148 (MQ=255)                           gAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTAc                                                                                                                      >  6:317517/1‑149 (MQ=255)                             aGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTAcgc                                                                                                                    <  2:52252/149‑1 (MQ=255)                                  ccATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGtgat                                                                                                               >  6:289318/1‑149 (MQ=255)                                    aTCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGtgatga                                                                                                             <  7:171964/149‑1 (MQ=255)                                      cAAGTGATCCCACAGTCCGACTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCCCGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGAcc                                                                                                           >  5:28856/1‑149 (MQ=255)                                        aGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGAcc                                                                                                           >  6:38521/1‑147 (MQ=255)                                         gTGATCCCACAGTCCGCCTCGTCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAg                                                                                                        <  2:321142/149‑1 (MQ=255)                                         gTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAg                                                                                                        >  8:214545/1‑149 (MQ=255)                                         gTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAg                                                                                                        >  6:96024/1‑149 (MQ=255)                                         gTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAg                                                                                                        <  2:152669/149‑1 (MQ=255)                                                cacaGTCCGCCTCGGCGTGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGggcg                                                                                                                           <  2:247912/123‑1 (MQ=255)                                                cacaGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGggcg                                                                                                                           >  1:247912/1‑123 (MQ=255)                                                 acaGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTAcgc                                                                                                                    <  3:48265/129‑1 (MQ=255)                                                 acaGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTAcgc                                                                                                                    >  4:48265/1‑129 (MQ=255)                                                  caGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTAcgcg                                                                                                                   <  5:256588/129‑1 (MQ=255)                                                  caGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTAcgcg                                                                                                                   >  6:256588/1‑129 (MQ=255)                                                    gTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTaa                                                                                             <  8:140944/149‑1 (MQ=255)                                                     tCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAg                                                                                            <  5:36260/149‑1 (MQ=255)                                                       cGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAgcg                                                                                          <  7:58353/149‑1 (MQ=255)                                                       cGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAgcg                                                                                          >  8:35540/1‑149 (MQ=255)                                                       cGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGATTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAgcg                                                                                          >  8:229369/1‑149 (MQ=255)                                                             ggcggGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCttt                                                                                      <  2:30334/147‑1 (MQ=255)                                                               cggGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCgt                                                                                  <  7:214545/149‑1 (MQ=255)                                                                gggCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCgtg                                                                                 <  8:196495/149‑1 (MQ=255)                                                                   cTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCAGAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCa                                                                              <  5:289318/149‑1 (MQ=255)                                                                    ttCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGc                                                                                    <  6:49538/142‑1 (MQ=255)                                                                    ttCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGc                                                                                    >  5:49538/1‑142 (MQ=255)                                                                    ttCACCATGCATTCACACTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCgtgt                                                                                >  5:307691/1‑146 (MQ=255)                                                                      cACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCttt                                                                                      >  7:4873/1‑138 (MQ=255)                                                                          aTGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAg                                                                                                        <  1:319133/116‑1 (MQ=255)                                                                          aTGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAg                                                                                                        >  2:319133/1‑116 (MQ=255)                                                                          aTGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCAt                                                                        >  5:261785/1‑148 (MQ=255)                                                                             cATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCaa                                                                    <  5:25149/149‑1 (MQ=255)                                                                                  aCGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGaaa                                                               >  8:204383/1‑149 (MQ=255)                                                                                  aCGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATGGTCAAGTGCGGCGTTTTATTTGGGCGGGTAGCCGTGATGACCCAGGTACACACTAACCGCTTTGCGTTTCAGATCATACAAgg                                                                  >  1:102234/1‑148 (MQ=255)                                                                                     cTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGGTGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATAc                                                            <  3:267860/149‑1 (MQ=255)                                                                                        ttttCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACtt                                                          <  6:28856/148‑1 (MQ=255)                                                                                         tttCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACttt                                                         <  7:229369/148‑1 (MQ=255)                                                                                           tCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGc                                                       <  6:261785/148‑1 (MQ=255)                                                                                              cGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTtgg                                                   <  5:96024/149‑1 (MQ=255)                                                                                               gCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTtggt                                                  <  7:204383/149‑1 (MQ=255)                                                                                                     cccTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCt                                              <  5:33427/147‑1 (MQ=255)                                                                                                            caGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTg                                     <  3:171967/149‑1 (MQ=255)                                                                                                             aGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTgg                                    <  5:317517/149‑1 (MQ=255)                                                                                                                 aaGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGc                                >  6:312516/1‑149 (MQ=255)                                                                                                                      tCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAg                             >  1:155295/1‑147 (MQ=255)                                                                                                                      tCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAg                             <  2:155295/147‑1 (MQ=255)                                                                                                                      tCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGt                           >  5:243771/1‑149 (MQ=255)                                                                                                                        aTCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGt                           <  2:102234/147‑1 (MQ=255)                                                                                                                            gCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGt                     >  3:149737/1‑149 (MQ=255)                                                                                                                            gCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGt                     >  6:314010/1‑149 (MQ=255)                                                                                                                                cccTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTg                 >  1:333268/1‑149 (MQ=255)                                                                                                                                      gTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGTGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGATGGACGCTAGGTGCAGGTCGTGCCCTAc           >  1:288452/1‑149 (MQ=255)                                                                                                                                      gTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTAc           <  6:307691/149‑1 (MQ=255)                                                                                                                                        aaGCCCGCACTG‑C‑AATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTAccc         >  7:54130/2‑149 (MQ=255)                                                                                                                                           cccGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTca      <  4:149737/149‑1 (MQ=255)                                                                                                                                            ccGCACTG‑C‑AATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTca      >  7:13815/1‑148 (MQ=255)                                                                                                                                              gCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACAt   >  2:172869/1‑149 (MQ=255)                                                                                                                                                aTTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACAtt  >  7:61895/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                                     GTAGATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTG‑C‑AGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACATT  >  NC_002947/4414888‑4415178

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