Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,136,937	G→C	6.2%	intergenic (‑187/‑204)	fadE ← / → gstB	medium‑long chain acyl‑CoA dehydrogenase/glutathione S‑transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,136,937	0	G	C	6.2%	210.6 / 4.0	69	intergenic (‑187/‑204)	fadE/gstB	medium‑long chain acyl‑CoA dehydrogenase/glutathione S‑transferase
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base C (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 = 2.51e-01

CTTAGACGCAATGCTGTGAACGGCATGTTCAGAACAGTGACAGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTA  >  NC_002947/2136795‑2137083                                                                                                                                                 |                                                                                                                                                   cTTAGACGCAATGCTGTGAACGGCATGTTCAGAACAGTGACAGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCg                                                                                                                                              >  3:176073/1‑149 (MQ=255) cTTAGACGCAATGCTGTGAACGGCATGTTCAGAACAGTGACAGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCg                                                                                                                                              >  3:170668/1‑149 (MQ=255)    aGACGCAATGCTGTGAACGGCATGTTGAGAACAGTGACAGGGCGTGGCGGTGGGAGTTTCAGCTCAGTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑TTTGATGAATGCGGTg                                                                                                                                           >  6:362167/1‑149 (MQ=255)         cAATGCTGTGAACGGCATGTTCAGAACAGTGACAGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGACG‑CCACCGCATTCATCAAT‑CGGTggcatc                                                                                                                                      <  3:209716/149‑4 (MQ=255)                     cGGCATGTTCAGAACAGTGACAGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGaa                                                                                                                          <  3:166098/149‑1 (MQ=255)                      ggCATGTTCAGAACAGTGACAGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAAt                                                                                                                         <  4:351896/149‑1 (MQ=255)                        cATGTTCAGAACAGTGACAGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGa                                                                                                                       >  7:132358/1‑149 (MQ=255)                          tGTTCAGAACAGTGACAGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGaa                                                                                                                                                  <  4:236863/120‑1 (MQ=255)                          tGTTCAGAACAGTGACAGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGaa                                                                                                                                                  >  3:236863/1‑120 (MQ=255)                                 aaCAGTGACAGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGt                                                                                                              >  8:180624/1‑149 (MQ=255)                                    aGTGACAGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTtgg                                                                                                           >  6:177902/1‑149 (MQ=255)                                     gTGACAGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTg                                                                                                                                           <  2:99546/116‑1 (MQ=255)                                     gTGACAGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTg                                                                                                                                           >  1:99546/1‑116 (MQ=255)                                     gTGACAGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTtggc                                                                                                          <  7:152676/149‑1 (MQ=255)                                       gACAGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAAt                                                                                                                         >  7:114747/1‑132 (MQ=255)                                       gACAGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAAt                                                                                                                         <  8:114747/132‑1 (MQ=255)                                       gACAGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAAAGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTtggctg                                                                                                        >  4:241046/1‑149 (MQ=255)                                          aGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCCTCCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTtggctggc                                                                                                      <  1:279300/149‑1 (MQ=255)                                            ggCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAgg                                                                                                   >  6:306605/1‑149 (MQ=255)                                               gTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCGCATTCATCAAT‑CGGTggcatc                                                                                                                                      >  5:317792/1‑108 (MQ=38)                                               gTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCGCATTCATCAAT‑CGGTggcatc                                                                                                                                      <  6:317792/111‑4 (MQ=38)                                                   cGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGcc                                                                                            <  2:174225/149‑1 (MQ=255)                                                    ggTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCt                                                                                           <  4:170668/149‑1 (MQ=255)                                                    ggTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCt                                                                                           <  4:176073/149‑1 (MQ=255)                                                            tttCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTgg                                                                                   <  5:74373/149‑1 (MQ=255)                                                                   tCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGAcc                                                                                                                     >  2:284614/1‑108 (MQ=255)                                                                   tCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGAcc                                                                                                                     <  1:284614/108‑1 (MQ=255)                                                                   tCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAAc                                                                            >  6:307952/1‑149 (MQ=255)                                                                     aTTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAAGCCCTCCTAAAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCTGTGGCGGCCGCGACGGTGAATGACCGTGGGGTTGGCGGGCAGGCTCGGCCTATGTATTGTTaaaa                                                                             <  6:143220/146‑1 (MQ=255)                                                                     aTTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCggggg                                                                                                                >  8:136554/1‑111 (MQ=255)                                                                     aTTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCggggg                                                                                                                <  7:136554/111‑1 (MQ=255)                                                                             gCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTAtgtgt                                                                                     <  8:181107/130‑1 (MQ=255)                                                                             gCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTAtgtgt                                                                                     >  7:181107/1‑130 (MQ=255)                                                                              ccTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGaa                                                                 >  7:261401/1‑149 (MQ=255)                                                                                  tACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGt                                                                                                              >  6:231792/1‑100 (MQ=255)                                                                                     cGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGt                                                          <  5:307952/149‑1 (MQ=255)                                                                                     cGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGt                                                          <  8:261401/149‑1 (MQ=255)                                                                                     cGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGt                                                          >  8:320586/1‑149 (MQ=255)                                                                                       gCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATTAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGa                                                        >  7:129823/1‑149 (MQ=255)                                                                                       gCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTg                                                         >  8:319015/1‑148 (MQ=255)                                                                                        ccTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGAc                                                       >  7:327168/1‑149 (MQ=255)                                                                                          tATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCT‑GCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGt                                                     >  4:102702/1‑148 (MQ=255)                                                                                             tCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGTCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTg                                                    >  3:202546/1‑147 (MQ=255)                                                                                                    tAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATTAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACggtg                                           <  8:129823/149‑1 (MQ=255)                                                                                                       aCCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACggtggt                                         >  3:306351/1‑148 (MQ=255)                                                                                                       aCCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCT‑GCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGa                                       <  3:102702/149‑1 (MQ=255)                                                                                                        ccccTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGa                                       <  7:319015/149‑1 (MQ=255)                                                                                                               tACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCtgtg                                >  8:103564/1‑149 (MQ=255)                                                                                                               tACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCATATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCtgtg                                <  7:180624/149‑1 (MQ=255)                                                                                                                aCAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTgg                               >  2:62978/1‑149 (MQ=255)                                                                                                                 cAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCGTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACggtg                                           <  3:308764/136‑1 (MQ=255)                                                                                                                 cAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGCCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGAAGTGAAGACggag                                           >  4:308764/1‑134 (MQ=255)                                                                                                                     ggAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGt                                                                                 >  5:298052/1‑94 (MQ=255)                                                                                                                     ggAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGt                                                                                 <  6:298052/94‑1 (MQ=255)                                                                                                                       aaCGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGTCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCggg                        <  4:202546/149‑1 (MQ=255)                                                                                                                       aaCGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCa                                                                     >  1:160109/1‑104 (MQ=255)                                                                                                                       aaCGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCa                                                                     <  2:160109/104‑1 (MQ=255)                                                                                                                       aaCGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCggg                        <  5:306605/149‑1 (MQ=255)                                                                                                                        aCGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGc                       <  3:241046/149‑1 (MQ=255)                                                                                                                                ccg‑CCACCGCATTCATCAAT‑CGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCcgc               >  2:307541/3‑149 (MQ=255)                                                                                                                                  g‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCtg                                  <  6:270136/128‑1 (MQ=255)                                                                                                                                  g‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCtg                                  >  5:270136/1‑128 (MQ=255)                                                                                                                                     cACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGTTCAATCACCGGGGGGTTGGCTGGCAGGCTCGGCCTAGGGGTGCGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAag           <  5:362167/149‑1 (MQ=255)                                                                                                                                      aCCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGaa            >  5:92662/1‑147 (MQ=255)                                                                                                                                      aCCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGaa            >  6:19409/1‑147 (MQ=255)                                                                                                                                       ccG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAagag         <  5:177902/149‑1 (MQ=255)                                                                                                                                        cG‑ATTGATGAATGCGGGGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCCGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCAAGCCcgctcccaca        <  3:310031/149‑8 (MQ=255)                                                                                                                                           aTTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGTGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGcc      >  2:152543/1‑149 (MQ=255)                                                                                                                                            ttGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGccc     <  7:103564/149‑1 (MQ=255)                                                                                                                                               atgaATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTa  <  8:327168/149‑1 (MQ=255)                                                                                                                                                 gaATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTa  >  8:341571/1‑147 (MQ=255)                                                                                                                                                 |                                                                                                                                                   CTTAGACGCAATGCTGTGAACGGCATGTTCAGAACAGTGACAGGGCGTGGCGGTGGGAGTTTCAGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATG‑CCACCG‑ATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTA  >  NC_002947/2136795‑2137083

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