Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,950,194:1	+C	10.7%	coding (150/969 nt)	fliM ←	flagellar biosynthesis switching/energizing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,950,194	1	.	C	10.7%	265.4 / 29.0	75	coding (150/969 nt)	fliM	flagellar biosynthesis switching/energizing protein
Reads supporting (aligned to +/- strand):  ref base . (32/35);  new base C (4/4);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01

GTACAGCGAGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGACAGCAGGTCCTGTACGG  >  NC_002947/4950053‑4950339                                                                                                                                               |                                                                                                                                                  gTACAGCGAGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTC‑GCATGCCGAccct                                                                                                                                            <  6:244068/149‑2 (MQ=255)   aCAGCGAGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATg                                                                                                                                                    <  3:300115/140‑1 (MQ=255)   aCAGCGAGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATg                                                                                                                                                    >  4:300115/1‑140 (MQ=255)     aGCGAGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCAc                                                                                                                                        <  7:235957/149‑1 (MQ=255)     aGCGAGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTC‑GCATGCCGAccc                                                                                                                                             <  4:209177/144‑1 (MQ=255)     aGCGAGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTC‑GCATGCCGAccc                                                                                                                                             >  3:209177/1‑144 (MQ=255)      gCGAGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACg                                                                                                                                       >  4:197242/1‑149 (MQ=255)         aGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTC‑GCATGCCGAccct                                                                                                                                            <  2:71672/141‑2 (MQ=255)         aGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTC‑GCATGCCGAccct                                                                                                                                            >  1:71672/1‑140 (MQ=255)           tGCACGTATTCGTCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACg                                                                                                                                  <  7:230690/149‑1 (MQ=255)             cACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGt                                                                                                                                >  7:152286/1‑149 (MQ=255)             cACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGt                                                                                                                                >  4:222657/1‑149 (MQ=255)              aCGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTc                                                                                                                               >  3:15773/1‑149 (MQ=255)              aCGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTc                                                                                                                               >  4:67185/1‑149 (MQ=255)                gTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCt                                                                                                                             <  3:222657/149‑1 (MQ=255)                  aTTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGa                                                                                                                           >  4:52020/1‑149 (MQ=255)                  aTTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTAGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGa                                                                                                                           >  5:29158/1‑149 (MQ=255)                   ttCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGa                                                                                                                           <  3:67185/148‑1 (MQ=255)                      gCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACgg                                                                                                                                 >  4:211691/1‑139 (MQ=255)                      gCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACgg                                                                                                                                 <  3:211691/139‑1 (MQ=255)                      gCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTgg                                                                                                                       >  1:225341/1‑149 (MQ=255)                       ccGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGt                                                                                                                      >  4:226514/1‑149 (MQ=255)                       cagAACTTCATGACCTGCACGCCACCGACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGATGATACGGGTATAGCGGGCGAAGCGCCCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGt                                                                                                                      <  1:250812/147‑1 (MQ=255)                         gAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCa                                                                                                                    <  3:197242/149‑1 (MQ=255)                           aCTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAgg                                                                                                                  <  3:226514/149‑1 (MQ=255)                            cTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGa                                                                                                                           <  4:241280/139‑1 (MQ=255)                            cTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGa                                                                                                                           >  3:241280/1‑139 (MQ=255)                                aTGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCg                                                                                                               >  6:46229/1‑147 (MQ=255)                                 tGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACCAGTCAggaccgta                                                                                                            <  6:176014/149‑7 (MQ=255)                                     cTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCttt                                                                                                        >  5:38945/1‑149 (MQ=255)                                     cTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCATTGATCATTTCCAGGGTC‑GCATGCCGAccct                                                                                                                                            >  3:289149/1‑112 (MQ=255)                                     cTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCATTGATCATTTCCAGGGTC‑GCATGCCGAccct                                                                                                                                            <  4:289149/113‑2 (MQ=255)                                      tGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCtttt                                                                                                       >  3:57463/1‑149 (MQ=255)                                       gCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCtttt                                                                                                       >  6:3805/1‑148 (MQ=255)                                        cACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCtttt                                                                                                       <  5:112289/147‑1 (MQ=255)                                        cACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGa                                                                                                     <  2:225341/149‑1 (MQ=255)                                           gCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCt                                                                                                                             >  4:102076/1‑122 (MQ=255)                                           gCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCt                                                                                                                             <  3:102076/122‑1 (MQ=255)                                            ccacccacGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGGCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGGCAGGTCGTAGCTTTTGATACt                                                                                                 <  3:220287/149‑1 (MQ=255)                                             cacccacGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTg                                                                                                >  5:256324/1‑149 (MQ=255)                                                     gCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTc                                                                                        >  7:57265/1‑149 (MQ=255)                                                      ccACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCg                                                                                       >  3:51537/1‑149 (MQ=255)                                                             agcggagcggCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGc                                                                                  >  5:240944/1‑147 (MQ=255)                                                             agcggagcggCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGc                                                                                  <  6:240944/147‑1 (MQ=255)                                                               cggagcggCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGCCCCCGCACGAGACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCAAGGCTCGGATGCActc                                                                              >  5:292143/1‑149 (MQ=255)                                                               cggagcggCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTg                                                                                                <  5:228995/131‑1 (MQ=255)                                                               cggagcggCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTg                                                                                                >  6:228995/1‑131 (MQ=255)                                                                           gcagGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCAc                                                                   >  7:194872/1‑148 (MQ=255)                                                                            cagGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCa                                                                 <  2:196728/149‑1 (MQ=255)                                                                                  gAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCcg                                                            <  8:57265/148‑1 (MQ=255)                                                                                      aTGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATTCTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCgtcgtc                                                       <  4:57463/149‑1 (MQ=255)                                                                                       tGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGtcgtcg                                                      <  6:38945/149‑1 (MQ=255)                                                                                        gCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGa                                                     >  8:84351/1‑149 (MQ=255)                                                                                                gggTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGc                                             <  7:84351/149‑1 (MQ=255)                                                                                                       gCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAggg                                      <  4:15773/149‑1 (MQ=255)                                                                                                        cGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAggg                                      >  6:166073/1‑148 (MQ=255)                                                                                                        cGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAgg                                       <  8:247993/147‑1 (MQ=255)                                                                                                        cGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAgg                                       >  7:247993/1‑147 (MQ=255)                                                                                                          ggCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCAt                                   <  5:3805/149‑1 (MQ=255)                                                                                                             gAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATg                                                                                   >  8:128440/1‑98 (MQ=255)                                                                                                             gAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATg                                                                                   <  7:128440/98‑1 (MQ=255)                                                                                                               aGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGAtt                              <  3:52020/149‑1 (MQ=255)                                                                                                                       ttGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTg                      <  4:51537/149‑1 (MQ=255)                                                                                                                         gATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTgg                     <  5:46229/148‑1 (MQ=255)                                                                                                                         gATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTgg                     >  4:109252/1‑148 (MQ=255)                                                                                                                         gATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTgg                     <  6:256324/148‑1 (MQ=255)                                                                                                                           tCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGa                   >  5:20373/1‑148 (MQ=255)                                                                                                                           tCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGa                   >  1:209210/1‑148 (MQ=255)                                                                                                                           tCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGAc                  <  8:194872/149‑1 (MQ=255)                                                                                                                           tCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGCCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGa                   >  2:264987/1‑148 (MQ=255)                                                                                                                               ttCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGAcag                >  5:282796/1‑147 (MQ=255)                                                                                                                                 ccATGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCAc                                                                   <  8:210627/94‑1 (MQ=255)                                                                                                                                   ggggTC‑GCATGCCGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGACAGCAGGTc          >  4:63014/2‑149 (MQ=255)                                                                                                                                   aGGTTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGCCGCCATGCAACAGGGCATCGATTTCATCCTGGGACAGCAGGTc          <  6:292143/149‑1 (MQ=255)                                                                                                                                    gggTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCGACAGGGCATCGATTTCATCCTGGGACAGCAGGTcc         <  8:152286/149‑1 (MQ=255)                                                                                                                                          gCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGACAGCAGGTCCTGTACg   <  2:209210/149‑1 (MQ=255)                                                                                                                                           cATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGACAGCAGGTCCTGTACgg  <  6:20373/149‑1 (MQ=255)                                                                                                                                              g‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGACAGCAGGTCCTGTa     <  4:89248/143‑1 (MQ=255)                                                                                                                                              g‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGACAGCAGGTCCTGTa     >  3:89248/1‑143 (MQ=255)                                                                                                                                               |                                                                                                                                                  GTACAGCGAGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGACAGCAGGTCCTGTACGG  >  NC_002947/4950053‑4950339

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