Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,787,767	G→A	59.7%	Q290* (CAG→TAG)	PP_3350 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,787,767	0	G	A	59.7%	40.2 / 77.0	72	Q290* (CAG→TAG)	PP_3350	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (13/16);  new base A (19/24);  total (32/40)
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.23e-01

GACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTTCCA  >  NC_002947/3787621‑3787911                                                                                                                                                   |                                                                                                                                                 gACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTAcc                                                                                                                                                <  1:118688/149‑1 (MQ=255) gACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTAcc                                                                                                                                                <  4:157425/149‑1 (MQ=255)  aCTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTAccc                                                                                                                                               >  1:249672/1‑149 (MQ=255)        ccGTGGCCCAGGTCACTGTAGCCGCTTGCGGCATGGATGCGCATGGACAGGCGCGGTTGGCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGAGTATCGACGCGTTGCCCAGGTGGCCGGTCTGCCcgacga                                                                                                                                          <  8:61915/148‑1 (MQ=255)        ccGTCGCCCAGGTCGCTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCcgacga                                                                                                                                          <  4:200567/148‑1 (MQ=255)            cGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTc                                                                                                                                     <  8:174891/149‑1 (MQ=255)               ccAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAg                                                                                                                                  >  6:219484/1‑149 (MQ=255)                cAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGc                                                                                                                                 <  6:167723/149‑1 (MQ=255)                 aGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCt                                                                                                                                >  6:243010/1‑149 (MQ=255)                  ggTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAg                                                                                                                                         <  4:302154/139‑1 (MQ=255)                  ggTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAg                                                                                                                                         >  3:302154/1‑139 (MQ=255)                   gTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTcc                                                                                                                                    >  8:87724/1‑143 (MQ=255)                   gTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTcc                                                                                                                                    <  7:87724/143‑1 (MQ=255)                    tCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCTACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCtcca                                                                                                                             <  2:249672/149‑1 (MQ=255)                     cACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCtccag                                                                                                                            >  7:342569/1‑149 (MQ=255)                          tCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCtccagtccag                                                                                                                       >  6:189094/1‑149 (MQ=255)                              cgcTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGcc                                                                                                                   <  3:336247/149‑1 (MQ=255)                                   gCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGAt                                                                                                               <  2:215655/148‑1 (MQ=255)                                   gCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGa                                                                                                                <  1:169830/147‑1 (MQ=255)                                   gCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATg                                                                                                              >  8:120826/1‑149 (MQ=255)                                      gCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATggcg                                                                                                           >  5:233756/1‑149 (MQ=255)                                            aTGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGc                                                                                                                    <  1:136146/134‑1 (MQ=255)                                            aTGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGc                                                                                                                    >  2:136146/1‑134 (MQ=255)                                            aTGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGa                                                                                                      >  6:114449/1‑148 (MQ=255)                                                 cATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGcc                                                                                                <  5:219484/149‑1 (MQ=255)                                                  aTGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCcg                                                                                               >  2:259684/1‑149 (MQ=255)                                                    ggCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCcgc                                                                                              <  8:342569/148‑1 (MQ=255)                                                     gCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGAACGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGt                                                                                            >  6:346102/1‑149 (MQ=255)                                                      ccAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCtc                                                                                                                               <  6:31620/113‑1 (MQ=255)                                                      ccAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCtc                                                                                                                               >  5:31620/1‑113 (MQ=255)                                                          gcgcGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCg                                                                                        >  6:36033/1‑148 (MQ=255)                                                            gcgGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGaaa                                                                                     <  1:259684/149‑1 (MQ=255)                                                                   ccAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAg                                                                               <  6:201652/148‑1 (MQ=255)                                                                   ccAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAg                                                                               <  8:311003/148‑1 (MQ=255)                                                                    cAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAgcg                                                                             >  6:101710/1‑149 (MQ=255)                                                                     aGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAgcgc                                                                            >  6:307824/1‑149 (MQ=255)                                                                        ggtgggtgAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCgctgc                                                                         <  6:233756/149‑1 (MQ=255)                                                                          tgggtgAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTc                                                                       <  5:101710/149‑1 (MQ=255)                                                                               gaaggtgaagCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCAc                                                                  >  7:177522/1‑149 (MQ=255)                                                                               gaaggtgaagCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCAc                                                                  >  8:193513/1‑149 (MQ=255)                                                                                 aggtgaagCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCAccc                                                                <  3:21215/149‑1 (MQ=255)                                                                                 aggtgaagCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCAccc                                                                >  2:599/1‑149 (MQ=255)                                                                                  ggtgaagCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCg                                                               >  7:320092/1‑149 (MQ=255)                                                                                   gtgaagCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGa                                                              >  3:181369/1‑149 (MQ=255)                                                                                        ggcGCTGTCGAGCGACAGCGGCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGCCGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGc                                                         >  3:133333/3‑149 (MQ=255)                                                                                         gcGCTGTCGAGCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGcc                                                        >  3:159172/2‑149 (MQ=255)                                                                                                 gCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAgg                                                 <  5:243010/148‑1 (MQ=255)                                                                                                   tCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTg                                              >  1:196862/1‑149 (MQ=255)                                                                                                   tCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTg                                              >  5:225908/1‑149 (MQ=255)                                                                                                      aCAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCgg                                           <  5:307824/149‑1 (MQ=255)                                                                                                      aCAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCgg                                           <  4:159172/149‑1 (MQ=255)                                                                                                       cAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGgt                                          >  7:48082/1‑149 (MQ=255)                                                                                                       cAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGgt                                          >  6:21391/1‑149 (MQ=255)                                                                                                          cgcCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGgtcgt                                       <  7:120826/149‑1 (MQ=255)                                                                                                          cgcCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGgtcgt                                       <  7:193513/149‑1 (MQ=255)                                                                                                           gcCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGGCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTc                                                                    <  4:152371/119‑1 (MQ=255)                                                                                                           gcCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTc                                                                    >  3:152371/1‑119 (MQ=255)                                                                                                            gccAGGCGCGGATCGACGCGTTGCCCAGGGGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGtcgtcg                                     <  4:133333/148‑1 (MQ=255)                                                                                                            cccAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGtcgtcg                                     <  8:320092/149‑1 (MQ=255)                                                                                                            cccAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGtcgtcg                                     <  5:189094/149‑1 (MQ=255)                                                                                                                     ggATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAgg                            <  5:21391/149‑1 (MQ=255)                                                                                                                            cgcgTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAgg                     >  3:288778/1‑149 (MQ=255)                                                                                                                                ttGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAg                 <  8:177522/149‑1 (MQ=255)                                                                                                                                   cccAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAgg                     <  5:107839/142‑1 (MQ=255)                                                                                                                                   cccAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAgg                     >  6:107839/1‑142 (MQ=255)                                                                                                                                      aGGTGGCCGGTCTGCCCGACGAGGATCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTc           <  8:48082/149‑1 (MQ=255)                                                                                                                                        gTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGtcgtcg                                     >  5:275038/1‑121 (MQ=255)                                                                                                                                        gTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGtcgtcg                                     <  6:275038/121‑1 (MQ=255)                                                                                                                                         tGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCAt        <  5:114449/149‑1 (MQ=255)                                                                                                                                             cGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTTc    <  6:225908/149‑1 (MQ=255)                                                                                                                                               gTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTTCCa  <  5:346102/149‑1 (MQ=255)                                                                                                                                               gTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTTCCa  <  4:181369/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTGCGGTCGTCGGTTTCCAGGCCAAAGGCTAGCAGGTCCATGTTCCA  >  NC_002947/3787621‑3787911

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