Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,370,020	G→C	6.6%	A56G (GCC→GGC)	ramA ←	(R)‑stereoselective amidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,370,020	0	G	C	6.6%	185.2 / 4.0	61	A56G (GCC→GGC)	ramA	(R)‑stereoselective amidase
Reads supporting (aligned to +/- strand):  ref base G (23/34);  new base C (2/2);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

TCAATGCGATCCCTTGCGGGGTGACCAGCAAGGAGGTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCGGTATCGCCATCGCGACCGGAGAGTT  >  NC_002947/4369874‑4370167                                                                                                                                                     |                                                                                                                                                    tCAATGCGATCCCTTGCGGGGTGACCAGCAAGGAGGTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGc                                                                                                                                                    >  5:36221/1‑148 (MQ=255)   aaTGCGATCCCTTGCGGGGTGACCAGCAAGGAGGTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTGTCACACTCGGCCCGTTAACGGGCTCGGc                                                                                                                                                    <  1:172067/147‑1 (MQ=255)      gCGATCCCTTGCGGGGTGACCAGCAAGGAGGTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTc                                                                                                                                              >  5:225643/1‑149 (MQ=255)                 cGGGGTGACCAGCAAGGAGGTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGt                                                                                                                                    <  3:31828/148‑1 (MQ=255)                         ccAGGAAGGAGGTGTTGTATACGGTGCCGTAGGCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCTT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCTTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTg                                                                                                                           <  8:106007/149‑1 (MQ=255)                          cAGCAAGGAGGTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGGTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTgg                                                                                                                          >  1:97020/1‑149 (MQ=255)                               aggaggTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAAt                                                                                                                     >  1:195061/1‑149 (MQ=255)                               aggaggTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAAt                                                                                                                     <  2:149546/149‑1 (MQ=255)                                ggaggTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATc                                                                                                                    <  3:227754/149‑1 (MQ=255)                                 gaggTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATcc                                                                                                                   <  5:192662/149‑1 (MQ=255)                                    gtgtTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGCCGaaat                                                                                                                                               <  6:3785/118‑4 (MQ=255)                                    gtgtTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGCCGaaat                                                                                                                                               >  5:3785/1‑115 (MQ=255)                                      gttgtAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGTGTCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCAt                                                                                                                >  1:93143/1‑148 (MQ=255)                                      gttgtAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGAttt                                                                                                                                               >  5:329573/1‑116 (MQ=255)                                      gttgtAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGAttt                                                                                                                                               <  6:329573/116‑1 (MQ=255)                                       ttgtAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAg                                                                                                             >  4:202564/1‑149 (MQ=255)                                         gtAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGc                                                                                                            <  6:225643/148‑1 (MQ=255)                                         gtAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCt                                                                                                           >  2:118151/1‑149 (MQ=255)                                            aaCGGTGCCTGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTgg                                                                                                         <  5:316773/148‑1 (MQ=255)                                            aaCGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCGCGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTggg                                                                                                        >  7:144721/1‑149 (MQ=255)                                            aaCGGTGCCGGAGTCGTTCTCAACCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTgg                                                                                                         <  4:243136/148‑1 (MQ=255)                                             aCGGTGCCGGATTCGTGCTCTGGCATACGTATGAGCGGTGAGACGTGTCGT‑TCTTGTACGGCGCGTTGCAGGGGCTG‑GTCACTGGGGGCGTTAACGGGCTCGGCGATTTCGGCGTGGTGTTGGGCGCTGGCTAATCCCATCAGGTggtt                                                                                                       <  8:154934/149‑3 (MQ=255)                                                  gCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTg                                                                                                  >  3:32675/1‑149 (MQ=255)                                                   ccGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTgg                                                                                                 >  6:218898/1‑149 (MQ=255)                                                   ccGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTg                                                                                                  <  2:195061/148‑1 (MQ=255)                                                    cGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTggg                                                                                                <  4:32675/149‑1 (MQ=255)                                                               tcAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCAt                                                                                                                <  1:134020/122‑1 (MQ=255)                                                               tcAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCAt                                                                                                                >  2:134020/1‑122 (MQ=255)                                                                 atcccgACCAATCAACGCTGAGACGTTTCGT‑TCGTGTCCGGCGCGTTGCACGGCCTG‑CACCCTCGGCCCGTTAACGGGCGCGGCGATTTCGGCTAGCGGTTGGGCGCTGGCAACTCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGAt                                                                                   <  6:89195/144‑1 (MQ=255)                                                                  gCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGCCGAAATCGCCGAGccc                                                                                                                                     <  8:192435/98‑3 (MQ=255)                                                                  gCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGCCGAAATCGCCGAGccc                                                                                                                                     >  7:192435/1‑96 (MQ=255)                                                                   ccAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCTAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATcag                                                                                        <  4:963/142‑1 (MQ=255)                                                                   ccAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCg                                                                                 <  3:202564/149‑1 (MQ=255)                                                                   ccAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCg                                                                                 <  2:93143/149‑1 (MQ=255)                                                                   ccAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCAAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATcag                                                                                        >  3:963/1‑142 (MQ=255)                                                                        ccAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATc                                                                            <  2:97020/149‑1 (MQ=255)                                                                            tCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCGCGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACg                                                                        <  8:144721/149‑1 (MQ=255)                                                                                cGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGgcg                                                                     >  4:299593/1‑148 (MQ=255)                                                                                     agaCGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTa                                                               <  1:225074/149‑1 (MQ=255)                                                                                       aCGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGc                                                             >  2:238145/1‑149 (MQ=255)                                                                                              gt‑tcgtGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGc                                                       <  5:218898/148‑1 (MQ=255)                                                                                               t‑tcgtGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCt                                                     <  3:154084/149‑1 (MQ=255)                                                                                                       aCGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATg                                               >  5:158562/1‑148 (MQ=255)                                                                                                        cGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTa                                             >  7:281371/1‑149 (MQ=255)                                                                                                          gcgcGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACg                                           >  3:127802/1‑149 (MQ=255)                                                                                                               ttGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCa                                      <  3:299593/149‑1 (MQ=255)                                                                                                                   aCGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGtt                                  <  6:15897/149‑1 (MQ=255)                                                                                                                    cGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTa                                             <  8:119543/137‑1 (MQ=255)                                                                                                                    cGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTa                                             <  8:115169/137‑1 (MQ=255)                                                                                                                    cGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTa                                             >  7:115169/1‑137 (MQ=255)                                                                                                                    cGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTa                                             >  7:119543/1‑137 (MQ=255)                                                                                                                    cGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTta                                 <  6:158562/149‑1 (MQ=255)                                                                                                                        cTG‑CACACTCGGCCAGTTAACGGGCTCGGCGATTTCGGCTAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGCGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAg                             <  7:120871/149‑1 (MQ=255)                                                                                                                            cacaCTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGc                            <  4:127802/147‑1 (MQ=255)                                                                                                                                   ggCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCGGTATCg                     <  1:238145/147‑1 (MQ=255)                                                                                                                                    gCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCGGTATCGCCa                  >  5:243529/1‑149 (MQ=255)                                                                                                                                    gCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCGGTATCGCCa                  >  6:328828/1‑149 (MQ=255)                                                                                                                                     cccGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCGGTATCGCCAt                 <  1:118151/149‑1 (MQ=255)                                                                                                                                      ccGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCGGTATCGCCATc                <  5:328828/149‑1 (MQ=255)                                                                                                                                          tAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCACATGGGTCTCTGGGAAAATCAGCAGATCGGTATCAACGGCGCAGGTAGCGATGGCCTCCAATGTACGTTCCAGGTAATAACCGGAATCGCCATCGCGa            >  7:74931/1‑149 (MQ=255)                                                                                                                                               ggCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCGGTATCGCCATCGCGACCGga       <  6:243529/149‑1 (MQ=255)                                                                                                                                                    ggCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCGGTATCGCCATCGCGACCGGAGAGtt  >  8:225155/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TCAATGCGATCCCTTGCGGGGTGACCAGCAAGGAGGTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGT‑TCGTGTACGGCGCGTTGCACGGCCTG‑CACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCGGTATCGCCATCGCGACCGGAGAGTT  >  NC_002947/4369874‑4370167

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