Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,631,094	T→G	10.5%	intergenic (+28/‑153)	PP_1430 → / → lepA	DegP‑like serine endoprotease/elongation factor 4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,631,094	0	T	G	10.5%	187.4 / 13.9	67	intergenic (+28/‑153)	PP_1430/lepA	DegP‑like serine endoprotease/elongation factor 4
Reads supporting (aligned to +/- strand):  ref base T (29/31);  new base G (4/3);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

CATCGGTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGAGT  >  NC_002947/1630946‑1631238                                                                                                                                                     |                                                                                                                                                  cATCGGTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCtt                                                                                                                                                   >  8:89264/1‑149 (MQ=255) cATCGGTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCtt                                                                                                                                                   >  7:9345/1‑149 (MQ=255)      gTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCt                                                                                                                                              <  6:34857/149‑1 (MQ=255)       ttCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTg                                                                                                                                             <  2:73870/149‑1 (MQ=255)       atccgCCAAGGAATTCACCGAAATCGCCAAGACACTGACGAAGAATCGTTCCGTGTCGATCCGCATACTGAGTCAGGGCCGTTCAAGCTTCATCACTTTCAAGCTAGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTg                                                                                                                                             >  6:20277/5‑149 (MQ=255)             cAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCttt                                                                                                                                       <  3:86038/149‑1 (MQ=255)             cAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCttt                                                                                                                                       >  6:299060/1‑149 (MQ=255)             cAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCttt                                                                                                                                       >  7:87790/1‑149 (MQ=255)             cAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCttt                                                                                                                                       >  2:237239/1‑149 (MQ=255)              aaGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCCGAAGCTGCCCttt                                                                                                                                       <  1:120859/148‑1 (MQ=255)              aaGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCCGAAGCTGCCCttt                                                                                                                                       >  2:120859/1‑148 (MQ=255)               aGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCttttt                                                                                                                                     <  7:89264/149‑1 (MQ=255)                ggAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTc                                                                                                                                    <  2:103801/149‑1 (MQ=255)                  aaTTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCAAAAGGGCAGCCGAAGCTGCCCtttcgcct                                                                                                                                  <  5:123714/149‑6 (MQ=255)                   aTTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATg                                                                                                                                 <  6:10615/149‑1 (MQ=255)                    ttCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGa                                                                                                                                <  8:87790/149‑1 (MQ=255)                    ttCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCCGAAGCTGCCCtttcgccttg                                                                                                                                <  2:175759/149‑8 (MQ=255)                       aCCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATc                                                                                                                             <  5:299060/149‑1 (MQ=255)                         cGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCtt                                                                                                                           >  6:266993/1‑149 (MQ=255)                                  cAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑c                                                                                                                 >  2:18371/1‑149 (MQ=255)                                   aaGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAAtt                                                                                                                    >  8:123753/1‑146 (MQ=255)                                   aaGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAAtt                                                                                                                    <  7:123753/146‑1 (MQ=255)                                      gAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTca                                                                                                                         >  3:162990/1‑138 (MQ=255)                                      gAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTca                                                                                                                         <  4:162990/138‑1 (MQ=255)                                            acGAAGAATCGTTCCGTGTCGATGCTCGTACTGCGTCAGGTCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTAGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGAt                                                                                                        <  8:9345/147‑1 (MQ=255)                                               aagaaTCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCt                                                                                                    >  5:198468/1‑149 (MQ=255)                                                   aTCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTACATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTATTTCATGAATCTTCACAATTG‑CCGGAATGATTGCtgtg                                                                                                 >  6:135782/1‑148 (MQ=255)                                                    tCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAAATG‑CCGGAATGATTGCTgtgt                                                                                                >  7:220399/1‑148 (MQ=255)                                                    tCGTTCCGGGTCGATGCGCGTACTGCGTCAGGGCCGTGCCATCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTa                                                                                               <  7:294780/149‑1 (MQ=255)                                                           gtgtCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCAAAATTG‑CAGGAATGATTGCTGTGTACATTCg                                                                                         >  6:316520/1‑148 (MQ=255)                                                           gtgtCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCAAAATTG‑ACTGAATGATTGCTGTGTACAGTCg                                                                                         >  6:311947/1‑148 (MQ=255)                                                            tgtCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGa                                                                                       <  4:33368/149‑1 (MQ=255)                                                                   gcgcGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAgg                                                                                 >  4:298715/1‑148 (MQ=255)                                                                   gcgcGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAgg                                                                                 >  5:77818/1‑148 (MQ=255)                                                                        tACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAAtctc                                                                           <  6:77818/149‑1 (MQ=255)                                                                           tGCGTCAGGGCCGTGCCAGCTTCATCACTTTTAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCa                                                                        >  7:77525/1‑149 (MQ=255)                                                                           tGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCa                                                                        <  3:298715/149‑1 (MQ=255)                                                                             cGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGCCAGAGGAATCTCCCAt                                                                       >  6:101769/1‑148 (MQ=255)                                                                               tCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAAc                                                                    >  4:163943/1‑149 (MQ=255)                                                                                   ggCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCca                                                                 <  1:237239/149‑1 (MQ=255)                                                                                           cAGCTTCATCACTTTCAATCTTGATGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCACTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTaa                                                        <  5:316520/149‑1 (MQ=255)                                                                                           cAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTaa                                                        <  2:100701/149‑1 (MQ=255)                                                                                           cAGCTTAATCACTTTCAAGCTTTCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCACTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTaa                                                        <  5:311947/149‑1 (MQ=255)                                                                                                   tcaCTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATGTCACAAAACCCAGCAGCTTAAGGTACAAt                                                >  3:271324/1‑149 (MQ=255)                                                                                                          cAAGCTTGATGAATAAGCAGTTTCGCTAGGCTAAAGGGCAGCTTCGGCTGCCCTTTTTCATTAATGTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGTCAGAGGAATCTCCCATAACCCAGCAGGTTAAGGTACAATTCCCgg                                          <  5:20277/148‑1 (MQ=255)                                                                                                           aaGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCt                                        >  6:56971/1‑149 (MQ=255)                                                                                                             gCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTAt                                      >  7:232604/1‑149 (MQ=255)                                                                                                                    gAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑GCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTcgg                               >  8:204968/1‑149 (MQ=255)                                                                                                                    gAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTcgg                               >  3:79413/1‑149 (MQ=255)                                                                                                                     aaTAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTcg                                <  6:198468/147‑1 (MQ=255)                                                                                                                     aaTAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTcg                                <  7:99915/147‑1 (MQ=255)                                                                                                                     aaTAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTATTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTcg                                <  5:135782/147‑1 (MQ=255)                                                                                                                          gCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCggcgggcg                         >  4:194697/1‑149 (MQ=255)                                                                                                                           cAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTcggcg                             >  4:902/1‑144 (MQ=255)                                                                                                                           cAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTcggcg                             <  3:902/144‑1 (MQ=255)                                                                                                                           cAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGCCTATTTTTCGGCGGGCGt                        >  4:37581/1‑149 (MQ=255)                                                                                                                            aGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTc                       <  5:56971/149‑1 (MQ=255)                                                                                                                                   cAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCg                >  5:39389/1‑149 (MQ=255)                                                                                                                                      tgaaaaaGGGCAGCCGAAGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAg             >  7:126876/5‑149 (MQ=255)                                                                                                                                       gCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGc            <  4:79413/149‑1 (MQ=255)                                                                                                                                        aaaaaGGGCAGCCGAAGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGcc           >  3:177990/3‑149 (MQ=255)                                                                                                                                         gAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTc                       >  3:293553/1‑136 (MQ=255)                                                                                                                                         gAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTc                       <  4:293553/136‑1 (MQ=255)                                                                                                                                         gAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCt          <  3:37581/149‑1 (MQ=255)                                                                                                                                         aaaaGGGCAGCCGAAGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCt          >  6:135376/2‑149 (MQ=255)                                                                                                                                             gggCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGt                        <  3:148377/131‑1 (MQ=255)                                                                                                                                             gggCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGt                        >  4:148377/1‑131 (MQ=255)                                                                                                                                                cAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGAg   <  8:232604/149‑1 (MQ=255)                                                                                                                                                 aGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGAgt  <  3:194697/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CATCGGTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTG‑CCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGAGT  >  NC_002947/1630946‑1631238

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