Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,319,933	A→C	17.4%	M62R (ATG→AGG)	PP_4681 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,319,933	0	A	C	17.4%	149.9 / 24.8	63	M62R (ATG→AGG)	PP_4681	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/31);  new base C (4/7);  total (25/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GCGGAAGGAGCGATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGG  >  NC_002947/5319786‑5320079                                                                                                                                                    |                                                                                                                                                   tcGGAAGGAGCGATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCTGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCAt                                                                                                                                                   <  6:250459/148‑1 (MQ=255)     aaGGAGCGATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCgg                                                                                                                                               >  1:201577/1‑149 (MQ=255)      aGGAGCGATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCggg                                                                                                                                              <  3:155846/149‑1 (MQ=255)        gAGCGATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGc                                                                                                                                            <  5:94624/149‑1 (MQ=255)             aTCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAgg                                                                                                                                         >  4:50703/1‑147 (MQ=255)                aaTTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCCTGCCCCGCATGAAGa                                                                                                                                     <  7:230371/148‑1 (MQ=255)                aaTTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCACCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGAt                                                                                                                                    <  1:226055/149‑1 (MQ=255)                    gATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCCTGCCCCGCATGAAGATcat                                                                                                                                 <  4:157477/148‑3 (MQ=255)                     aTCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCCTGCCCCGCATGAAGATcatc                                                                                                                                <  7:56059/148‑4 (MQ=255)                     aTCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGGGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCAt                                                                                                                                <  5:249520/148‑1 (MQ=255)                      tCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATACATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTg                                                                                                                              >  8:42163/1‑149 (MQ=255)                       cAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGa                                                                                                                                     >  1:273037/1‑141 (MQ=255)                       cAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGa                                                                                                                                     <  2:273037/141‑1 (MQ=255)                       cAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGAt                                                                                                                                    >  7:101798/1‑142 (MQ=255)                       cAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGAt                                                                                                                                    <  8:101798/142‑1 (MQ=255)                         aaaGCTGCCTACAGGCGGCGTAGCTCACGACTTTGCTCCTGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGc                                                                                                                                            <  4:45567/132‑1 (MQ=255)                                ccAACACGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGc                                                                                                                    <  8:125863/149‑1 (MQ=255)                                 cAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTTGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCCTGCCCCGCATGAAGATcatca                                                                                                                               <  8:122168/137‑5 (MQ=255)                                 cAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCCTGCCCCGCATGAAGATcatca                                                                                                                               >  7:122168/1‑133 (MQ=255)                                 cAACAGGCGGCGTAGATCAAGACTTTGCTCCGGAAGGTCGGTAAACACGTTTTCCGCCATGACAATGATGATCGATGAATTCGGCAGCGGGTGCCCCAGCTCGATGAAGATCTTCATCAGGGGCAGGAAGATCGATTGCAGCAACTGCt                                                                                                                   >  6:22058/1‑149 (MQ=255)                                          gcgTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGa                                                                                                          <  3:50703/149‑1 (MQ=255)                                           cgTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAg                                                                                                         <  5:301049/149‑1 (MQ=255)                                           cgTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAg                                                                                                         <  2:246860/149‑1 (MQ=255)                                            gTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAgg                                                                                                        <  8:80675/149‑1 (MQ=255)                                             tAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAggg                                                                                                       <  8:175459/149‑1 (MQ=255)                                                        tttGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGc                                                                                             >  1:239836/1‑148 (MQ=255)                                                            cTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAAc                                                                                         <  2:201577/148‑1 (MQ=255)                                                            cTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACg                                                                                        <  8:273166/149‑1 (MQ=255)                                                             tCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACgg                                                                                       >  5:138395/1‑149 (MQ=255)                                                               cggccggTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAACCTGAGGGTATCCACGGCGAACGGCa                                                                                     >  6:235793/1‑149 (MQ=255)                                                                             acacGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATcc                                                                                                                                  >  4:277561/1‑90 (MQ=255)                                                                             acacGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATcc                                                                                                                                  <  3:277561/90‑1 (MQ=255)                                                                             acacGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGt                                                                       <  2:239836/149‑1 (MQ=255)                                                                                cGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGcc                                                                    <  5:85903/149‑1 (MQ=255)                                                                                               cAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTggcggc                                                     >  6:178880/1‑149 (MQ=255)                                                                                                 aGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTc                                                   <  5:178880/149‑1 (MQ=255)                                                                                                  ggATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCg                                                  >  7:129356/1‑149 (MQ=255)                                                                                                   gATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGt                                                 >  3:219061/1‑149 (MQ=255)                                                                                                       cTCGATACATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCAt                                              <  7:42163/148‑1 (MQ=255)                                                                                                        tCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATc                                             >  6:27986/1‑148 (MQ=255)                                                                                                          gATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATc                                             <  5:27986/146‑1 (MQ=255)                                                                                                          gATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCcca                                          >  3:66900/1‑149 (MQ=255)                                                                                                              cATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCTTCCCACCaa                                      >  3:154852/1‑149 (MQ=255)                                                                                                                   ggCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTg                                 >  4:84782/1‑149 (MQ=255)                                                                                                                   ggCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTg                                 <  3:84782/149‑1 (MQ=255)                                                                                                                      aGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCAc                              <  6:138395/149‑1 (MQ=255)                                                                                                                                ccAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACgg                                                                                              <  5:280064/75‑1 (MQ=255)                                                                                                                                ccAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACgg                                                                                              >  6:280064/1‑75 (MQ=255)                                                                                                                                     tCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTc               >  5:159059/1‑149 (MQ=255)                                                                                                                                       gatgatgatCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTcg                    >  4:193878/1‑142 (MQ=255)                                                                                                                                       gatgatgatCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTcg                    <  3:193878/142‑1 (MQ=255)                                                                                                                                       gatgatgatCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGa             <  5:211731/149‑1 (MQ=255)                                                                                                                                         aatggatCTTCCTGCCCCGCATGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCgcgc                 >  3:205922/5‑143 (MQ=255)                                                                                                                                         aatggatCTTCCTGCCCCGCATGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGAt            >  1:142426/5‑148 (MQ=255)                                                                                                                                         aatggatCTTCCTGCCCCGCATGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGAt            >  1:147224/5‑148 (MQ=255)                                                                                                                                         aatggatCTTCCTGCCCCGCATGAAGATCAATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCgcgc                 <  4:205922/139‑1 (MQ=255)                                                                                                                                           atgatCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCaa         <  8:129356/149‑1 (MQ=255)                                                                                                                                             gatCTTCCTGCCCCGCATGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAaga       <  2:142426/149‑1 (MQ=255)                                                                                                                                             gatCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCaa         >  4:139087/1‑147 (MQ=255)                                                                                                                                              atCTTCCTGCCCCGCATGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAaga       <  2:147224/148‑1 (MQ=255)                                                                                                                                                   cATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAagag      <  6:76719/144‑1 (MQ=255)                                                                                                                                                   cATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAgg  <  6:159059/148‑1 (MQ=255)                                                                                                                                                   cATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAgg  <  4:66900/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCGGAAGGAGCGATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGG  >  NC_002947/5319786‑5320079

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