Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,845,237	C→A	10.4%	R50R (CGG→AGG)	coaD →	phosphopantetheine adenylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,845,237	0	C	A	10.4%	215.5 / 12.5	74	R50R (CGG→AGG)	coaD	phosphopantetheine adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (31/34);  new base A (3/5);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

GATGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTCCGACTTCGAGT  >  NC_002947/5845089‑5845378                                                                                                                                                       |                                                                                                                                              gATGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGaacaac                                                                                                                                               >  1:73898/1‑149 (MQ=255) gATGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGaacaac                                                                                                                                               <  1:274417/149‑1 (MQ=255)  aTGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACg                                                                                                                                              <  6:91854/149‑1 (MQ=255)     aaCCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCTAGAAAAACCCGTTGTTCCAGGG‑GGAACAAAggg                                                                                                                                            <  1:231998/148‑1 (MQ=255)     aaCCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACggg                                                                                                                                            <  8:83256/148‑1 (MQ=255)     aaCCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACggg                                                                                                                                            <  2:138886/148‑1 (MQ=255)     aaCCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCGTTGTTCCAGGG‑GGAACAAAggg                                                                                                                                            <  1:333219/148‑1 (MQ=255)      aCCGAGTGTTGTACCCGGTTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGATCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGGCCGAAGAAAAACCCTTTGTTCGCCCT‑GTAACAACGGGtt                                                                                                                                          <  4:201226/149‑1 (MQ=255)            tgttgtACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTc                                                                                                                                    >  2:101544/1‑149 (MQ=255)            tgttgtACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCGTTGTTCCAGGG‑GGAACAAAGGGtttttctt                                                                                                                                    <  2:20380/149‑7 (MQ=38)                gtACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGccc                                                                                                                                >  1:298127/1‑149 (MQ=255)                 tACCAGGGTACATTCGACCCCAATAACAAAGGCCATGGCGATCTGGTCGAGCTCGCATCGCTGCTATTCGACAACGAGATCATAGCGGTGGC‑CGACATCCCAAAGAAAATCCCATTCTTACCCCT‑GGAACAACCGGTTGAGCAAGCCCg                                                                                                                               >  8:276173/1‑149 (MQ=255)                    cctggTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGa                                                                                                                            <  4:90873/146‑1 (MQ=255)                    ccGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGt                                                                                                                              >  8:177146/1‑147 (MQ=255)                    ccGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGa                                                                                                                            <  2:73898/149‑1 (MQ=255)                            tttCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCAcc                                                                                                                     >  1:9879/1‑148 (MQ=255)                                 aacccATTACCAAAGGCCATGGCTATCTGGTCGAGCGCGCTTCGCGTCTATTCGACCACGTTATCATCGCTGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCAc                                                                                                               <  7:276173/147‑1 (MQ=255)                                 aCCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GTAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCAc                                                                                                               <  1:195680/149‑1 (MQ=255)                                   cccATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCt                                                                                                             >  5:325663/1‑149 (MQ=255)                                           cAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGtt                                                                                                                                          >  4:87480/1‑112 (MQ=255)                                           cAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGtt                                                                                                                                          <  3:87480/112‑1 (MQ=255)                                           cAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATg                                                                                                     >  3:77342/1‑149 (MQ=255)                                           cAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATg                                                                                                     >  2:267408/1‑149 (MQ=255)                                           cAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATg                                                                                                     >  2:180918/1‑149 (MQ=255)                                           cAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCGTTGTTCCAGGG‑GGAACAAAGGGtttttc                                                                                                                                      <  1:110149/116‑5 (MQ=37)                                           cAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCGTTGTTCCAGGG‑GGAACAAAGGGtttttc                                                                                                                                      >  2:110149/1‑112 (MQ=37)                                           cAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCGTTGTTCCAGGG‑GGAACAAAGGGtt                                                                                                                                          <  4:310625/112‑1 (MQ=38)                                           cAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCGTTGTTCCAGGG‑GGAACAAAGGGtt                                                                                                                                          >  3:310625/1‑112 (MQ=38)                                             aaGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTc                                                                                                   >  6:344055/1‑149 (MQ=255)                                              aGGCCATGGCGATCTGGTCGAGCGCGCTTAGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCg                                                                                                  <  4:251449/149‑1 (MQ=255)                                                gCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCg                                                                                                  >  1:282586/1‑147 (MQ=255)                                                gCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCg                                                                                                  >  2:35365/1‑147 (MQ=255)                                                gCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGaa                                                                                                <  2:282586/149‑1 (MQ=255)                                                gCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGaa                                                                                                <  4:77342/149‑1 (MQ=255)                                                 ccATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAg                                                                                               <  7:176897/149‑1 (MQ=255)                                                      gCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGATCAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATc                                                                                          <  1:101544/149‑1 (MQ=255)                                                      gCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATc                                                                                          <  7:177146/149‑1 (MQ=255)                                                              gTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATc                                                                                          >  5:19839/1‑141 (MQ=255)                                                              gTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATc                                                                                          <  6:19839/141‑1 (MQ=255)                                                              gTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTtctc                                                                                  >  1:269258/1‑149 (MQ=255)                                                                 gAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAAAGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTc                                                                               >  7:179806/1‑149 (MQ=255)                                                                  aGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTcc                                                                              >  2:116688/1‑149 (MQ=255)                                                                   gcgcgcTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTccc                                                                             <  2:298127/149‑1 (MQ=255)                                                                       gcTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGtt                                                                         >  2:116958/1‑149 (MQ=255)                                                                           cgcgGCTATTCGAAGACGTGATGATCGTGGTGGC‑GGCCAGCCCCAAGAAAAAGCCTTTGTTCGCCCT‑GGAACTACGGGTTGAGCTCGCCCGTGAGGTCACCAATCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGgc                                                                      <  3:254064/148‑1 (MQ=255)                                                                            gcgGCTATTCGACCACGTGATCATCGCGGTGGCGGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGgcg                                                                     >  1:54686/1‑149 (MQ=255)                                                                             cgGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCa                                                                   <  1:180918/149‑1 (MQ=255)                                                                               gCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTATGAAAAAC‑GGTTTAGCTCGCCCCTGAGGTCACCAAGCACCAGCCCAATTTAGAAGTCAAAGGCGTATATTACCAGTTTGCGCaat                                                                 >  1:110460/1‑147 (MQ=255)                                                                                      gACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTcc                                                                              <  2:207021/129‑1 (MQ=255)                                                                                      gACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTGGCCaa                                                          <  2:9879/149‑1 (MQ=255)                                                                                       aCCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAg                                                         <  2:352798/149‑1 (MQ=255)                                                                                            gTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGaa                                                      >  5:117847/1‑147 (MQ=255)                                                                                               atcatcGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAgg                                                  <  2:269258/148‑1 (MQ=255)                                                                                                            c‑ggcCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTc                                    <  5:344055/149‑1 (MQ=255)                                                                                                               gcCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCt                                  <  2:110460/149‑1 (MQ=255)                                                                                                                      ccAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGt                              <  5:316682/146‑1 (MQ=255)                                                                                                                      ccAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCAGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCACTGTTGGCGCATTTCGCCAAGGAACAGGGCCCCAACCTCTTCCTGCGt                              >  6:316682/1‑146 (MQ=255)                                                                                                                         aGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTcc                                                                              >  4:287993/1‑95 (MQ=255)                                                                                                                         aGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTcc                                                                              <  3:287993/95‑1 (MQ=255)                                                                                                                          gAAAAACCCTTTGTTCCCCCT‑GGAACAAGGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCtgc                       <  8:179806/149‑1 (MQ=255)                                                                                                                               aCCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCtgcgtgcg                  >  2:29945/1‑149 (MQ=255)                                                                                                                                 ccTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGgt                <  6:325663/149‑1 (MQ=255)                                                                                                                                   tttGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGgtgt              >  7:302616/1‑149 (MQ=255)                                                                                                                                      gTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTcc            >  5:329663/1‑148 (MQ=255)                                                                                                                                      gTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTcc            <  1:35365/148‑1 (MQ=255)                                                                                                                                       ttCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTcc            >  3:254373/1‑147 (MQ=255)                                                                                                                                       ttCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTcc            <  4:254373/147‑1 (MQ=255)                                                                                                                                         cccccT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTCCGAc         <  1:116958/148‑1 (MQ=255)                                                                                                                                          ccccT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGgt                <  3:128416/140‑1 (MQ=255)                                                                                                                                          ccccT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGgt                >  4:128416/1‑140 (MQ=255)                                                                                                                                          ccccT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTCCGACtt       >  3:156762/1‑149 (MQ=255)                                                                                                                                           cccT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTc             <  6:295923/142‑1 (MQ=255)                                                                                                                                           cccT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTc             >  5:295923/1‑142 (MQ=255)                                                                                                                                            ccT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTCCGACTTc      >  1:306377/1‑148 (MQ=255)                                                                                                                                             cT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTCCGACTTc      >  8:1386/1‑147 (MQ=255)                                                                                                                                                 gAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTCCGACTTCGAGt  <  1:29945/148‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                              GATGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGC‑GGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCT‑GGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTCCGACTTCGAGT  >  NC_002947/5845089‑5845378

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