Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,713,824	A→G	5.8%	I230I (ATT→ATC)	paaK ←	phenylacetate‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,713,824	0	A	G	5.8%	203.4 / 2.5	69	I230I (ATT→ATC)	paaK	phenylacetate‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (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 = 7.82e-01

GTGGCAGGGTCGATGATTTCGGGATAGAAGTGGTCTTCCCAGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTAGGACGGGGTGACCATGATGATGTCG  >  NC_002947/3713678‑3713965                                                                                                                                                   |                                                                                                                                               gTGGCAGGGTCGATGATTTCGGGATAGAAGTGGTCTTCCCAGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAAt                                                                                                                                              <  7:106639/149‑1 (MQ=255) gTGGCAGGGTCGATGATTTCGGGATAGAAGTGGTCTTCCCAGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCCGCTCGAt                                                                                                                                              >  2:153863/1‑149 (MQ=255)  ttgCAGGGTCGATGATTTCGTGATAGAAGTGGTCTTCCCAGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATc                                                                                                                                             <  2:239444/147‑1 (MQ=255)  tGGCAGGGTCGATGATTTCGGGATAGAAGTGGTCTTCCCAGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATc                                                                                                                                             >  3:317190/1‑149 (MQ=255)      aGGGTGGATTATTTCTGGATATAAGTGGTCTTCCCATATGGTCGGGCCGTCCTTTGTCTCTATGCACTCCATTGCCACCCATGTCCCCATGATTTCCGAAATGCCGTAGATATCGAGGGCATTGATGCCCATACTCTGCTCAATCGAgc                                                                                                                                         <  5:226621/149‑1 (MQ=255)         gTCGATGATTTCGGGATAGAAGTGGTCTTCCCAGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAgcggc                                                                                                                                      <  4:220352/149‑1 (MQ=255)          tCGATGATTTCGGGATAGAAGTGGTCTTCCCAGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAgcggcg                                                                                                                                     >  6:180997/1‑149 (MQ=255)           cGATGATTTCGGGATAGAAGTGGTCTTCCCAGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAgcggcg                                                                                                                                     >  1:151668/1‑148 (MQ=255)           cGATGATTTCGGGATAGAAGTGGTCTTCCCAGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAgcggag                                                                                                                                     >  1:157385/1‑146 (MQ=255)             atgatTTCGGGCTAGAAGTGGTCTTCCCAGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCCGCTCGATTGAgcagcgtat                                                                                                                                  <  1:117259/149‑4 (MQ=255)              tgatTTCGGGATAGAAGTGGTCTTCCCAGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCCGCTCGATTGAgcagcgtatg                                                                                                                                 <  1:153863/149‑5 (MQ=255)                      ggATAGAAGTGGTCTTCCCAGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGt                                                                                                                         <  5:180997/149‑1 (MQ=255)                                 gTCTTCCCAGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCgg                                                                                                              >  4:122095/1‑149 (MQ=255)                                      cccAGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCAcc                                                                                                          >  1:397936/1‑148 (MQ=255)                                      cccAGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCAc                                                                                                           <  8:199570/147‑1 (MQ=255)                                        cAGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAgcgg                                                                                                                                       >  7:382140/1‑117 (MQ=255)                                        cAGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAgcgg                                                                                                                                       <  8:382140/117‑1 (MQ=255)                                         aGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCaaa                                                                                                       >  5:400760/1‑148 (MQ=255)                                         aGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACACCGGGCCCCATGATTTCCGAAAGGCCATAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCaaa                                                                                                       >  1:76861/1‑148 (MQ=255)                                              gTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑c                                                                                                 <  6:400760/148‑1 (MQ=255)                                                   gCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGgc                                                                                           <  4:317190/149‑1 (MQ=255)                                                   gCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGgc                                                                                           >  5:67145/1‑149 (MQ=255)                                                        ccTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAgg                                                                                      >  1:350271/1‑148 (MQ=255)                                                               ctctATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGt                                                                               >  1:159418/1‑149 (MQ=255)                                                                  tATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCAc                                                                                                                     >  6:111332/1‑109 (MQ=255)                                                                  tATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGt                                                                            >  2:364502/1‑149 (MQ=255)                                                                  tATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGt                                                                            >  2:57636/1‑149 (MQ=255)                                                                    tGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGc                                                                          <  5:364406/149‑1 (MQ=255)                                                                       aCTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCggg                                                                       >  1:405099/1‑149 (MQ=255)                                                                         tCCATTGCCACCCC‑GGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTc                                                                     >  2:44794/1‑148 (MQ=255)                                                                           cATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAAACGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGACGa                                                                   >  5:64027/1‑149 (MQ=255)                                                                               gCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGcc                                                               >  5:313494/1‑149 (MQ=255)                                                                                ccACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGccc                                                              >  7:152760/1‑149 (MQ=255)                                                                                ccACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATACCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGccc                                                              >  2:120759/1‑149 (MQ=255)                                                                                  aCCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTg                                                            >  3:383732/1‑149 (MQ=255)                                                                                           cccATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGAt                                                   >  8:238417/1‑149 (MQ=255)                                                                                             cATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGAttt                                                 <  2:397936/149‑1 (MQ=255)                                                                                             cATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGAGCAGCTTGAGGTCGTGCGAGTCGATGCCCTGGCGCTCGaat                                                  >  8:146697/1‑146 (MQ=255)                                                                                              aTGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTtc                                                <  8:152760/149‑1 (MQ=255)                                                                                               tGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTtca                                               <  1:44794/149‑1 (MQ=255)                                                                                                  tttCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTc                                                                              <  2:400546/115‑1 (MQ=255)                                                                                                  tttCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTc                                                                              >  1:400546/1‑115 (MQ=255)                                                                                                  tttCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTtcatc                                             <  1:57636/148‑1 (MQ=255)                                                                                                  tttCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTtcat                                              >  5:391224/1‑147 (MQ=255)                                                                                                    tCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGc                                          >  3:368127/1‑149 (MQ=255)                                                                                                        aaaGGCCGTAGATATCGAGGGCTTTGATTCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAgg                                      <  6:64027/149‑1 (MQ=255)                                                                                                           ggCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTg                                   <  6:391224/149‑1 (MQ=255)                                                                                                           ggCCATAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTACACGGATCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGGCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTg                                   <  2:76861/149‑1 (MQ=255)                                                                                                            gCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGa                                  <  4:383732/149‑1 (MQ=255)                                                                                                               gTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCa                               >  7:292400/1‑149 (MQ=255)                                                                                                               gTAGATATCGAGGGCATTGATACCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATGCCCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGc                                <  1:120759/149‑1 (MQ=255)                                                                                                                  gATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGGCGTGCGGGGCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATg                             <  5:408473/148‑1 (MQ=255)                                                                                                                    tatCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTa                           <  1:364502/148‑1 (MQ=255)                                                                                                                          gggCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTAGGACggg                    <  2:151668/149‑1 (MQ=255)                                                                                                                          gggCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTAGGACggg                    <  6:313494/149‑1 (MQ=255)                                                                                                                          gggCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTAGGACggg                    <  6:67145/149‑1 (MQ=255)                                                                                                                           ggCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTAGGACggg                    <  2:157385/148‑1 (MQ=255)                                                                                                                              aTTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTAGGACGGGGt                  >  3:285565/1‑147 (MQ=255)                                                                                                                              aTTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTAGGACGGGGt                  <  4:285565/147‑1 (MQ=255)                                                                                                                                     cccATACGCTGCTCAATCGAGCGGCGCAGCTCATGGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTAGGACGGGGTGACCatgat         <  7:146697/149‑1 (MQ=255)                                                                                                                                     cccATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTAGGACGGGGTGACCatgat         <  8:292400/149‑1 (MQ=255)                                                                                                                                     cccATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGGCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTAGGACGGGGTGACCatga          <  5:392403/148‑1 (MQ=255)                                                                                                                                      ccgtACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTAGGACGGGGTGAcc              <  8:251010/140‑1 (MQ=255)                                                                                                                                      ccATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCATGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTAGGAGGGGGTGACCaaga          >  7:251010/1‑147 (MQ=255)                                                                                                                                       cATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTAGGACGGGGTGACCatgatga       >  3:363790/1‑149 (MQ=255)                                                                                                                                         tACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTAGGACGGGGTGACCatgatgatg     >  8:293809/1‑149 (MQ=255)                                                                                                                                          gcgccgcTCGATTGAGCAGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTAGGACGGGGTGACCATGATGATGt    >  2:22620/2‑149 (MQ=255)                                                                                                                                            gcttctCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTAGGACGGGGTGACCATGATGATGTCg  <  3:122095/149‑1 (MQ=255)                                                                                                                                            gctgctCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTAGGACGGGGTGACCATGATGATGTCg  >  3:233547/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                               GTGGCAGGGTCGATGATTTCGGGATAGAAGTGGTCTTCCCAGATGGTCGGGCCGTCCTTGGTCTCTATGCACTCCATTGCCACCCCGGGCCCCATGATTTCCGAAAGGCCGTAGATATCGAGGGCATTGATGCCCATACGCTGCTCAATCGAGCGGCGCAGCTCATCGGTCCACGGCTCGGCACCAAAAATG‑CCCAGGCGCAGCTTGAGGTCGTGCGGGTCGATGCCCTGGCGCTCGATTTCATCAGCCAGGTTGAGCATGTAGGACGGGGTGACCATGATGATGTCG  >  NC_002947/3713678‑3713965

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑