Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,402,317	Δ1 bp	7.6%	coding (50/330 nt)	PP_2103 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,402,317	0	T	.	7.6%	206.7 / 10.7	66	coding (50/330 nt)	PP_2103	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/29);  new base . (2/3);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

ATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGGCGCACGCCGCAGTA  >  NC_002947/2402170‑2402464                                                                                                                                                     |                                                                                                                                                    aTGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCtt                                                                                                                                                    >  8:199721/1‑149 (MQ=255)    ccccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgcg                                                                                                                                                 <  8:142971/149‑1 (MQ=255)     cccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgcgc                                                                                                                                                <  2:234554/149‑1 (MQ=255)     cccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgcg                                                                                                                                                 >  4:179825/1‑148 (MQ=255)      ccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgcg                                                                                                                                                 >  7:44709/1‑147 (MQ=255)      ccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCa                                                                                                                                               <  3:53024/149‑1 (MQ=255)        acaGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGCAAGGCGCTGCTGCGCC‑TGCGCatg                                                                                                                                             <  1:182297/149‑3 (MQ=255)         caGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGTTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGCAAGGCGCTGCTGCGCC‑TGCGCatgg                                                                                                                                            <  5:82257/149‑4 (MQ=255)              gcgcATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTcgcgc                                                                                                                                       >  3:248767/1‑149 (MQ=255)                gcATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGCCGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCgg                                                                                                                                     <  5:221739/149‑1 (MQ=255)                gcATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCgg                                                                                                                                     >  5:158045/1‑149 (MQ=255)                     cTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTAc                                                                                                                                <  4:220838/149‑1 (MQ=255)                        gATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTAcgcg                                                                                                                             >  4:119621/1‑149 (MQ=255)                             ccAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTcgcgc                                                                                                                                       <  4:248368/134‑1 (MQ=255)                             ccAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTcgcgc                                                                                                                                       >  3:248368/1‑134 (MQ=255)                             ccAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGttt                                                                                                                        >  6:186791/1‑149 (MQ=255)                             ccAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGttt                                                                                                                        >  3:206914/1‑149 (MQ=255)                               aGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTgg                                                                                                                      <  3:119621/149‑1 (MQ=255)                               aGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTgg                                                                                                                      <  4:248767/149‑1 (MQ=255)                                ggCCTTCATGTCTGGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGc                                                                                                                     <  1:132381/149‑1 (MQ=255)                                 gCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCa                                                                                                                    <  4:47922/149‑1 (MQ=255)                                   cTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGt                                                                                                                  >  3:47707/1‑149 (MQ=255)                                    ttCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGtt                                                                                                                 >  7:175271/1‑149 (MQ=255)                                      cATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGtt                                                                                                                 <  2:218036/147‑1 (MQ=255)                                        tGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTgg                                                                                                             <  2:77729/149‑1 (MQ=255)                                            gCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAgg                                                                                                         <  3:181223/149‑1 (MQ=255)                                                  tGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATg                                                                                                   >  6:276826/1‑149 (MQ=255)                                                   gCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATgg                                                                                                  >  1:120325/1‑149 (MQ=255)                                                    cccATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATgg                                                                                                  <  8:175271/148‑1 (MQ=255)                                                    cccATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATgg                                                                                                  >  2:315133/1‑148 (MQ=255)                                                       aTACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGACaa                                                                                               >  1:41779/1‑148 (MQ=255)                                                         aCCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgc                                                                                                                                                  >  5:322279/1‑95 (MQ=255)                                                         aCCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgc                                                                                                                                                  <  6:322279/95‑1 (MQ=255)                                                           cgcgCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGgc                                                                                          >  7:200695/1‑149 (MQ=255)                                                           cgcgCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAAGCTCATGGCAACAGgc                                                                                          >  5:360047/1‑149 (MQ=255)                                                                aCGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTc                                                                                      >  1:140619/1‑148 (MQ=255)                                                                aCGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTc                                                                                      <  2:140619/148‑1 (MQ=255)                                                                         aGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTg                                                                            <  4:47707/149‑1 (MQ=255)                                                                           cGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGc                                                                          <  3:12254/149‑1 (MQ=255)                                                                             gTTCAAGCAACTGACCGGACTCATGACGGATTTCCTTGCGGTGCATATACATGCGC‑AGGCGCAGCATCGCCTAGCGAAGTTCCCGCGGGTTAAGCGTGTTTGGCAGTTCTTTCAGGCTCATGGCAACAGGCGCTCCCTGACCTTAGCac                                                                        >  8:68349/1‑148 (MQ=255)                                                                                 aaGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTtc                                                                    >  6:174083/1‑149 (MQ=255)                                                                                   gCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTtcgt                                                                  <  4:206914/149‑1 (MQ=255)                                                                                            cGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGa                                                               <  8:154783/143‑1 (MQ=255)                                                                                            cGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGa                                                               >  7:154783/1‑143 (MQ=255)                                                                                             ggACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGgc                                                        >  3:277824/1‑149 (MQ=255)                                                                                                    tGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCATTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGACGaat                                                 >  7:175993/1‑148 (MQ=255)                                                                                                     gACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAgg                                                >  5:147461/1‑149 (MQ=255)                                                                                                     gACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAgg                                                >  5:335264/1‑149 (MQ=255)                                                                                                      aCGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAggg                                               <  1:315133/149‑1 (MQ=255)                                                                                                       cGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATAGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGc                                              <  6:158045/149‑1 (MQ=255)                                                                                                        ggATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGc                                              >  1:53770/1‑148 (MQ=255)                                                                                                        ggATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGc                                              >  7:183713/1‑148 (MQ=255)                                                                                                           tttCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCgat                                           >  2:39271/1‑148 (MQ=255)                                                                                                               cTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGGTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATtc                                      >  3:190336/1‑149 (MQ=255)                                                                                                                     gTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGaa                                <  2:120325/149‑1 (MQ=255)                                                                                                                      tGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGaa                                >  8:188895/1‑148 (MQ=255)                                                                                                                      tGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAac                               <  4:133967/149‑1 (MQ=255)                                                                                                                      tGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAac                               <  7:199721/149‑1 (MQ=255)                                                                                                                                tGCGCAAGGCGCTGCTGCGCC‑TGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCt                     >  4:316955/1‑149 (MQ=255)                                                                                                                                   gcaaggcgctgCTGCGCC‑TGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCtt                    >  4:218563/4‑147 (MQ=255)                                                                                                                                   gc‑aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCtt                    >  3:125791/1‑147 (MQ=255)                                                                                                                                      aggcgctgCTGCGCC‑TGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGgc               <  3:218563/149‑1 (MQ=255)                                                                                                                                             gcagcGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGGCGCAcgc         <  5:174083/149‑1 (MQ=255)                                                                                                                                               agcGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGGCGCAcgccg       <  7:68349/149‑1 (MQ=255)                                                                                                                                               agcGCCTTGCGCAGTTCGCGCGGGTTACGCGTGATTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGGCGCAcgcag       >  1:322867/1‑147 (MQ=255)                                                                                                                                                 cgcCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGGCGCACGCCGCa     <  5:186791/149‑1 (MQ=255)                                                                                                                                                    cTTGCGCAGTTCGCGCGGGTTACGGGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGGCGCACGCCGCAGTa  <  2:322867/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    ATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGGCGCACGCCGCAGTA  >  NC_002947/2402170‑2402464

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