Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,740,842	T→A	19.5%	intergenic (+79/+674)	gltA → / ← yeiW	citrate synthase/putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,740,842	0	T	A	19.5%	124.1 / 29.3	67	intergenic (+79/+674)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (27/26);  new base A (6/7);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.55e-01

GGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCAC  >  NC_002947/4740695‑4740988                                                                                                                                                     |                                                                                                                                                     ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGc                                                                                                                                                    >  5:344767/1‑149 (MQ=255) ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGc                                                                                                                                                    >  5:4907/1‑149 (MQ=255)  gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGcc                                                                                                                                                   >  8:348315/1‑149 (MQ=255)        cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGaa                                                                                                                                           <  5:314787/149‑1 (MQ=255)         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAACGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGaa                                                                                                                                           <  1:328402/149‑1 (MQ=255)         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAag                                                                                                                                          <  2:292235/149‑1 (MQ=255)          cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAaga                                                                                                                                         <  3:83793/149‑1 (MQ=255)          cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCACTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAaga                                                                                                                                         <  5:13022/149‑1 (MQ=255)           gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGAagag                                                                                                                                        <  5:58332/149‑1 (MQ=38)           gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAagag                                                                                                                                        <  2:69553/149‑1 (MQ=255)           gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAagag                                                                                                                                        <  7:108256/149‑1 (MQ=255)              aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGcc                                                                                                                                     >  4:87251/1‑149 (MQ=255)                     cacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGcaca                                                                                                                               <  6:344767/148‑1 (MQ=255)                      acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGaa                                                                                                                                           <  8:145575/135‑1 (MQ=38)                      acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGaa                                                                                                                                           >  7:145575/1‑135 (MQ=38)                       ccGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGAagag                                                                                                                                        >  5:356868/1‑137 (MQ=38)                       ccGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGAagag                                                                                                                                        <  6:356868/137‑1 (MQ=38)                                aGACGCGCGAAATCGTGGCCGTTAAACCCCCCACTAACGCGGCAGTCCCCCAACCCAAC‑ACGGGTGCCCGCGGGCAG‑CTTTTTTTATTTGGCCTGGGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACaaaa                                                                                                                   <  2:351734/149‑1 (MQ=255)                                aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGAAGAGGc                                                                                                                                      >  6:374177/1‑130 (MQ=37)                                aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGAAGAGGc                                                                                                                                      <  5:374177/130‑1 (MQ=37)                                 gAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCCTTTTTTT‑TTTGTCCTGTGCCGGCCTCTTCGCGGCCATGCC‑‑CGCGAAGCGGCCGGCCCCGACACCCCAAAAc                                                                                                                  >  3:74225/1‑149 (MQ=255)                                     cgcgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCt                                                                                                                >  5:419855/1‑147 (MQ=255)                                        gACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGACAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGtt                                                                                                           >  1:328995/1‑149 (MQ=255)                                            tCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGAaga                                                                                                                                         <  2:386189/115‑1 (MQ=37)                                            tCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGAaga                                                                                                                                         >  1:386189/1‑115 (MQ=37)                                                   ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCACAGATTTTATccgccc                                                                                                >  8:161000/1‑149 (MQ=255)                                                      gAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAAc                                                                                                                  <  5:195303/128‑1 (MQ=39)                                                      gAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAAc                                                                                                                  >  6:195303/1‑128 (MQ=39)                                                         aGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGGGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCAc                                                                                          <  7:161000/149‑1 (MQ=255)                                                                      gCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGCCCCGCACAGACCCCACAAAACCTCac                                                                                                             >  5:420508/1‑116 (MQ=255)                                                                        tGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGGCCAGCCACCACAAACCCTCAGTTTTTCTCCGCCCCCTCACGCAGCCCTGTAAggg                                                                           >  3:300727/1‑149 (MQ=255)                                                                         gCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCCGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGCCACCCCAAAACCTCAGTTTTTATCCCCCCGCTCACGCAGCCCGTTCAggg                                                                           >  6:157603/1‑148 (MQ=255)                                                                            gTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                                                                       <  6:419855/149‑1 (MQ=255)                                                                              aGCCGAACGCAAG‑AAGGCTGCCCTCGGGCATCC‑TTTTTTATTTGTCCTGTGCCGGCATCTTCGCGGCCTTGCCCTCGCTATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                                                                       <  1:126550/149‑1 (MQ=255)                                                                                  gAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCGCACGCAGCCCGTTCAGGGAATTGAACgg                                                                 >  5:374007/1‑149 (MQ=255)                                                                                  gAACGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACATAGCGGGCATGGC‑‑CGCGAATAGGGCGGGCCCGAAAACCCAAAACCACACTTTTTATACCCCCCCGCACCCCCCCCGTTAAAGGTATTTAAAgg                                                                 >  2:340994/1‑149 (MQ=255)                                                                                   aaCGCAAG‑AAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTGGCGGCCTTGCC‑‑CGCGATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                                                                <  7:348315/149‑1 (MQ=255)                                                                                       taaaaaagGCTGCCCGAGGGCAGCC‑TTCTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCa                                                             <  2:95878/144‑1 (MQ=255)                                                                                        aag‑aagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCACAGTTTTTATACGCCCGCTCACGCAGCCCTTTCAGGGTATTTAACGGCGCATc                                                           >  6:363325/1‑149 (MQ=255)                                                                                             agGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                                                              <  6:169664/126‑1 (MQ=39)                                                                                             agGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                                                              >  5:169664/1‑126 (MQ=39)                                                                                               gcggCCCTCGGGAAG‑CTTTTTTTATTTGGCCTGGGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCcaccac                                                     <  6:374007/146‑1 (MQ=255)                                                                                                   cccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAATGGCGCATCCACCACGAAc                                                 >  6:179483/1‑149 (MQ=255)                                                                                                     cTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAATGGCGCATCCACCACGAACtt                                               <  5:179483/149‑1 (MQ=255)                                                                                                       ctatctgCCTTTTTTT‑TTTGTCCTGTGCGGGCCTCTTCGCGGCCTTGCC‑‑CGGGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGTTATTGAACGGCTCATCCACCACGAACttg                                              <  1:263214/142‑1 (MQ=255)                                                                                                         ggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg                                           <  2:63652/149‑1 (MQ=255)                                                                                                                   ttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATc                                  >  1:308012/1‑149 (MQ=255)                                                                                                                               tgTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTa                      >  5:300724/2‑149 (MQ=255)                                                                                                                                 tGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTac                     <  6:406545/148‑1 (MQ=255)                                                                                                                                   ccGGCCTCTACGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACTGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTaccac                  <  7:192418/149‑1 (MQ=255)                                                                                                                                    cgcccTCTTCGCGGCCTGCCC‑‑CGCGATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTaccacc                 <  3:87251/146‑1 (MQ=255)                                                                                                                                         tctTCGCGGGCAAGGC‑‑CGCGAAGAGGCCGGCACAGACCCCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGCACTTGTTGGACACCATCACCGGCAAGCTACCACCTGGCt            >  1:168214/1‑149 (MQ=255)                                                                                                                                         tctTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCCGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt            >  3:311544/1‑149 (MQ=255)                                                                                                                                         tctTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt            >  5:156556/1‑149 (MQ=255)                                                                                                                                          ctTCGCGGGCAAGGC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTac                     >  4:100715/1‑139 (MQ=255)                                                                                                                                          ctTCGCGGGCAAGGC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTac                     <  3:100715/139‑1 (MQ=255)                                                                                                                                          ctTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc           >  1:268226/1‑149 (MQ=255)                                                                                                                                           ttCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTg               >  3:332430/1‑144 (MQ=255)                                                                                                                                           ttCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTg               <  4:332430/144‑1 (MQ=255)                                                                                                                                           ttCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCg          >  3:285681/1‑149 (MQ=255)                                                                                                                                             cgcgGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCg          <  4:300727/147‑1 (MQ=255)                                                                                                                                              gcgGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAATACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgt        <  1:340994/148‑1 (MQ=255)                                                                                                                                              gcgGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgtg       >  3:395375/1‑149 (MQ=255)                                                                                                                                              gcgGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgt        <  5:363325/148‑1 (MQ=255)                                                                                                                                               cgGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                                                                        <  3:193144/83‑1 (MQ=255)                                                                                                                                               cgGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                                                                        >  4:193144/1‑83 (MQ=255)                                                                                                                                                   cTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCAc  >  3:380480/1‑149 (MQ=255)                                                                                                                                                    ttGCC‑‑CGCGAAGAGGCCTGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACtt                                               <  5:256271/103‑1 (MQ=255)                                                                                                                                                    ttGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACtt                                               >  6:256271/1‑103 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAG‑AAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCAC  >  NC_002947/4740695‑4740988

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