Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,489,483	A→G	8.6%	intergenic (+40/‑36)	PP_4824 → / → PP_4825	two‑component system sensor histidine kinase/response regulator/MarC family membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,489,483	0	A	G	8.6%	162.4 / 4.9	59	intergenic (+40/‑36)	PP_4824/PP_4825	two‑component system sensor histidine kinase/response regulator/MarC family membrane protein
Reads supporting (aligned to +/- strand):  ref base A (31/22);  new base G (2/3);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATGGACGGACACATGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCAAGGAGCGCA  >  NC_002947/5489336‑5489630                                                                                                                                                    |                                                                                                                                                     aTGGACGGACACATGGCCAAACCGGTGGAGTTGTCGCAGTTGTGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACaa                                                                                                                                                     <  6:363064/149‑1 (MQ=255)    gacggacACATGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAgg                                                                                                                                                   <  5:240871/148‑1 (MQ=255)       ggacACATGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCcgc                                                                                                                                               >  2:309571/1‑149 (MQ=255)         acacATGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGCGCGGCCTTGTGTCGCGAGGCcgctc                                                                                                                                             <  8:54064/149‑3 (MQ=255)             aTGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGc                                                                                                                                         <  1:360380/149‑1 (MQ=255)              tGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCa                                                                                                                                        >  2:117711/1‑149 (MQ=255)                 ccAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATa                                                                                                                                      >  3:319716/1‑148 (MQ=255)                  cAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATaaa                                                                                                                                    >  3:300529/1‑149 (MQ=255)                       cGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAAccc                                                                                                                                 >  1:378996/1‑147 (MQ=255)                          tGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTTACCGATATCCAATGCTATTGCAGGAGCGGCCTCGCGACACAAGTACGCGCCGGCATAAAACAAAgt                                                                                                                             >  5:111986/1‑147 (MQ=255)                                tGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTg                                                                                                                      >  6:200467/1‑149 (MQ=255)                                  tCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGc                                                                                                                     <  5:304406/148‑1 (MQ=255)                                          gCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATCCTc                                                                                                             >  4:255541/1‑148 (MQ=255)                                            gTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCAt                                                                                                          >  2:382582/1‑149 (MQ=255)                                                  tGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAgttg                                                                                                    >  5:381256/1‑149 (MQ=255)                                                     ttCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCCAGGAACACTGCTCATGCTCCATGAgttgtt                                                                                                  >  2:297408/1‑148 (MQ=255)                                                     ttCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAgttgtt                                                                                                  <  6:196958/148‑1 (MQ=255)                                                      tCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTTGTGTCGCGAGGCcgc                                                                                                                                               >  8:174256/1‑102 (MQ=255)                                                      tCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTTGTGTCGCGAGGCcgc                                                                                                                                               <  7:174256/102‑1 (MQ=255)                                                             tGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTa                                                                                         >  6:387710/1‑149 (MQ=255)                                                              gggCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTAc                                                                                        >  4:217387/1‑149 (MQ=255)                                                                 ccAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCt                                                                                      >  8:205570/1‑148 (MQ=255)                                                                   aaTCACCGGGAAGCCAAGGTGGACCCGGCACTAACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACCAGGCCGCGCCGGCATAAACCCAAGGAACACAGCTGATGCTCCATGAGTAGTTCAGCGTCTACCTGa                                                                                    >  7:18221/1‑148 (MQ=255)                                                                         cGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCt                                                                             >  1:60805/1‑149 (MQ=255)                                                                         cGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCt                                                                             <  3:217387/149‑1 (MQ=255)                                                                          gggAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTc                                                                            <  1:382582/149‑1 (MQ=255)                                                                               gCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGc                                                                              >  3:301819/1‑142 (MQ=255)                                                                               gCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGc                                                                              <  4:301819/142‑1 (MQ=255)                                                                                ccAAGGTGGACCCGGCACATAAGTCCTGACCGATATACAATGCTCTTGCAGGAGCGGCCTATCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑Gctc                                                                     <  5:53056/149‑1 (MQ=255)                                                                                  aaGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTa                                                                   >  8:41846/1‑149 (MQ=255)                                                                                        gACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGccc                                                             >  8:87496/1‑149 (MQ=255)                                                                                         aCCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCg                                                            <  5:200467/149‑1 (MQ=255)                                                                                             ggCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑g                                                                        <  4:212176/133‑1 (MQ=255)                                                                                             ggCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑g                                                                        >  3:212176/1‑133 (MQ=255)                                                                                             ggCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGttct                                                        >  1:263517/1‑149 (MQ=255)                                                                                               cacaTACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGttcttc                                                      <  5:387710/149‑1 (MQ=255)                                                                                                   tACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGc                                                  >  6:220679/1‑149 (MQ=255)                                                                                                    aCGTCCTGACCGATATCCAATGATCTTGCAGGAGCGGACTCGCGACACAAGGCCGAGCCGGGATAACCCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCATCGTCTACCTGTAAATGCTCGT‑GCTCTACAGAACGTTCTTCGTGCt                                                 >  6:259521/1‑149 (MQ=255)                                                                                                       tCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTc                                              >  3:20875/1‑149 (MQ=255)                                                                                                         cTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCgt                                            >  7:157742/1‑149 (MQ=255)                                                                                                           gACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCgtgt                                          >  8:34072/1‑149 (MQ=255)                                                                                                           gACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCgt                                            <  5:308096/147‑1 (MQ=255)                                                                                                           gACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCgt                                            >  6:308096/1‑147 (MQ=255)                                                                                                                    ccAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGc                                 >  3:333578/1‑149 (MQ=255)                                                                                                                     cAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGcc                                >  6:66635/1‑149 (MQ=255)                                                                                                                             ttGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGt                        >  8:332071/1‑149 (MQ=255)                                                                                                                               gCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTgg                      >  5:19658/1‑149 (MQ=255)                                                                                                                                cAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGc                     <  7:34072/149‑1 (MQ=255)                                                                                                                                cAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGc                     <  2:60805/149‑1 (MQ=255)                                                                                                                                 aGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGcc                    >  6:314019/1‑149 (MQ=255)                                                                                                                                    cgCGGCCTTGTGTCGCGAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCg                                                     <  3:324719/112‑1 (MQ=255)                                                                                                                                    cgCGGCCTTGTGTCGCGAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCg                                                     >  4:324719/2‑113 (MQ=255)                                                                                                                                        gCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAg             <  1:138334/149‑1 (MQ=255)                                                                                                                                         ccTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTCGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAg             <  1:309571/149‑1 (MQ=255)                                                                                                                                         ccTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGc            <  5:220679/149‑1 (MQ=255)                                                                                                                                            cgcgACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCAAg         <  4:333578/149‑1 (MQ=255)                                                                                                                                             gcgACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCaa          <  7:205570/147‑1 (MQ=255)                                                                                                                                             gcgACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCAAg         <  2:378996/148‑1 (MQ=255)                                                                                                                                              cgACACAAGGCCGCGCCGGCATAAACCCAAGGAACTCTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCAAGGa       <  3:255541/149‑1 (MQ=255)                                                                                                                                              cgACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCAAGGa       <  4:319716/149‑1 (MQ=255)                                                                                                                                                  acaAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCAAgg        >  2:179453/1‑144 (MQ=255)                                                                                                                                                  acaAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCAAgg        <  1:179453/144‑1 (MQ=255)                                                                                                                                                    aaGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGATTCATCAGCCTGCCCCGTGGCCACTCCAGCAAGGAGCgca  >  2:324393/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ATGGACGGACACATGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGT‑GCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCAAGGAGCGCA  >  NC_002947/5489336‑5489630

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