Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,251,506	G→A	6.8%	D112D (GAC→GAT)	PP_1094 ←	electron transport complex protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,251,506	0	G	A	6.8%	172.6 / 2.7	59	D112D (GAC→GAT)	PP_1094	electron transport complex protein
Reads supporting (aligned to +/- strand):  ref base G (31/24);  new base A (2/2);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

ACCCTGGTCCCGCTTCAGCGCCCAGTGATCGCCATGATCGCGGTTGTTGAACTGACTCAGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGG  >  NC_002947/1251362‑1251652                                                                                                                                                 |                                                                                                                                                     aCCCTGGTCCCGCTTCAGCGCCCAGTGATCGCCATGATCGCGGTTGTTGAACTGACTCAGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGccgtccg                                                                                                                                                  >  1:16480/1‑149 (MQ=255)   ccTGGTCCCGCTTCAGCGCCCAGTGATCGCCATGATCGCGGTTGTTGAACTGACTCAGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCgtg                                                                                                                                                <  6:239009/149‑1 (MQ=255)        tCCCGCTTCAGCGCCCAGTGATCGCCATGATCGCGGTTGTTGAACTGACTCAGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAt                                                                                                                                           >  1:94215/1‑149 (MQ=255)         cccGCTTCAGCGCCCAGTGATCGCCATGATCGCGGTTGTTGAACTGACTCAGCCAGTGCAATCGCGGATCCCCC‑GGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑gg                                                                                                                                        >  5:61027/1‑149 (MQ=255)          ccGCTTCAGCGCCCAGTGATCGCCATGATCGCGGTTGTTGAACTGACTCAGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑gg                                                                                                                                        >  3:61020/1‑149 (MQ=255)             cTTCAGCGCCCAGTGATCGCCATGATCGCGGTTGTTGAACTGACTCAGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCTGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGa                                                                                                                                     >  8:28244/1‑149 (MQ=255)              ttCAGCGCCCAGTGATCGCCATGATCGCGGTTGTTGAACTGACTCAGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGAc                                                                                                                                    >  3:137598/1‑149 (MQ=255)                 ctcgcCCAGTGATCGCCATGATCGCGGTTGTTGAACTGACTCAGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAG‑GGCGACGCt                                                                                                                                 <  5:134995/147‑1 (MQ=255)                    gcCCAGTGATCGCCATGATCGCGGTTGTTGAACTGACTCAGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCt                                                                                                                                 >  6:134995/1‑146 (MQ=255)                       cAGTGATCGCCATGATCGCGGTTGTTGAACTGACTCAGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCGGGGCTCTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCAc                                                                                                                           >  7:234612/1‑149 (MQ=255)                                              gttgAACTGACTCAGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTgg                                                                                                     <  6:118984/148‑1 (MQ=255)                                                tgAACTGACTCAGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGcc                                                                                                  >  1:9150/1‑149 (MQ=255)                                                  aaCTGACTCAGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCgg                                                                                                                <  4:48870/133‑1 (MQ=255)                                                  aaCTGACTCAGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCgg                                                                                                                >  3:48870/1‑133 (MQ=255)                                                   aCTGACTCAGCCAGTGCAATCGCGGAT‑CCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCtgt                                                                                               <  6:61027/148‑1 (MQ=255)                                                    cTGACTCAGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCtgtg                                                                                              <  4:61020/149‑1 (MQ=255)                                                         tCAGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGAt                                                                                          >  4:113576/1‑148 (MQ=255)                                                         tCAGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGAt                                                                                          >  1:201490/1‑148 (MQ=255)                                                           aGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCTGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAg                                                                                       <  7:28244/149‑1 (MQ=255)                                                                gTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATaa                                                                                  <  2:94215/149‑1 (MQ=255)                                                                          ggATCCCCCAGGCGATCGCCCAGCCCGGGGCTCTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTgg                                                                        <  8:234612/149‑1 (MQ=255)                                                                             tCCCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCtt                                                                     <  5:89584/149‑1 (MQ=255)                                                                                cccAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTAcc                                                                  >  4:19089/1‑149 (MQ=255)                                                                                   aGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGc                                                              <  1:268368/149‑1 (MQ=255)                                                                                       gATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCa                                                            >  7:145199/1‑147 (MQ=255)                                                                                       gATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCATCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCa                                                            >  3:248840/1‑147 (MQ=255)                                                                                       gATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCATCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATa                                                          >  3:252807/1‑149 (MQ=255)                                                                                         tCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGt                                                        <  2:9150/149‑1 (MQ=255)                                                                                         tCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGt                                                        >  1:137993/1‑149 (MQ=255)                                                                                             ccAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCa                                                    <  4:137598/149‑1 (MQ=255)                                                                                                gCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACgg                                                 <  8:145199/149‑1 (MQ=255)                                                                                                 cccGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACgg                                                 >  6:216414/1‑148 (MQ=255)                                                                                                    ggggCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTg                                                       >  4:10662/1‑139 (MQ=255)                                                                                                    ggggCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTg                                                       <  3:10662/139‑1 (MQ=255)                                                                                                      ggCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAgcc                                           >  4:243741/1‑149 (MQ=255)                                                                                                        cTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCGTGCCACTATGGCA‑GGTAAgccaa                                         <  2:16480/149‑3 (MQ=255)                                                                                                           tCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCATCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCa                                       <  4:248840/148‑1 (MQ=255)                                                                                                           tCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCATCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCa                                       <  4:252807/148‑1 (MQ=255)                                                                                                             cTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGa                                    >  4:128090/1‑149 (MQ=255)                                                                                                                ttgttCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATAc                                 >  3:88194/1‑149 (MQ=255)                                                                                                                  gttCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACgg                               >  7:99071/1‑149 (MQ=255)                                                                                                                  gttCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACgg                               >  3:150021/1‑149 (MQ=255)                                                                                                                   ttCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGc                              >  8:15757/1‑149 (MQ=255)                                                                                                                      gACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTAt                           >  8:54005/1‑149 (MQ=255)                                                                                                                       aCAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTAtt                          <  7:54005/149‑1 (MQ=255)                                                                                                                         aaCCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTgg                        <  5:216414/149‑1 (MQ=255)                                                                                                                         aaCCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTgg                        <  3:243741/149‑1 (MQ=255)                                                                                                                          aCCTGGCTGCCGATCACCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCATGGGCATGTGTGATCAGGATAATCGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATAAGGCTATAggg                       >  4:206700/1‑149 (MQ=255)                                                                                                                           ccTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGc                      >  4:199054/1‑149 (MQ=255)                                                                                                                            cTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATAc                                 <  5:202385/137‑1 (MQ=255)                                                                                                                            cTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATAc                                 >  6:202385/1‑137 (MQ=255)                                                                                                                            cTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGcc                     >  4:247295/1‑149 (MQ=255)                                                                                                                            cTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGcc                     >  2:60904/1‑149 (MQ=255)                                                                                                                               tcggCCGATGAGCGGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCctct                  <  8:74782/146‑1 (MQ=255)                                                                                                                               gCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCTCGGTAAGCCGCCAGGATACGGCTATTGGGCctct                  <  1:172734/149‑1 (MQ=255)                                                                                                                                cTGCCGATCAGCCGGCCGTCCGTGGTGATGGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCctct                  >  1:165014/1‑149 (MQ=255)                                                                                                                                 tGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAAt                                                                                 >  8:215913/1‑85 (MQ=255)                                                                                                                                    cGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTgcg             <  2:137993/149‑1 (MQ=255)                                                                                                                                     gATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTgcgg            >  4:126384/1‑149 (MQ=255)                                                                                                                                      aTCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTgcggc           <  7:234848/149‑1 (MQ=255)                                                                                                                                           ccggccgTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCtgt                >  6:240895/1‑139 (MQ=255)                                                                                                                                           ccggccgTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCtgt                <  5:240895/139‑1 (MQ=255)                                                                                                                                              gccgTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTg   <  3:199054/149‑1 (MQ=255)                                                                                                                                               ccgtccgtGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTgg  <  3:126384/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     ACCCTGGTCCCGCTTCAGCGCCCAGTGATCGCCATGATCGCGGTTGTTGAACTGACTCAGCCAGTGCAATCGCGGATCCCCCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGAT‑GGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACC‑TGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGG  >  NC_002947/1251362‑1251652

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