Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	107,981:1	+G	5.8%	coding (460/642 nt)	PP_5426 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	107,981	1	.	G	5.8%	284.0 / 10.3	69	coding (460/642 nt)	PP_5426	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (33/32);  new base G (2/2);  total (35/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 = 9.87e-01

GCGACAGGCACGGCCAGCGCTTCCACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAGGCAGTC  >  NC_002947/107834‑108128                                                                                                                                                      |                                                                                                                                                    gcgACAGGCACGGCCAGCGCTTCCACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGc                                                                                                                                                      >  5:137428/1‑148 (MQ=255)     cAGGCACGGCCAGCGCTTCCACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGc                                                                                                                                                <  5:304819/149‑1 (MQ=255)        gCACGGCCAGCGCTTCCACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCttt                                                                                                                                             >  5:187692/1‑149 (MQ=255)        gCACGGCCAGCGCTTCCACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCttt                                                                                                                                             >  6:30450/1‑149 (MQ=255)                aGCGCTTCCACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAAt                                                                                                                                      >  8:235494/1‑148 (MQ=255)                      tCCACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCACCGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCAGCACCTCATCCACGAGGTTGGC‑AAAGATTTCGCCAATc                                                                                                                                     >  4:38296/1‑143 (MQ=255)                      tCCACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATc                                                                                                                                     <  3:38296/143‑1 (MQ=255)                        cACGCCGCTCCAGTTTCCTGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCTATTCAACTGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGTTGATTGGCGAAAGCTTT‑GCCAACCAGCtgga                                                                                                                              <  3:35815/148‑2 (MQ=37)                        cACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACtggtg                                                                                                                             >  1:292537/1‑149 (MQ=255)                        cACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACtggtg                                                                                                                             >  6:255318/1‑149 (MQ=255)                        cACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACtggtg                                                                                                                             <  3:192283/149‑1 (MQ=255)                         aCGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACtggtgg                                                                                                                            >  6:217116/1‑149 (MQ=255)                          cgccgcGCCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTggg                                                                                                                           <  4:126400/149‑1 (MQ=255)                            ccgcGCGAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGTCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCAGCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCa                                                                                                                         <  7:133915/149‑1 (MQ=255)                              gcgcCAGTTGCCGGGCTGCTCAGCTTCTCTTCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCaa                                                                                                                        >  1:25401/1‑149 (MQ=255)                               cgcCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATc                                                                                                                      <  5:217116/149‑1 (MQ=255)                                 ccAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAg                                                                                                                    <  2:25401/149‑1 (MQ=255)                                   aGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGt                                                                                                                  >  4:307227/1‑149 (MQ=255)                                        ccGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGc                                                                                                              >  6:100885/1‑148 (MQ=255)                                         cGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGTTGATTGGCGAAAGCTTT‑GCCAACCAGCtggatgagg                                                                                                                         >  7:233945/1‑130 (MQ=37)                                         cGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGTTGATTGGCGAAAGCTTT‑GCCAACCAGCtggatgagg                                                                                                                         <  8:233945/136‑7 (MQ=37)                                         cGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGa                                                                                                            >  2:65173/1‑149 (MQ=255)                                           ggCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGAcc                                                                                                          <  1:65173/149‑1 (MQ=255)                                                ctCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGATGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCgtggt                                                                                                     >  7:143869/1‑149 (MQ=255)                                                 tCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCgtggtg                                                                                                    >  1:193835/1‑149 (MQ=255)                                                   aGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGgtgt                                                                                                   <  5:30450/148‑1 (MQ=255)                                                          c‑tcGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCgg                                                                                            >  4:135276/1‑148 (MQ=255)                                                          c‑tcGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCgg                                                                                            >  5:281770/1‑148 (MQ=255)                                                              gATAAGCCTCAACCGATTCAACGGCAAAGCACACATCATTAAGCATTTCGTCGACCGACTGGAGAACCACCACATCCAGCTCGTTGGC‑AAAGCTTTCGCCAAACAACTGGTGGGCAATCAGGTTGGCGACAGTGCTGTCGACCGGGAt                                                                                         >  7:16198/1‑148 (MQ=255)                                                                    ccTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGgc                                                                                  >  1:273657/1‑149 (MQ=255)                                                                    ccTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGgc                                                                                  <  5:255318/149‑1 (MQ=255)                                                                         aCCGATTCAACGGCAAAGCTCTCAACATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGctgctg                                                                             >  7:200480/1‑149 (MQ=255)                                                                          ccGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATc                                                                                                                                     >  1:267922/1‑92 (MQ=255)                                                                          ccGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATc                                                                                                                                     <  2:267922/92‑1 (MQ=255)                                                                          ccGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTgg                                                                            <  6:187692/149‑1 (MQ=255)                                                                               tCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCGTCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCTTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTg                                                                       >  7:117730/1‑149 (MQ=255)                                                                                 aaCGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTgg                                                                      <  2:193835/148‑1 (MQ=255)                                                                                    ggCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTggagga                                                                  >  7:101479/1‑149 (MQ=255)                                                                                       aaaGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGt                                                               <  3:135276/149‑1 (MQ=255)                                                                                         aGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTa                                                              >  5:29234/1‑148 (MQ=255)                                                                                           ctctCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACaa                                                           <  2:273657/149‑1 (MQ=255)                                                                                           ctctCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACaa                                                           <  2:267107/149‑1 (MQ=255)                                                                                           ctctCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACaa                                                           <  6:137428/149‑1 (MQ=255)                                                                                              tcatcaTTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTa                                                              <  7:128116/143‑1 (MQ=255)                                                                                              tcatcaTTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTa                                                              >  8:128116/1‑143 (MQ=255)                                                                                              tcatcaTTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGt                                                         >  6:196512/1‑148 (MQ=255)                                                                                               catcatTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACg                                                                                    <  5:308408/120‑1 (MQ=255)                                                                                               catcatTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACg                                                                                    >  6:308408/1‑120 (MQ=255)                                                                                                atcatTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCAt                                                      >  3:297880/1‑149 (MQ=255)                                                                                                   atTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGaa                                                                          >  3:271459/1‑126 (MQ=255)                                                                                                   atTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGaa                                                                          <  4:271459/126‑1 (MQ=255)                                                                                                     tAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCAcc                                                 >  3:111547/1‑149 (MQ=255)                                                                                                          aTTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTc                                            <  7:266609/149‑1 (MQ=255)                                                                                                           tttCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTc                                            >  5:233692/1‑148 (MQ=255)                                                                                                                   aCCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCg                                   <  1:292472/149‑1 (MQ=255)                                                                                                                    ccGAGTGCTGCACCACCTCATCCACCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCACGTTGGCGACAGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCCGATGTACAAGTCATTCACCTCTTCCACCAGGCgg                                  >  2:134460/1‑149 (MQ=255)                                                                                                                       aGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTgc                               <  4:111547/149‑1 (MQ=255)                                                                                                                           gtgCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCa                           <  6:29234/149‑1 (MQ=255)                                                                                                                                ccaccTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAg                      >  1:226431/1‑149 (MQ=255)                                                                                                                                  accTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGc                    <  8:143869/149‑1 (MQ=255)                                                                                                                                   ccTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCt                   <  7:270532/149‑1 (MQ=255)                                                                                                                                       aTCCAGCTGGTTGGC‑AAAGCTTTCGCAAATCAACTGGTTGGCAATCATGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCa                <  5:100885/148‑1 (MQ=255)                                                                                                                                        tCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATc              <  4:297880/149‑1 (MQ=255)                                                                                                                                        accAGTTGATTGGCGAAAGCTTT‑GCCAACCAGCTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATc              >  3:226586/2‑149 (MQ=38)                                                                                                                                          cAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATcc             <  2:226431/148‑1 (MQ=255)                                                                                                                                            gCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAg           >  1:159616/1‑148 (MQ=255)                                                                                                                                                gTTTGC‑AAAGCTTTCGCCAATCAACTGGTGGTGAATCAGGTTGGCGACCGTGGTGTCGAGCGTGATCAACGGCTGCTTGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAgg       <  8:16198/148‑1 (MQ=255)                                                                                                                                                gTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCTTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAgg       <  8:117730/148‑1 (MQ=255)                                                                                                                                                gTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAgg       <  7:235494/148‑1 (MQ=255)                                                                                                                                                    gC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAGGCAGTc  >  7:172089/1‑149 (MQ=255)                                                                                                                                                     c‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAGGCa     <  6:74649/145‑1 (MQ=255)                                                                                                                                                     c‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAGGCa     >  5:74649/1‑145 (MQ=255)                                                                                                                                                      |                                                                                                                                                    GCGACAGGCACGGCCAGCGCTTCCACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAGGCAGTC  >  NC_002947/107834‑108128

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