Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	994,947	T→A	7.3%	intergenic (+26/+99)	PP_0859 → / ← PP_0860	ketoglutaramate omega‑amidase/sulfite reductase flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	994,947	0	T	A	7.3%	196.2 / 4.7	68	intergenic (+26/+99)	PP_0859/PP_0860	ketoglutaramate omega‑amidase/sulfite reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (33/29);  new base A (3/2);  total (36/32)
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.79e-01

CCTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAAC  >  NC_002947/994801‑995085                                                                                                                                                   |                                                                                                                                            ccTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTACAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTg                                                                                                                                           >  1:357190/1‑149 (MQ=255)    gCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCACTCTGGATGCCGATACATTAGAGCTGCACTGACAGCAGGCCCGGCTTCTTCTCttgcttgcac                                                                                                                                        >  1:248956/1‑147 (MQ=255)        aGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCcgcg                                                                                                                                     >  7:138761/1‑148 (MQ=255)         gTGCCGGGGAGGCAGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                   <  6:122531/149‑1 (MQ=255)            ccGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAg                                                                                                                                 <  7:80745/148‑1 (MQ=255)             cGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAggg                                                                                                                               >  1:207235/1‑149 (MQ=255)               gggAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggcc                                                                                                                             >  7:170263/1‑149 (MQ=255)                     cGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGGCCGGCCCTTTCGCGGGCAAGCCCGCGAAGAAGCCGGGCC‑t                                                                                                                      <  4:116005/149‑1 (MQ=37)                       gACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑Tgc                                                                                                                    >  4:92696/1‑149 (MQ=255)                         cGGGGTGTTCACTACAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑Tgcca                                                                                                                  <  2:357190/149‑1 (MQ=255)                                ttCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAAt                                                                                                            <  8:170263/148‑1 (MQ=255)                                   aCTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTcc                                                                                                        <  4:113988/149‑1 (MQ=255)                                       cAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCg                                                                                                    <  1:254249/149‑1 (MQ=255)                                            cTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCg                                                                                                <  5:300126/148‑1 (MQ=255)                                             tGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAGTTTCCACCGGCCGt                                                                                               >  4:309463/1‑148 (MQ=255)                                               aGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCa                                                                                            >  2:123456/1‑149 (MQ=255)                                               aGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCa                                                                                            >  6:327594/1‑149 (MQ=255)                                                  gcggcgGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCg                                                                                         >  7:103620/1‑149 (MQ=255)                                                   cggcggAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGt                                                                                        >  3:307015/1‑149 (MQ=255)                                                    ggcggAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGtt                                                                                       <  2:8205/149‑1 (MQ=255)                                                      cggAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTGCACCGTTgc                                                                                     >  7:302129/1‑149 (MQ=255)                                                        gAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTgccg                                                                                   >  1:53559/1‑149 (MQ=255)                                                         aaCTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTgccgc                                                                                  >  1:357491/1‑149 (MQ=255)                                                           cTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTgccgccg                                                                                <  4:309701/149‑1 (MQ=255)                                                                ggcgTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCAcc                                                                           <  2:53559/149‑1 (MQ=255)                                                                       cGGGCCAAGTTCCCTGCCAATCTGGATGCCGATACATTCGAGCTGAACTGACAGCAGGACCGGCTTCTTCGCGGGCTTGCCCCCGAAAGGGCCGGACC‑TGCCATACAATTTCCACCGGCCGATCCCCGTTGACGCCGCCACCGTAttag                                                                    >  4:55294/1‑147 (MQ=255)                                                                         ggCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGTCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTcgcg                                                                  >  1:263204/1‑149 (MQ=255)                                                                          gCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                 <  3:92696/149‑1 (MQ=255)                                                                          gCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                 <  8:138761/149‑1 (MQ=255)                                                                          gCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                 <  8:103620/149‑1 (MQ=255)                                                                          gCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                 <  1:169633/149‑1 (MQ=255)                                                                          gCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                 <  2:127961/149‑1 (MQ=255)                                                                          gCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                 <  2:357491/149‑1 (MQ=255)                                                                            cAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCggg                                                                >  3:331354/1‑148 (MQ=255)                                                                                ttCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAAc                                                           <  3:93017/149‑1 (MQ=255)                                                                                 gcccGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAAcc                                                          <  8:68044/148‑1 (MQ=255)                                                                                   ccGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCg                                                        <  4:307015/149‑1 (MQ=255)                                                                                       ccAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTcc                                                    <  2:359217/149‑1 (MQ=255)                                                                                              ggATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggccg                                                                                                                            >  6:344521/1‑71 (MQ=255)                                                                                               gaggCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGTCTTCTTCGCGGGCTTGCCCGCGAAAGGTCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTacacac                                             <  2:248956/145‑1 (MQ=255)                                                                                               gATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTacaca                                              <  6:97917/147‑1 (MQ=255)                                                                                                    cGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATc                                        >  8:217351/1‑148 (MQ=255)                                                                                                        aCCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATa                                    >  6:372812/1‑148 (MQ=255)                                                                                                         ccTTCGAGCTGCACTGACAGCAGGTCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATac                                   <  2:263204/148‑1 (MQ=255)                                                                                                          cTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTcgcg                                                                  <  7:234264/116‑1 (MQ=255)                                                                                                          cTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTcgcg                                                                  >  8:234264/1‑116 (MQ=255)                                                                                                          cTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATaca                                  >  7:274446/1‑148 (MQ=255)                                                                                                                   gCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCATGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGcc                         >  1:98164/1‑149 (MQ=255)                                                                                                                     aCTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGGAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGa                       >  7:99081/1‑148 (MQ=255)                                                                                                                      cTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGata                     <  5:372812/149‑1 (MQ=255)                                                                                                                       tGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGcc                         >  3:173845/1‑144 (MQ=255)                                                                                                                       tGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGcc                         <  4:173845/144‑1 (MQ=255)                                                                                                                       tGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGatat                    >  5:281136/1‑149 (MQ=255)                                                                                                                            gcagGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATc                                        <  7:106376/124‑1 (MQ=255)                                                                                                                            gcagGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATc                                        >  8:106376/1‑124 (MQ=255)                                                                                                                                gCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATGAATACACATCTCGCCGATATCGCCCATcc           >  8:89864/1‑149 (MQ=255)                                                                                                                                gCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATcc           >  1:22132/1‑149 (MQ=255)                                                                                                                                gCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATcc           >  6:383523/1‑149 (MQ=255)                                                                                                                                 cccGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCt          <  7:217351/149‑1 (MQ=255)                                                                                                                                 cccGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCt          <  4:331354/149‑1 (MQ=255)                                                                                                                                  ccGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                 <  8:261367/93‑1 (MQ=255)                                                                                                                                   cGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTc         >  8:287550/1‑148 (MQ=255)                                                                                                                                    ggccctttcGCGGGCAAGCCCGCGAAGAAGCCGGGCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGa       >  6:142401/7‑149 (MQ=21)                                                                                                                                       ttcttcGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATc     >  1:25659/1‑148 (MQ=255)                                                                                                                                       ttcttcGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATc     >  4:245249/1‑148 (MQ=255)                                                                                                                                        tcttcGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCaa   >  2:10570/1‑149 (MQ=255)                                                                                                                                        ctttcGCGGGCAAGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCaa   >  4:139888/3‑149 (MQ=255)                                                                                                                                         tttcGCGGGCAAGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAAc  >  4:348419/2‑149 (MQ=255)                                                                                                                                         cttcGCGGGCGTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCaa   <  2:207235/148‑1 (MQ=255)                                                                                                                                                   ttGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCTTACACACCGATCAATACACATCTCGCCGATATCGc                 >  4:318644/1‑124 (MQ=255)                                                                                                                                                   ttGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCTTACACACCGATCAATACACATCTCGCCGATATCGc                 <  3:318644/124‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                            CCTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCC‑TGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAAC  >  NC_002947/994801‑995085

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