Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	796,816	T→A	6.5%	intergenic (+18/‑172)	PP_0681 → / → PP_0682	hypothetical protein/inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	796,816	0	T	A	6.5%	179.5 / 3.1	62	intergenic (+18/‑172)	PP_0681/PP_0682	hypothetical protein/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (29/29);  new base A (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

ATGGTTATCAAGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTT  >  NC_002947/796677‑796961                                                                                                                                             |                                                                                                                                                   atgGTTATCAAGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCAAAGCGGcccc                                                                                                                                           <  5:119826/149‑2 (MQ=255)     ttATCAATTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATcgc                                                                                                                                       <  6:274921/149‑1 (MQ=255)      tATCAAGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATcgc                                                                                                                                       >  2:84813/1‑148 (MQ=255)      tATCAAGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATcgc                                                                                                                                       >  2:135906/1‑148 (MQ=255)       aTCAAGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATcgc                                                                                                                                       <  6:210834/147‑1 (MQ=255)       aTCAAGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGCTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCAAcgc                                                                                                                                       >  5:210834/1‑147 (MQ=255)          aaGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGac‑ac                                                                                                                                 <  2:102249/149‑1 (MQ=255)          aaGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGac‑a                                                                                                                                  <  3:106391/148‑1 (MQ=255)          aaGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCAGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGac‑a                                                                                                                                  >  2:62287/1‑148 (MQ=255)               gCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGc                                                                                                                            <  2:57663/149‑1 (MQ=255)                        ccGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCcac                                                                                                                   <  1:171060/149‑1 (MQ=255)                        ccGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCca                                                                                                                    >  4:13259/1‑148 (MQ=255)                        ccGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCca                                                                                                                    >  3:314095/1‑148 (MQ=255)                         cGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCca                                                                                                                    >  5:267831/1‑147 (MQ=255)                                cTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAg                                                                                                           >  5:21845/1‑149 (MQ=255)                                   gtgtCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCAt                                                                                                        >  5:279725/1‑149 (MQ=255)                                   gtgtCTGCATTTGCAGGCGCAGTTGCAGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACa                                                                                                            <  1:62287/146‑1 (MQ=255)                                        tGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCTCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCAAAGCGGccccg                                                                                                                                          <  5:332183/111‑3 (MQ=255)                                        tGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCt                                                                                                   >  3:301422/1‑149 (MQ=255)                                        tGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCt                                                                                                   >  1:71997/1‑149 (MQ=255)                                        tGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCAAAGCGGccccg                                                                                                                                          >  6:332183/1‑109 (MQ=255)                                          cATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCtt                                                                                                  <  6:21845/148‑1 (MQ=255)                                           aTTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGc                                                                                                >  5:99083/1‑149 (MQ=255)                                           aTTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGc                                                                                                <  4:119002/149‑1 (MQ=255)                                            tttGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCt                                                                                               >  1:125477/1‑149 (MQ=255)                                                cAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCtgtg                                                                                            <  4:157349/148‑1 (MQ=255)                                                cAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTgg                                                                                           <  1:135906/149‑1 (MQ=255)                                                  ggCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGa                                                                                         <  2:125477/149‑1 (MQ=255)                                                         ttGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGc                                                                                  <  6:99083/149‑1 (MQ=255)                                                             ggcggATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGACTACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCtt                                                                               <  1:84813/149‑1 (MQ=255)                                                             ggcggATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTtg                                                                              >  6:143676/1‑149 (MQ=255)                                                              gcggATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTtgt                                                                             <  6:267831/149‑1 (MQ=255)                                                                   tCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTcg                                                                         <  4:314095/148‑1 (MQ=255)                                                                           gAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGTGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGc                                                                >  2:177960/1‑149 (MQ=255)                                                                           gAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTgg                                                                                           <  1:45619/122‑1 (MQ=255)                                                                           gAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTgg                                                                                           >  2:45619/1‑122 (MQ=255)                                                                           gAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCAACGGCGTGCTGTGGGa                                                                                         >  8:31790/1‑124 (MQ=255)                                                                                gtgtCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCg                                                            >  1:90092/1‑148 (MQ=255)                                                                                    cGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGagca                                                       <  5:143676/149‑1 (MQ=255)                                                                                           tGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCtt                                                                               >  2:215971/1‑118 (MQ=255)                                                                                           tGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCtt                                                                               <  1:215971/118‑1 (MQ=255)                                                                                                     ccAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCt                                      <  2:71997/149‑1 (MQ=255)                                                                                                     ccAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCt                                      >  6:294666/1‑149 (MQ=255)                                                                                                             aTCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAgg                              >  1:222953/1‑149 (MQ=255)                                                                                                             aTCCTGACACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAACCCCGCTCCCACAGGTACAGCATCGGCATGCTGAGGGAGTGTGGCCATGTCTCGCGATAAGGCTGCCGAGCAGCCCCCCCCCATTTGCTGTAACAg                               >  4:260846/1‑148 (MQ=37)                                                                                                              tCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGTGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGt                             <  1:177960/149‑1 (MQ=255)                                                                                                              tCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGt                             >  8:30061/1‑149 (MQ=255)                                                                                                                cTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGtgg                           <  4:301422/149‑1 (MQ=255)                                                                                                                cTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGtgg                           >  3:84182/1‑149 (MQ=255)                                                                                                                   agaCCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGtggc                          <  3:260846/147‑1 (MQ=255)                                                                                                                     aCCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGtggctggc                      <  7:30061/149‑1 (MQ=255)                                                                                                                      ccTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGtggctggc                      <  4:84182/148‑1 (MQ=255)                                                                                                                       cTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCta                    <  5:294666/149‑1 (MQ=255)                                                                                                                        tGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGCTCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGAGGGCTGGCta                    >  1:273270/1‑148 (MQ=255)                                                                                                                        tGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCGTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCta                    >  5:191295/1‑148 (MQ=255)                                                                                                                        tGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCCCTCCCACACGTACAGCCTCGGCGTGCAGTGGGAGGGGGGCCTTGTGTCGCCATCAGGCTGCGGCGCAGCCCCCGGGCTTTTGCTGTAAACGGTGGCTGGCta                    >  8:23246/1‑148 (MQ=37)                                                                                                                         gaagaaTGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGc                                       <  5:307646/128‑1 (MQ=255)                                                                                                                         gaagaaTGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGc                                       >  6:307646/1‑128 (MQ=255)                                                                                                                         cacgaaGGCCGGGC‑CCGCTTGCCGGCCCAGCGCCAC‑ACACGGCCGCTCCCACCGGTACACCAGCGGCTGGCTGTGGGAGGGGGGCCTTGTGGCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCtata                  <  6:191295/146‑1 (MQ=255)                                                                                                                          aagaaTGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGGGTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTGTTGCTGTAACAGGTGGCTGGCTATAc                 >  2:86257/1‑149 (MQ=255)                                                                                                                                   tggg‑CCGCAAAGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCg        >  3:102113/2‑149 (MQ=255)                                                                                                                                         cGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCtttt   <  2:222953/149‑1 (MQ=255)                                                                                                                                            tttGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCCGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCttttt  >  6:94239/1‑147 (MQ=255)                                                                                                                                             |                                                                                                                                                   ATGGTTATCAAGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGG‑CCGCTTTGCGGCCCATCGCGAC‑ACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTT  >  NC_002947/796677‑796961

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