Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,080,989	T→A	14.4%	intergenic (+36/+26)	PP_5332 → / ← PP_5333	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,080,989	0	T	A	14.4%	133.8 / 15.2	62	intergenic (+36/+26)	PP_5332/PP_5333	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/25);  new base A (6/3);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

GCAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGG  >  NC_002947/6080845‑6081130                                                                                                                                                 |                                                                                                                                                gCAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGccc                                                                                                                                             >  4:48320/1‑149 (MQ=255)  cAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCcg                                                                                                                                            >  4:292296/1‑149 (MQ=255)  cAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCcg                                                                                                                                            >  5:130035/1‑149 (MQ=255)  cAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCcg                                                                                                                                            >  3:174620/1‑149 (MQ=255)    gCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCcgc                                                                                                                                           >  2:334072/1‑148 (MQ=255)      gAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑aa                                                                                                                                       <  3:292296/149‑1 (MQ=255)      gAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCCGCG‑aa                                                                                                                                       <  4:174620/149‑1 (MQ=255)         tGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑Aagag                                                                                                                                    <  3:154489/149‑1 (MQ=255)           cTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGc                                                                                                                                  <  1:347832/149‑1 (MQ=255)              tgCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCg                                                                                                                                <  1:135533/148‑1 (MQ=255)                  tCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGc                                                                                                                            >  2:211954/1‑148 (MQ=255)                     agaCAGCGAATGTGGGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAAtt                                                                                                                        <  7:41073/149‑1 (MQ=255)                      gaCAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCATCGCGGGTTTGCCCGCG‑AAGAGGCCGGGGCAATTa                                                                                                                       >  3:74597/1‑149 (MQ=255)                         aGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATc                                                                                                                     >  3:158105/1‑148 (MQ=255)                         aGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATc                                                                                                                     >  2:176468/1‑148 (MQ=255)                          gCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCa                                                                                                                    >  6:379499/1‑148 (MQ=255)                            gAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAg                                                                                                                 <  2:362767/149‑1 (MQ=255)                            gAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAg                                                                                                                 <  3:48320/149‑1 (MQ=255)                              aTGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCCGCG‑a                                                                                                                                        >  7:184050/1‑124 (MQ=255)                              aTGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCCGCG‑a                                                                                                                                        <  8:184050/124‑1 (MQ=255)                                               gtggAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCcaccac                                                                                               <  6:130035/148‑1 (MQ=255)                                                 ggAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTATTCCCCCCCCCt                                                                                            >  6:299087/1‑149 (MQ=255)                                                    gtgcaaGCCTGGGGCTATGCATGGGCGGGGCGTCCAGTTCGGTTGTTAAAATATCTGAGCCCTGCTCTTTTTTTCCGGGCCTCTTCGCGGGTTTGCCCGCG‑AAGTGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCAc                                                                                         <  4:74597/145‑1 (MQ=255)                                                       gAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCtt                                                                                      >  4:49656/1‑149 (MQ=255)                                                        aaGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTg                                                                                     <  6:73934/149‑1 (MQ=255)                                                        aaGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTAACGGCCTCTTCGCGGGCTTGCCCGGG‑CAGAGGCCGGGGCTATTTTCAAAGGTCTTCGTGTTGCACCACCCTGACCTTg                                                                                     >  5:42299/1‑149 (MQ=255)                                                          gCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTc                                                                                   <  5:379499/149‑1 (MQ=255)                                                           ccTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCt                                                                                  >  3:200590/1‑149 (MQ=255)                                                               gggCGAGTCATTCGCTCGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCa                                                                               <  5:299087/148‑1 (MQ=255)                                                                      tCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAgcgc                                                                       <  1:334072/149‑1 (MQ=255)                                                                        aTTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTGCCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCa                                                                      >  6:195781/1‑148 (MQ=255)                                                                        aTTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCg                                                                                                                                <  7:304555/90‑1 (MQ=255)                                                                        aTTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCg                                                                                                                                >  8:304555/1‑90 (MQ=255)                                                                        aTTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAt                                                                     <  1:176468/149‑1 (MQ=255)                                                                         ttCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑a                                                                   <  1:211954/149‑1 (MQ=255)                                                                                       aGATCGGTTGTTAAAATATCTGAGCCCTGCTCTTTTTTTCCCGGCCTCTTCGCGGGGTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTc                                                                                                     <  7:137506/102‑1 (MQ=255)                                                                                         aTCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCCGCG‑aa                                                                                                                                       >  3:25892/1‑66 (MQ=255)                                                                                         aTCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCCGCG‑aa                                                                                                                                       <  4:25892/66‑1 (MQ=255)                                                                                             gTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCAt                                               >  1:66210/1‑149 (MQ=255)                                                                                                gtgAAAATATCTGAGCCCTGCTCTTTTTTTGCCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTg                                             <  5:195781/148‑1 (MQ=255)                                                                                                gtgAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTg                                             <  7:358865/148‑1 (MQ=255)                                                                                                gtgAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTg                                             <  3:49656/148‑1 (MQ=255)                                                                                                    aaaTATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTATCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTGCTGAc                                        >  6:370772/1‑149 (MQ=255)                                                                                                           tGAGCCCTACTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATCAGGCCGCATCGGCCAGGTCATTGCTg                                          <  1:5828/141‑1 (MQ=255)                                                                                                           tGAGCCCTACTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTGCTg                                          >  2:5828/1‑140 (MQ=255)                                                                                                                   gCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTGCTGACCTtc                                    <  6:208632/138‑1 (MQ=255)                                                                                                                   gCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTGCTGACCTtc                                    >  5:208632/1‑138 (MQ=255)                                                                                                                   gCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAg                         >  6:152859/1‑149 (MQ=255)                                                                                                                       tttttttACCGGCCTCTTCGCGGGCTTGCCCGCGAAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTc                                                 >  1:267682/1‑122 (MQ=255)                                                                                                                       tttttttACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTc                                                 <  2:267682/121‑1 (MQ=255)                                                                                                                       tttttttACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTGCTGAcc                                       >  4:66289/1‑131 (MQ=255)                                                                                                                       tttttttACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTGCTGAcc                                       <  3:66289/131‑1 (MQ=255)                                                                                                                       tttttttACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTg                     >  1:293293/1‑149 (MQ=255)                                                                                                                                cGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGAcc            >  7:58432/1‑149 (MQ=255)                                                                                                                                cGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCCCTGAcc            >  7:190962/1‑149 (MQ=255)                                                                                                                                 ggCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCGAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGAcc            >  1:331719/1‑148 (MQ=255)                                                                                                                                  gCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCca          >  3:114672/1‑149 (MQ=255)                                                                                                                                       ttCGCGGGCAAGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCg     >  6:82975/1‑149 (MQ=255)                                                                                                                                       ttCGCGGGCAAGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCg     >  1:350290/1‑149 (MQ=255)                                                                                                                                       ttCGCGGGCAAGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCg     >  8:162805/1‑149 (MQ=255)                                                                                                                                          gcgGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCggg   <  4:200590/148‑1 (MQ=255)                                                                                                                                          gcgGGCTTGCCCGCG‑AAGAGGCCGGGGCAAGTATCAAAGGTCTGCGTCTTCCACGACCCTCACCTTGTCTGCAGCCAGCGCAT‑CGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGAGCCACAGCgggg  >  7:266342/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                GCAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCG‑AAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAT‑AGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGG  >  NC_002947/6080845‑6081130

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑