Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,319,943	G→T	17.5%	L59M (CTG→ATG)	PP_4681 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,319,943	0	G	T	17.5%	152.9 / 21.4	62	L59M (CTG→ATG)	PP_4681	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/29);  new base T (4/7);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGGTGGTCGGCAATA  >  NC_002947/5319798‑5320091                                                                                                                                                  |                                                                                                                                                     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                                                                                                                                                <  8:101798/142‑1 (MQ=255)           cAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGAt                                                                                                                                                >  7:101798/1‑142 (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                                                                                                                                              <  3:277561/90‑1 (MQ=255)                                                                 acacGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATcc                                                                                                                                              >  4:277561/1‑90 (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                                                                                                          >  6:280064/1‑75 (MQ=255)                                                                                                                    ccAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACgg                                                                                                          <  5:280064/75‑1 (MQ=255)                                                                                                                         tCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTc                           >  5:159059/1‑149 (MQ=255)                                                                                                                           gatgatgatCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTcg                                <  3:193878/142‑1 (MQ=255)                                                                                                                           gatgatgatCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTcg                                >  4:193878/1‑142 (MQ=255)                                                                                                                           gatgatgatCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGa                         <  5:211731/149‑1 (MQ=255)                                                                                                                             aatggatCTTCCTGCCCCGCATGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCgcgc                             >  3:205922/5‑143 (MQ=255)                                                                                                                             aatggatCTTCCTGCCCCGCATGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGAt                        >  1:147224/5‑148 (MQ=255)                                                                                                                             aatggatCTTCCTGCCCCGCATGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGAt                        >  1:142426/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)                                                                                                                                           cGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGGTGGTc         <  4:219061/149‑1 (MQ=255)                                                                                                                                           cGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGGTGGTc         >  3:20137/1‑149 (MQ=255)                                                                                                                                             gggCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGGTGGTCgg       >  6:102923/1‑149 (MQ=255)                                                                                                                                              tgCAGGAAGATCCATTGCAGCCACTGCTAAAAACTGAGGGTATCCACGGCTAACGGCACCGTACTGGCCAGTGCCTCGGCGCGTGGCGGCTCGTCATCCCACCAACCTTGCACCTTCAGCTCTCGCTCGATCAAGAGCAGGTGGTCGGc      <  5:22058/148‑1 (MQ=255)                                                                                                                                              ggCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGGTGGTCgg       >  3:231831/1‑148 (MQ=255)                                                                                                                                                  ggAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGGTGGTCGGCAata  <  3:139087/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     ATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGGTGGTCGGCAATA  >  NC_002947/5319798‑5320091

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