Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,319,938	G→C	17.4%	P60P (CCC→CCG) ‡	PP_4681 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,319,938	0	G	C	17.4%	146.2 / 18.5	62	P60P (CCC→CCG) ‡	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 = 7.48e-01

AAGGAGCGATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGGTGGTCGG  >  NC_002947/5319790‑5320086                                                                                                                                                     |                                                                                                                                                     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                                        <  3:84782/149‑1 (MQ=255)                                                                                                               ggCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTg                                        >  4:84782/1‑149 (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)                                                                                                                                                     gggCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGGTGGTCgg  >  6:102923/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AAGGAGCGATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCCCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGCCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGGTGGTCGG  >  NC_002947/5319790‑5320086

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