Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	107,996	T→C	5.5%	I149V (ATT→GTT)	PP_5426 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	107,996	0	T	C	5.5%	218.9 / 2.3	74	I149V (ATT→GTT)	PP_5426	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/34);  new base C (2/2);  total (38/36)
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.03e-01

AGCGCTTCCACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAGGCAGTCCAGGCCTACATGCCCT  >  NC_002947/107849‑108144                                                                                                                                                      |                                                                                                                                                     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                                                                                                                                             <  3:192283/149‑1 (MQ=255)         cACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACtggtg                                                                                                                                             >  1:292537/1‑149 (MQ=255)         cACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACtggtg                                                                                                                                             >  6:255318/1‑149 (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                                                                                                            >  5:281770/1‑148 (MQ=255)                                           c‑tcGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCgg                                                                                                            >  4:135276/1‑148 (MQ=255)                                               gATAAGCCTCAACCGATTCAACGGCAAAGCACACATCATTAAGCATTTCGTCGACCGACTGGAGAACCACCACATCCAGCTCGTTGGC‑AAAGCTTTCGCCAAACAACTGGTGGGCAATCAGGTTGGCGACAGTGCTGTCGACCGGGAt                                                                                                         >  7:16198/1‑148 (MQ=255)                                                     ccTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGgc                                                                                                  <  5:255318/149‑1 (MQ=255)                                                     ccTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGgc                                                                                                  >  1:273657/1‑149 (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:267107/149‑1 (MQ=255)                                                                            ctctCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACaa                                                                           <  2:273657/149‑1 (MQ=255)                                                                            ctctCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACaa                                                                           <  6:137428/149‑1 (MQ=255)                                                                               tcatcaTTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTa                                                                              >  8:128116/1‑143 (MQ=255)                                                                               tcatcaTTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTa                                                                              <  7:128116/143‑1 (MQ=255)                                                                               tcatcaTTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGt                                                                         >  6:196512/1‑148 (MQ=255)                                                                                catcatTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACg                                                                                                    >  6:308408/1‑120 (MQ=255)                                                                                catcatTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACg                                                                                                    <  5:308408/120‑1 (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)                                                                                                                                        aaaGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAGGCAGTc                  >  5:21345/1‑147 (MQ=255)                                                                                                                                           gCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAGGCAGTCCAGGc             <  6:132027/149‑1 (MQ=255)                                                                                                                                             tttCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAGGCAGTCCAGGcc            >  6:39558/1‑148 (MQ=255)                                                                                                                                                cGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAGGCAGTCCAGGCCTAc         <  8:172089/148‑1 (MQ=255)                                                                                                                                                 gCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAGGCAGTCCAGGCCTACAt       >  8:178602/1‑149 (MQ=255)                                                                                                                                                  ccAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAGGCAGTCCAGGCCTACATg      >  4:204540/1‑149 (MQ=255)                                                                                                                                                   cAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAGGAAGTCCAGGCCTACCTGc     >  8:66349/1‑149 (MQ=255)                                                                                                                                                      tCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAGGCAGTCCAGGCCTACATGCCCt  >  5:244938/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     AGCGCTTCCACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTC‑TCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGC‑AAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAGGCAGTCCAGGCCTACATGCCCT  >  NC_002947/107849‑108144

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