Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,867,738	Δ1 bp	7.5%	intergenic (‑33/+19)	thyA ← / ← lgt	thymidylate synthase/prolipoprotein diacylglyceryl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,867,738	0	G	.	7.5%	207.7 / 10.6	67	intergenic (‑33/+19)	thyA/lgt	thymidylate synthase/prolipoprotein diacylglyceryl transferase
Reads supporting (aligned to +/- strand):  ref base G (28/33);  new base . (3/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.23e-01

GCAGGTCGAAACGCAGCATCGCGCCCGGCAGG‑CTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTTGT  >  NC_002947/5867591‑5867883                                                                                                                                                      |                                                                                                                                                  gCAGGTCGAAACGCAGCATCGCGCCCGGCAGGACTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCg                                                                                                                                                   >  1:40280/1‑149 (MQ=255)  cAGGTCGAAACGCAGCATCGCGCCCGGCAGG‑CTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGc                                                                                                                                                  <  6:120839/148‑1 (MQ=255)      tCGAAACGCAGCATCGCGCCCGGCAGG‑CTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGCGATC‑CCCCTg                                                                                                                                             <  1:188204/149‑1 (MQ=255)       cGAAACGCAGCATCGCGCCCGGCAGG‑CTGATGTTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGc                                                                                                                                            >  6:92436/1‑149 (MQ=255)           aCGCAGCATCGCGCCCGGCAGG‑CTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTg                                                                                                                                        <  8:300103/149‑1 (MQ=255)            cGCAGCATCGCGCCCGGCAGG‑CTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTgg                                                                                                                                       <  7:200004/149‑1 (MQ=255)                 caTCGCGCCCGGCAGG‑CTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGAtt                                                                                                                                   >  6:199172/1‑148 (MQ=255)                      cgcCCGGCAGG‑CTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAg                                                                                                                             >  1:326518/1‑149 (MQ=255)                        cccGGCAGG‑CTGATGGGGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGc                                                                                                                           <  4:191260/149‑1 (MQ=255)                         ccGGCAGG‑CTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCg                                                                                                                          <  2:68717/149‑1 (MQ=255)                             cAGG‑CTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGtt                                                                                                                      >  1:240976/1‑149 (MQ=255)                                        gTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAgg                                                                                                            <  7:321598/149‑1 (MQ=255)                                        gTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAg                                                                                                             <  2:40280/148‑1 (MQ=255)                                         tGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGc                                                                                                           >  8:122042/1‑149 (MQ=255)                                                     tGCGGTTTCCCTGCAACGTGCCGTTATCTATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATaa                                                                                               <  3:114271/149‑1 (MQ=255)                                                     tGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATaa                                                                                               <  7:122042/149‑1 (MQ=255)                                                      gCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAg                                                                                              <  5:199172/149‑1 (MQ=255)                                                       cGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGGAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAg                                                                                              <  5:92436/148‑1 (MQ=255)                                                         gTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGcc                                                                                            >  3:295486/1‑148 (MQ=255)                                                                  cAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATc                                                                                   <  4:332752/148‑1 (MQ=255)                                                                         ccGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTTTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAg                                                                           <  2:354633/149‑1 (MQ=255)                                                                          cGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAgg                                                                          <  7:108922/149‑1 (MQ=255)                                                                            ttATCAATGACGTCGCGGACAAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATAAAGACGAGGcc                                                                        <  2:326518/149‑1 (MQ=255)                                                                                cAATGACGCCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCAcc                                                                     >  2:77913/1‑148 (MQ=255)                                                                                 aaTGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGa                                                                   >  3:117488/1‑149 (MQ=255)                                                                                  aTGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATaa                                                                                               <  2:148871/120‑1 (MQ=255)                                                                                  aTGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATaa                                                                                               >  1:148871/1‑120 (MQ=255)                                                                                               ccAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATAATCGGAATg                                                     <  1:77913/149‑1 (MQ=255)                                                                                                cAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAg                                                                           >  6:174048/1‑126 (MQ=255)                                                                                                cAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAg                                                                           <  5:174048/126‑1 (MQ=255)                                                                                                  ggTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGca                                                   >  6:344461/1‑148 (MQ=255)                                                                                                   gTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGca                                                   >  4:71882/1‑147 (MQ=255)                                                                                                    tCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCgg                                                         >  6:361476/1‑140 (MQ=255)                                                                                                    tCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCgg                                                         <  5:361476/140‑1 (MQ=255)                                                                                                       aGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCa                                              >  8:208972/1‑148 (MQ=255)                                                                                                        gATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAAc                                            >  4:237752/1‑149 (MQ=255)                                                                                                          tACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGca                                                   <  8:98854/140‑1 (MQ=255)                                                                                                          tACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGca                                                   >  7:98854/1‑140 (MQ=255)                                                                                                           aCTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTgg                                         >  4:369495/1‑149 (MQ=255)                                                                                                                ttCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATaa                                                                                               >  4:205820/1‑90 (MQ=255)                                                                                                                ttCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATaa                                                                                               <  3:205820/90‑1 (MQ=255)                                                                                                                     ggCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtgagt                                 >  1:27215/1‑147 (MQ=255)                                                                                                                      gCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtgagtga                               <  2:27215/148‑1 (MQ=255)                                                                                                                         tACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGagc                             >  3:274385/1‑147 (MQ=255)                                                                                                                            cTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGagccagcc                        >  5:41455/1‑149 (MQ=255)                                                                                                                              gTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAg                                                                           <  5:346574/96‑1 (MQ=255)                                                                                                                              gTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAg                                                                           >  6:346574/1‑96 (MQ=255)                                                                                                                              gTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCg                       >  4:42433/1‑148 (MQ=255)                                                                                                                                 gCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACGAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAAc                    <  3:42433/148‑1 (MQ=255)                                                                                                                                 gCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAAcg                   <  3:71882/149‑1 (MQ=255)                                                                                                                                 gCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAAc                    <  3:237752/148‑1 (MQ=255)                                                                                                                                   aaaaaCGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAAcgc                  >  5:242860/1‑148 (MQ=255)                                                                                                                                    aaaaCGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGa                >  6:24391/1‑149 (MQ=255)                                                                                                                                     ccaCGGCAGGGGCGATC‑CCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtgagtga                               <  6:387521/131‑1 (MQ=255)                                                                                                                                     ccaCGGCAGGGGCGATC‑CCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtgagtga                               >  5:387521/3‑133 (MQ=255)                                                                                                                                     aaaCGGCAGGGG‑GATCTCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGAt               <  7:106550/149‑1 (MQ=255)                                                                                                                                      caCGGCAGGGGCGATC‑CCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATGAAGACGATGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATg              >  4:278/2‑149 (MQ=255)                                                                                                                                      aaCGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATg              >  2:179205/1‑149 (MQ=255)                                                                                                                                         ggCAGGGGCGATC‑CCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTaa           >  3:288881/1‑149 (MQ=255)                                                                                                                                         ggCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTaa           <  4:295486/149‑1 (MQ=255)                                                                                                                                           cAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAAcc         <  6:41455/149‑1 (MQ=255)                                                                                                                                             gggg‑gATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGAc                                                                              >  7:134980/1‑78 (MQ=255)                                                                                                                                             gggg‑gATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGAc                                                                              <  8:134980/78‑1 (MQ=255)                                                                                                                                                g‑gATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGtt    <  7:208972/149‑1 (MQ=255)                                                                                                                                                  gATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGc                                               >  6:48258/1‑105 (MQ=255)                                                                                                                                                  gATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGc                                               <  5:48258/105‑1 (MQ=255)                                                                                                                                                  gATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTtg   <  4:274385/149‑1 (MQ=255)                                                                                                                                                   aTCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTtgt  <  2:240976/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  GCAGGTCGAAACGCAGCATCGCGCCCGGCAGG‑CTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTTGT  >  NC_002947/5867591‑5867883

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