Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,319,937	G→C	17.7%	R61G (CGC→GGC)	PP_4681 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,319,937	0	G	C	17.7%	147.5 / 22.3	62	R61G (CGC→GGC)	PP_4681	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/30);  new base C (4/7);  total (25/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGAGCGATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGGTGGTC  >  NC_002947/5319790‑5320084                                                                                                                                                    |                                                                                                                                                    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                                                                                                                                          <  2:273037/141‑1 (MQ=255)                   cAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGa                                                                                                                                          >  1:273037/1‑141 (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                                                                                                              <  2:246860/149‑1 (MQ=255)                                       cgTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAg                                                                                                              <  5:301049/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: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  >  3:20137/1‑149 (MQ=255)                                                                                                                                                   cGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGGTGGTc  <  4:219061/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    AAGGAGCGATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGGTGGTC  >  NC_002947/5319790‑5320084

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