Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	3,845,705	+G	intergenic (‑72/+26)	PP_3394 ← / ← PP_3395	3‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,845,703	1	.	G	96.7%	200.8 / ‑0.1	61	intergenic (‑70/+28)	PP_3394/PP_3395	3‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (35/24);  total (36/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑C‑GG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCAC  >  NC_002947/3845557‑3845837                                                                                                                                                       |                                                                                                                                        gctgctGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑cg‑c‑g                                                                                                                                        <  2:184858/148‑1 (MQ=255)  ctgctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑C‑ggg                                                                                                                                      >  5:399704/1‑148 (MQ=255)          ggCATGATCACTTGCAGGCTCTGCAAGCCCTCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGCCGTGTTT‑CG‑CGGG‑CTGATCCa                                                                                                                              <  5:299087/149‑1 (MQ=255)               gATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGt                                                                                                                          >  4:279287/1‑148 (MQ=255)               gATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGAGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CT‑CGGG‑CTGATCCATCGt                                                                                                                          >  7:209305/1‑148 (MQ=255)                  cACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGc                                                                                                                      >  4:127602/1‑149 (MQ=255)                    cTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGccc                                                                                                                    >  3:242219/1‑148 (MQ=255)                     ttGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCAc                                                                                                                   >  3:102060/1‑149 (MQ=255)                        cAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGa                                                                                                                 <  6:399704/148‑1 (MQ=255)                         aGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATg                                                                                                               >  7:24297/1‑149 (MQ=255)                          ggCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATgg                                                                                                              <  3:279287/149‑1 (MQ=255)                           gCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGa                                                                                                             <  8:24297/149‑1 (MQ=255)                           gCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGa                                                                                                             >  1:231642/1‑149 (MQ=255)                            ctctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGAt                                                                                                            <  2:347475/149‑1 (MQ=255)                              ctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑ATGATCCATCGTCAGCAACGATGGATCa                                                                                                          >  7:8864/1‑149 (MQ=255)                                gCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGc                                                                                                        >  1:51247/1‑149 (MQ=255)                                 cAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGATGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCa                                                                                                          <  5:389754/146‑1 (MQ=255)                                      cATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCgcgc                                                                                                  >  1:192649/1‑149 (MQ=255)                                         cgcgCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGAAGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCt                                                                                               >  5:88328/1‑149 (MQ=255)                                            gcAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGcc                                                                                            <  5:401965/149‑1 (MQ=255)                                                  cGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCg                                                                                      >  3:337660/1‑149 (MQ=255)                                                    aCTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CGAG‑GG‑CTGATCCATAGTCAGCCAACATGGATCAGCCCGCGCCCTgagaaccg                                                                                       >  5:29671/1‑139 (MQ=255)                                                    aCTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGt                                                                                    <  3:127602/149‑1 (MQ=255)                                                     cTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTGCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGt                                                                                    >  1:59666/1‑149 (MQ=255)                                                      ttCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTg                                                                                   >  5:203941/1‑148 (MQ=255)                                                        cATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGAt                                                                                                                                 >  4:5163/1‑100 (MQ=255)                                                        cATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGAt                                                                                                                                 <  3:5163/100‑1 (MQ=255)                                                           tgattgTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGTTCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACg                                                                               >  3:52065/1‑147 (MQ=255)                                                                 ttattgttGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCaa                                                                         >  2:325380/1‑147 (MQ=255)                                                                   attgttGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAgg                                                                     >  5:109626/1‑149 (MQ=255)                                                                       ttgtCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATACATAGTCAGCCACGATGTATCAGACCGCGCCCTTCCAGGCCGGTGCACGCAGCAACAggtgg                                                                  >  5:392966/1‑148 (MQ=255)                                                                       ttgtCATCGCGGGCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggtgg                                                                  <  2:59666/148‑1 (MQ=255)                                                                        tgtCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggtggt                                                                 >  7:58851/1‑148 (MQ=255)                                                                          tCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCa                                                               <  2:51247/148‑1 (MQ=255)                                                                            aTCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAg                                                             >  8:279387/1‑148 (MQ=255)                                                                            aTCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAg                                                             <  1:325380/148‑1 (MQ=255)                                                                             tCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCATACTAAGAGCTGTATGACGGGCATAAAACACGTGTTT‑CG‑C‑GGCCTTACCCATCGTCACCCACTATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAAc                                                                        <  1:102974/136‑1 (MQ=255)                                                                             tCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAg                                                             >  3:72010/1‑147 (MQ=255)                                                                             tCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATAAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAAc                                                                        >  2:102974/1‑136 (MQ=255)                                                                               gcgGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAgg                                                         >  6:166502/1‑149 (MQ=255)                                                                                cgGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGTGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCTCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGgc                                                        <  1:260708/149‑1 (MQ=255)                                                                                 ggTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCCGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGgcg                                                       >  5:287504/1‑149 (MQ=255)                                                                                  gTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTCATCCATCGTCAGCCACTATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGgcgc                                                      <  5:106389/149‑1 (MQ=255)                                                                                   tctct‑cTGTCTGAAGCGTTGCGCCAGACTAAGAGTTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgcc                                                     >  5:344753/1‑149 (MQ=255)                                                                                      ct‑ctGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgag                                                  >  1:260031/1‑149 (MQ=255)                                                                                       t‑ctgtctgAAGCGTGCCGCAATACTAAGACCTGTAGTAGGGGCATAAAACACGTTTTT‑CG‑CGGG‑CTGATCCATCGTCACCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccga                                                   <  6:109626/147‑1 (MQ=255)                                                                                           gtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagccga                                              >  6:149842/1‑149 (MQ=255)                                                                                            tctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGCTCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagccga                                              >  7:260491/1‑148 (MQ=255)                                                                                                aaGCGTTGCGCCAGACTAAGAGCTGGAGGACGGGCAGAAAACCCGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAg                                           <  2:231642/147‑1 (MQ=255)                                                                                                 aGCGTTGCGCCAGACTAAGAGCGGGATGAGGGGCATAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCg                                         <  8:257367/148‑1 (MQ=255)                                                                                                     ttGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGagag                                     >  1:287056/1‑149 (MQ=255)                                                                                                       gcgcCAGACTAAGAGCTGGATGACGGGCACAAACGCCGTGTTT‑CG‑CGGG‑CTGATCCATCGTCACCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgt                                  <  3:230177/149‑1 (MQ=255)                                                                                                         gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgtgt                                <  6:203941/149‑1 (MQ=255)                                                                                                          ccACACTAAGAGCTGGATGACGGGCACAAAAGCCGTGTTT‑CG‑CGGG‑CTGATCCATCGTCACCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGtgtg                               <  5:101922/149‑1 (MQ=255)                                                                                                               cTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCgg                           <  2:260031/148‑1 (MQ=255)                                                                                                                tAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTc                         >  2:103595/1‑149 (MQ=255)                                                                                                                tAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTc                         >  3:334136/1‑149 (MQ=255)                                                                                                                             tcGGGCAGAAATGACGTGTTG‑CG‑CGGG‑CTGATCCATGTTCAGCCACGATGGTTCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATg            <  4:335090/148‑1 (MQ=255)                                                                                                                               gggCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATg            >  3:335090/1‑147 (MQ=255)                                                                                                                               gggCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGAc          >  1:132780/1‑149 (MQ=255)                                                                                                                               gggCAGAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATAGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGAc          >  5:78659/1‑149 (MQ=255)                                                                                                                                    gAAAAGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAAtt     <  5:166502/149‑1 (MQ=255)                                                                                                                                       aaGACGTGTTT‑CG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCAc  <  7:279387/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                        GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCT‑CTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTT‑CG‑C‑GG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCAC  >  NC_002947/3845557‑3845837

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