Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,867,733:1	+C	7.9%	intergenic (‑28/+24)	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,733	1	.	C	7.9%	238.3 / 14.2	63	intergenic (‑28/+24)	thyA/lgt	thymidylate synthase/prolipoprotein diacylglyceryl transferase
Reads supporting (aligned to +/- strand):  ref base . (28/30);  new base C (3/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

CTTCTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGG‑CTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTT  >  NC_002947/5867586‑5867881                                                                                                                                                      |                                                                                                                                                     cTTCTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGG‑CTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAgggg                                                                                                                                                       >  8:321598/1‑148 (MQ=255)  ttCTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGG‑CTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑Ga                                                                                                                                                    >  7:300103/1‑149 (MQ=255)    cTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGG‑CTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑Ga                                                                                                                                                    <  7:32194/147‑1 (MQ=255)      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                                                       <  5:361476/140‑1 (MQ=255)                                                                                                         tCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCgg                                                       >  6:361476/1‑140 (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                                                                         >  6:346574/1‑96 (MQ=255)                                                                                                                                   gTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAg                                                                         <  5:346574/96‑1 (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                                                                            <  8:134980/78‑1 (MQ=255)                                                                                                                                                  gggg‑gATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGAc                                                                            >  7:134980/1‑78 (MQ=255)                                                                                                                                                     g‑gATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGtt  <  7:208972/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CTTCTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGG‑CTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTT  >  NC_002947/5867586‑5867881

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