Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	556,986	A→C	10.2%	K12T (AAG→ACG)	rpsN →	30S ribosomal protein S14

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	556,986	0	A	C	10.2%	172.3 / 9.7	59	K12T (AAG→ACG)	rpsN	30S ribosomal protein S14
Reads supporting (aligned to +/- strand):  ref base A (29/24);  new base C (3/3);  total (32/27)
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.85e-01

CTCTGCGCGGTTTGGACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAGCGCCTCGCGCCTGCG  >  NC_002947/556841‑557133                                                                                                                                                  |                                                                                                                                                    ctctGCGCGGTTTGGACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGaa                                                                                                                                                     >  5:292451/1‑146 (MQ=255) ctctGCGCGGTTTGGACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGaa                                                                                                                                                     <  6:292451/146‑1 (MQ=255) ctctGCGCGGTTTGGACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGc                                                                                                                                                   <  8:281988/148‑1 (MQ=255)  tctGCGCGGTTTGGACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGt                                                                                                                                                 >  4:78783/1‑149 (MQ=255)   cggcgcGGTGAGGACATCACCCAGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCGGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTc                                                                                                                                                <  8:165082/147‑1 (MQ=255)     gcgcGGTTTGGACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGaa                                                                                                                                                     >  4:98661/1‑142 (MQ=255)     gcgcGGTTTGGACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGaa                                                                                                                                                     <  3:98661/142‑1 (MQ=255)      cgcgGTTTGGACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGTGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCa                                                                                                                                               <  6:629/147‑1 (MQ=255)       gcgGTTTGGACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCt                                                                                                                                            >  2:162305/1‑149 (MQ=255)           tttGGACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTAAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGTGAGCTGACGCTTCAGCTcgcg                                                                                                                                        <  8:71351/149‑1 (MQ=255)               gACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGTGAGCTGACGCTTCAGCTCGCGGttt                                                                                                                                    >  4:22172/1‑147 (MQ=255)               gACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTa                                                                                                                                     <  3:78783/148‑1 (MQ=255)                  aTCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGTGAGCTGACGCTTCAGCTcgcg                                                                                                                                        <  5:223592/142‑1 (MQ=255)                  aTCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGTGAGCTGACGCTTCAGCTcgcg                                                                                                                                        >  6:223592/1‑142 (MQ=255)                  aTCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCa                                                                                                                                  >  6:322736/1‑148 (MQ=255)                  aTCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCa                                                                                                                                  >  4:57682/1‑148 (MQ=255)                   tCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGTGAGCTGACGCTTCAGCTCGCGGtttttca                                                                                                                                <  3:22172/149‑7 (MQ=255)                   tCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCa                                                                                                                                  >  1:86686/1‑147 (MQ=255)                     aCCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAg                                                                                                                               <  1:165130/148‑1 (MQ=255)                      cccTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGtt                                                                                                                             >  6:132784/1‑149 (MQ=255)                       ccTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTc                                                                                                                            >  6:266377/1‑149 (MQ=255)                           accaccACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGc                                                                                                                          <  5:129699/147‑1 (MQ=255)                           accaccACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCt                                                                                                                         <  3:57682/148‑1 (MQ=255)                           accaccACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCt                                                                                                                         <  6:268096/148‑1 (MQ=255)                             caccacTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCt                                                                                                                         >  1:53937/1‑146 (MQ=255)                                       cGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTgct                                                                                                             >  5:69753/1‑148 (MQ=255)                                               cgacgaAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGaaa                                                                                                    <  7:261632/149‑1 (MQ=255)                                                 acgaAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGc                                                                                                  >  2:60114/1‑149 (MQ=255)                                                  cgaAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTAACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGc                                                                                                  >  5:262038/1‑148 (MQ=255)                                                   gaAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGa                                                                                                <  5:266377/149‑1 (MQ=255)                                                      ggCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCa                                                                                             >  6:222969/1‑149 (MQ=255)                                                       gCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTc                                                                                                                            <  1:2563/117‑1 (MQ=255)                                                       gCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTc                                                                                                                            >  2:2563/1‑117 (MQ=255)                                                        ccGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATc                                                                                           <  1:60114/149‑1 (MQ=255)                                                           cgcTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGaaa                                                                                                    >  7:139625/1‑137 (MQ=255)                                                           cgcTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGaaa                                                                                                    <  8:139625/137‑1 (MQ=255)                                                             ctctGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGCAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCaa                                                                                      <  5:290515/149‑1 (MQ=255)                                                                  ctgcGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAAcc                                                                                    <  5:222969/146‑1 (MQ=255)                                                                  ctgcGTGCTTTCAAGTTCCCGTTCCGCAACTGAATGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAAcc                                                                                    <  5:322736/146‑1 (MQ=255)                                                                     cGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGG‑AGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGc                                                                             >  1:170951/1‑149 (MQ=255)                                                                      gTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGc                                                                             <  5:132784/149‑1 (MQ=255)                                                                             cAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCtctc                                                                        >  8:225346/1‑147 (MQ=255)                                                                                     cGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCg                                                              >  3:110795/1‑149 (MQ=255)                                                                                       ttCCGGAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGtt                                                            <  6:93508/149‑1 (MQ=255)                                                                                          cGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTc                                                          <  6:325487/148‑1 (MQ=255)                                                                                             aaCTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTg                                                      >  5:347380/1‑149 (MQ=255)                                                                                                 gATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGc                                                     >  8:132883/1‑146 (MQ=255)                                                                                                 gATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGc                                                     <  7:132883/146‑1 (MQ=255)                                                                                                   ttGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTg                                                >  7:122350/1‑149 (MQ=255)                                                                                                      gAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCg                                             <  4:110795/149‑1 (MQ=255)                                                                                                            gAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTg                                        >  7:170702/1‑148 (MQ=255)                                                                                                            gAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGc                                       <  6:69753/149‑1 (MQ=255)                                                                                                               aaTGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAg                                     >  2:330748/1‑148 (MQ=255)                                                                                                                   gCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAagca                                <  1:162305/149‑1 (MQ=255)                                                                                                                    ccAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAagcag                               >  8:212391/1‑149 (MQ=255)                                                                                                                     cAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAagcagc                              >  5:321990/1‑149 (MQ=255)                                                                                                                        gaagaGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCAc                           >  5:22239/1‑149 (MQ=255)                                                                                                                                tGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGc                    >  8:316608/1‑148 (MQ=255)                                                                                                                                   aaaaCCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCg                                             >  1:101183/1‑120 (MQ=255)                                                                                                                                   aaaaCCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCg                                             <  2:101183/120‑1 (MQ=255)                                                                                                                                        cgtgAGCTGACGCTTCAGCTCGCGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAGCGCCTc           >  8:60703/4‑149 (MQ=255)                                                                                                                                                gAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAGCGCCTCGCGCCTgc   >  8:97816/1‑149 (MQ=255)                                                                                                                                                gAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAGCGCCTCGCGCCTgc   >  6:133172/1‑149 (MQ=255)                                                                                                                                                 aaGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAGCGCCTCGCGCCTgcg  <  2:86686/149‑1 (MQ=255)                                                                                                                                                 aaGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAGCGACTCGCGCCTgcg  >  7:109744/1‑149 (MQ=255)                                                                                                                                                  aGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCAGCAGAAGCAGCCACGCGACGCTAGCGCCTCGCGCCTgcg  >  1:113912/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CTCTGCGCGGTTTGGACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAGCGCCTCGCGCCTGCG  >  NC_002947/556841‑557133

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