Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,215,414	G→C	6.7%	F543L (TTC→TTG)	PP_3697 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,215,414	0	G	C	6.7%	187.6 / 3.9	60	F543L (TTC→TTG)	PP_3697	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/29);  new base C (2/2);  total (29/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 = 8.59e-01

CTCGAAGAAATGCTCCTTACCAGACATCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTGCAGGTGATTGTGTCTGGAGTTGATCCACTGAAGGGCGGCTTTGGCCACACCAGCAGGAACCAGAAC  >  NC_002947/4215270‑4215560                                                                                                                                                   |                                                                                                                                                   ctcGAAGAAATGCTCCTTACCAGACATCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGaa                                                                                                                                                  <  5:86541/147‑1 (MQ=255) ctcGAAGAAATGCTCCTTACCAGACATCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAgg                                                                                                                                                >  2:303363/1‑149 (MQ=255) ctcGAAGAAATGCTCCTTACCAGACATCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAg                                                                                                                                                 <  8:99761/148‑1 (MQ=255)  tcGAAGAAATGCTCCTTACCAGACATCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAgg                                                                                                                                                >  7:179902/1‑148 (MQ=255)   cGAAGAAATGCTCCTTACCAGACATCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGc                                                                                                                                               >  2:271283/1‑148 (MQ=255)    gaagaaATGCTCCTTACCAGACATCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGc                                                                                                                                               <  1:146759/147‑1 (MQ=255)        aaaTGCTCCTTACCAGACATCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCtt                                                                                                                                         <  6:62461/149‑1 (MQ=255)        aaaTGCTCCTTACCAGACATCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCtt                                                                                                                                         >  5:62461/1‑149 (MQ=255)                cTTACCAGACATCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTCAAGGCacacc                                                                                                                                          >  3:67072/1‑136 (MQ=255)                cTTACCAGACATCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTCAAGGCacacc                                                                                                                                          <  4:67072/140‑5 (MQ=255)                     cAGACATCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCATTGAAGGCAAGCTTCCGTGCGAAACCa                                                                                                                            >  7:162378/1‑149 (MQ=255)                       gACATCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCaa                                                                                                                                             >  2:253369/1‑130 (MQ=255)                       gACATCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCaa                                                                                                                                             <  1:253369/130‑1 (MQ=255)                          aTCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCt                                                                                                                        <  6:100137/148‑1 (MQ=255)                          aTCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCt                                                                                                                        >  5:100137/1‑148 (MQ=255)                          aTCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTCAAGGCacaccc                                                                                                                                         >  4:127045/1‑126 (MQ=255)                          aTCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTCAAGGCacaccc                                                                                                                                         <  3:127045/131‑6 (MQ=255)                                        ggTCATCGCATGTGGCATTGCATGCAAATCATACGCGTCAGACACGCAATTGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTTCCACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCAc                                                                                                            <  1:271283/148‑1 (MQ=255)                                        ggTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGt                                                                                                          <  8:179902/148‑1 (MQ=255)                                         gTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCa                                                                                                             <  7:119581/144‑1 (MQ=255)                                         gTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCa                                                                                                             >  8:119581/1‑144 (MQ=255)                                             tCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTc                                                                                                     >  7:61023/1‑148 (MQ=255)                                                cATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTg                                                                                                                                                    <  4:208864/98‑1 (MQ=255)                                                cATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTg                                                                                                                                                    >  3:208864/1‑98 (MQ=255)                                                   gtgGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCt                                                                                                >  2:143967/1‑147 (MQ=255)                                                    tgGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGa                                                                                              <  3:279317/148‑1 (MQ=255)                                                    tgGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGa                                                                                              <  2:166829/148‑1 (MQ=255)                                                    tgGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAg                                                                                             <  5:144143/149‑1 (MQ=255)                                                    tgGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAg                                                                                             >  7:256672/1‑149 (MQ=255)                                                        aTTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATc                                                                                         <  7:245603/149‑1 (MQ=255)                                                            catgcaAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTcctt                                                                                                                                     <  3:104665/101‑3 (MQ=255)                                                            catgcaAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGt                                                                                                                                     >  4:104665/1‑101 (MQ=255)                                                            catgcaAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAAcc                                                                                     >  5:206703/1‑149 (MQ=255)                                                            catgcaAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAAcc                                                                                     >  1:208591/1‑149 (MQ=255)                                                            catgcaAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAAc                                                                                      <  2:121030/148‑1 (MQ=255)                                                               gcaAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTc                                                                                   >  4:11813/1‑148 (MQ=255)                                                                caAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAAt                                                                                                                <  8:29619/118‑1 (MQ=255)                                                                caAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAAt                                                                                                                >  7:29619/1‑118 (MQ=255)                                                                 aaaGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTc                                                                                   >  6:166638/1‑146 (MQ=255)                                                                  aaGCATCCGCGTAAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCATT‑‑ACGGAACGGCTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCATCTACCTGAGGATCAACCTCgag                                                                                >  3:308890/1‑148 (MQ=255)                                                                          gcgTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGAc                                                                       <  8:162378/149‑1 (MQ=255)                                                                             tCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGt                                                                     >  7:59343/1‑148 (MQ=255)                                                                               aGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTg                                                                    >  6:63928/1‑147 (MQ=255)                                                                               aGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTGc                                                                   <  1:143967/148‑1 (MQ=255)                                                                                gACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTGCAg                                                                 <  2:208591/149‑1 (MQ=255)                                                                                   acGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTGCAGGt                                                               >  1:27635/1‑148 (MQ=255)                                                                                   acGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTGCAGGTg                                                              <  6:206703/149‑1 (MQ=255)                                                                                     gCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTGCAGGTGAt                                                            >  4:204515/1‑149 (MQ=255)                                                                                       aaTGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTGCAgg                                                                >  2:227512/1‑143 (MQ=255)                                                                                       aaTGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTGCAgg                                                                <  1:227512/143‑1 (MQ=255)                                                                                            cGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTGCAGGTGATTGTGTCt                                                     <  1:303363/149‑1 (MQ=255)                                                                                                      ccGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTGCAGGTGATTGTGTCTGGAGTTGATc                                           <  8:61023/149‑1 (MQ=255)                                                                                                         cgATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTGCAGGTGATTGTGTCTGGAGTTGATCCAc                                        <  3:204515/149‑1 (MQ=255)                                                                                                                  cTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTGCAGGTGATTGTGTCTGGAGTTGATCCACTGAAGggcg                               <  8:59343/149‑1 (MQ=255)                                                                                                                   ttCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGAc                                                                       >  7:37697/1‑108 (MQ=255)                                                                                                                   ttCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGAc                                                                       <  8:37697/108‑1 (MQ=255)                                                                                                                          gCTTT‑‑ACGGAACCGGTGTGCCTTGAAGGCAAGCTTCCGTTCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTGCAGGTGATTGTGTCTGGAGTTGATCCACTGAAGGGCGGCTTTGGc                       <  8:256672/149‑1 (MQ=255)                                                                                                                                          gtgtGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTGCAGGTGATTGTGTCTGGAGTTGATCCACTGAAGGGCGGCTTTGGCCACACCAGCAGGaa         >  3:297294/1‑149 (MQ=255)                                                                                                                                             tgCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTGCAGGTGATTGTGTCTGGAGTTGATCCACTGAAGGGCGGCTTTGGCCACACCAGCAGGAACCa      <  5:63928/149‑1 (MQ=255)                                                                                                                                               ccTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTAc                                                                                                  >  8:303277/1‑55 (MQ=255)                                                                                                                                                 ttGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTGCAGGTGATTGTGTCTGGAGTTGATCCACTGAAGGGCGGCTTTGGCCACACCAGCAgg           <  8:319753/140‑1 (MQ=255)                                                                                                                                                 ttGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTGCAGGTGATTGTGTCTGGAGTTGATCCACTGAAGGGCGGCTTTGGCCACACCAGCAGGaagatcgga  >  7:319753/1‑142 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CTCGAAGAAATGCTCCTTACCAGACATCAACCCTTCTTTGGTCATCGCATGTGGCATTGCATGCAAAGCATCCGCGTCAGACACGCAATGGCGGTCGTTGGCCGCGATGAGGGCTTCAACAGCTTT‑‑ACGGAACGGGTGTGCCTTGAAGGCAAGCTTCCGTGCGAAACCATCCTGGAAAAATGCACGTGCTTCTACCTGAGGATCAACCTCGAGAATTCTGACGTGCAGGTGATTGTGTCTGGAGTTGATCCACTGAAGGGCGGCTTTGGCCACACCAGCAGGAACCAGAAC  >  NC_002947/4215270‑4215560

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