Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,560,678	A→C	6.3%	S291A (TCC→GCC)	PP_4048 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,560,678	0	A	C	6.3%	217.0 / 2.5	72	S291A (TCC→GCC)	PP_4048	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/36);  new base C (3/2);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGCATGTACTGATCACTGACCACCTCCCCATAAATACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGC  >  NC_002947/4560534‑4560819                                                                                                                                                 |                                                                                                                                              cGAGCATGTACTGATCACTGACCACCTCCCCATAAATACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGTCCGAGGGTGCGCTCGGCCCGc                                                                                                                                           <  8:235437/149‑1 (MQ=255) cGAGCATGTACTGATCACTGACCACCTCCCCATAAATACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGc                                                                                                                                           >  3:56236/1‑149 (MQ=255)           cTGATCACTGTCCACCTCCCCATAAATACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTaag                                                                                                                                 <  4:181557/149‑2 (MQ=255)           cTGATCACTGACCACCTCCCCATAAATACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGCCCGCTCAAGGTAAc                                                                                                                                 >  2:262024/1‑149 (MQ=255)           cTGATCACTGACCACCTCCCCATAAATACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAAc                                                                                                                                 >  5:234231/1‑149 (MQ=255)            tGATCACTGACCACCTCCCCATAAATACGATTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGGGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACg                                                                                                                                <  8:283157/149‑1 (MQ=255)                 aCTGACCACCTCCCCATAAATACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGTCCGAGGGTGCGCTCGGCCCGc                                                                                                                                           >  4:144626/1‑133 (MQ=255)                 aCTGACCACCTCCCCATAAATACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGTCCGAGGGTGCGCTCGGCCCGc                                                                                                                                           <  3:144626/133‑1 (MQ=255)                       caccTCCCCATAAATACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGc                                                                                                                     <  4:56236/149‑1 (MQ=255)                         ccTCCCCATAAATACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACCGGCGGGCCGAGCGCACCCTCGGCCCGCTCAACGTAACCTGGGCGCAAGcc                                                                                                                    >  1:264429/1‑148 (MQ=255)                          cTCCCCATAAATACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGa                                                                                                                  <  4:159773/149‑1 (MQ=255)                            ccccATAAATACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATg                                                                                                                >  2:265062/1‑149 (MQ=255)                              ccATAAATACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGc                                                                                                               >  8:213420/1‑148 (MQ=255)                              ccATAAATACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGc                                                                                                               >  8:254374/1‑148 (MQ=255)                               cATAAATACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCa                                                                                                             <  4:2707/149‑1 (MQ=255)                                   aaTACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCagag                                                                                                          <  6:38838/148‑1 (MQ=255)                                    aTACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCaa                                                                                                                                       >  8:150852/1‑118 (MQ=255)                                    aTACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCaa                                                                                                                                       <  7:150852/118‑1 (MQ=255)                                    aTACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGt                                                                                                        >  5:34233/1‑149 (MQ=255)                                        ccTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATg                                                                                                                <  1:252146/137‑1 (MQ=255)                                         cTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCAccc                                                                                                   <  6:234231/149‑1 (MQ=255)                                            cTGCGCCCGGGTAGCCGTCTCGTGATAGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGcc                                                                                                >  7:64629/1‑149 (MQ=255)                                              gcgcCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGc                                                                                                               <  1:136172/132‑1 (MQ=255)                                              gcgcCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGc                                                                                                               >  2:136172/1‑132 (MQ=255)                                                          ccGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGc                                                                                   >  3:193868/1‑148 (MQ=255)                                                            gTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCAt                                                                                <  7:254374/149‑1 (MQ=255)                                                            gTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCAt                                                                                <  6:34233/149‑1 (MQ=255)                                                               gcgTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATggg                                                                             >  3:222621/1‑149 (MQ=255)                                                               gcgTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATggg                                                                             <  2:115624/149‑1 (MQ=255)                                                                 gTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGAACCCTAGGACCGCTCAAGGTAACGTGGGCGCATGACGATGCCAGAGGTCACACGCCCAGCACTCATTGCAATGGGCt                                                                           >  8:66975/1‑149 (MQ=255)                                                                    aTCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCg                                                                        >  5:299237/1‑149 (MQ=255)                                                                       gAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAg                                                                                                                                      <  5:169849/84‑1 (MQ=255)                                                                       gAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAg                                                                                                                                      >  6:169849/1‑84 (MQ=255)                                                                       gAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCt                                                                     >  2:219136/1‑149 (MQ=255)                                                                       gAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCt                                                                     >  6:92655/1‑149 (MQ=255)                                                                          cacGTGCACATTGGCGCCATGCGCAGGGATTGGCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGGCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCtttt                                                                  <  3:76559/149‑1 (MQ=255)                                                                           acGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTg                                                                 >  8:10989/1‑149 (MQ=255)                                                                                cacaTTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCAt                                                            <  6:299237/149‑1 (MQ=255)                                                                                  caTTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCg                                                                        >  6:170811/1‑135 (MQ=255)                                                                                  caTTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCg                                                                        <  5:170811/135‑1 (MQ=255)                                                                                    cggggtgcATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCAc                                                        <  2:243264/142‑1 (MQ=255)                                                                                       gcgcCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCAccc                                                      <  5:92655/148‑1 (MQ=255)                                                                                          ccATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGc                                                  >  8:43249/1‑149 (MQ=255)                                                                                          ccATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCAGGccggggtg                                                   >  1:212928/1‑141 (MQ=255)                                                                                              gcgcAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCtttt                                                                  >  4:14741/1‑129 (MQ=255)                                                                                              gcgcAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCtttt                                                                  <  3:14741/129‑1 (MQ=255)                                                                                                  aGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCaa                                          >  6:329545/1‑149 (MQ=255)                                                                                                    cGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAAcc                                        <  7:213420/149‑1 (MQ=255)                                                                                                         ctccgCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGTGCGCATGCAGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCAGTGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTgcg                                    <  5:192038/144‑1 (MQ=255)                                                                                                          tCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTgcgc                                   <  7:43249/148‑1 (MQ=255)                                                                                                            ggCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCggg                                >  4:127817/1‑149 (MQ=255)                                                                                                             gCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGgggg                                >  5:108263/1‑148 (MQ=255)                                                                                                             gCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCggg                                >  1:270503/1‑148 (MQ=255)                                                                                                             gCGACGTGGGCACGGGAGGGCCGTGCGCACCCTCGGACCGCTCAATGTAACGTGGGAGCATGCCGATGCAAGAGGTCACCCGCCCAGCACTCATGGCCATGGTCTCCGGCTTTTGCGCATGCACCCCTGCCTGCGCAACCTGCGCTGGc                               <  8:64629/149‑1 (MQ=255)                                                                                                                aCGTGGGCACGGGCTGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCgg                                 <  8:245523/144‑1 (MQ=255)                                                                                                                aCGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCgggagct                            >  7:245523/1‑145 (MQ=255)                                                                                                                aCGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTACGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAgg                            <  1:219136/149‑1 (MQ=255)                                                                                                                   tGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCg                         <  4:193868/149‑1 (MQ=255)                                                                                                                    gggCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGc                        >  1:221889/1‑149 (MQ=255)                                                                                                                     ggCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCt                       >  2:289523/1‑149 (MQ=255)                                                                                                                      gCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGGGGGCAGGTCGCTg                      <  6:108263/149‑1 (MQ=255)                                                                                                                       cACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTgg                     >  3:136112/1‑149 (MQ=255)                                                                                                                        aCGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGc                    <  4:125868/149‑1 (MQ=255)                                                                                                                              gggcCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTc              >  7:134512/1‑149 (MQ=255)                                                                                                                                gcCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCa            <  7:10989/149‑1 (MQ=255)                                                                                                                                    aGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCAc        <  8:134512/149‑1 (MQ=255)                                                                                                                                        cACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTGAGCGGCCCAGCACTCATGGCCATGGTCTCCGGCTTTTGCGCATGCACCGCGGCCTGCGGAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTg    <  2:264429/149‑1 (MQ=255)                                                                                                                                        cACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTg    <  4:136112/149‑1 (MQ=255)                                                                                                                                         acacTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCAATCATGGACATGGGATCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGTCCTGCTCCAGCACCTTgg   >  4:204782/4‑149 (MQ=255)                                                                                                                                         aCCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTgg   <  2:270503/149‑1 (MQ=255)                                                                                                                                          cccTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGc  <  1:289523/149‑1 (MQ=255)                                                                                                                                          cccTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGc  <  4:222621/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                              CGAGCATGTACTGATCACTGACCACCTCCCCATAAATACCCTTCTGCGCCCGGGTCGCCGTCGCGTGATCGAACACGTGCACATTGGCGCCATGCGCAGCGATTGTCGGCGACGTGGGCACGGGCGGGCCGAGCGCACCCTCGGACCGCTCAAGGTAACGTGGGCGCATGCCGATGCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGC  >  NC_002947/4560534‑4560819

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑