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	98.2%	186.6 / NA	57	intergenic (‑70/+28)	PP_3394/PP_3395	3‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base G (33/22);  minor base T (0/1);  total (33/24)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GC‑GGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAG‑CGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACC  >  NC_002947/3845557‑3845838                                                                                                                                                     |                                                                                                                                              gctgctGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTc‑gc‑g                                                                                                                                              <  3:188484/148‑1 (MQ=255) gctgctGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GC‑gg                                                                                                                                             >  7:289350/1‑149 (MQ=255)     ctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GC‑gggct                                                                                                                                          >  2:113641/1‑145 (MQ=255)     ctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GC‑gggc                                                                                                                                           <  3:203505/147‑3 (MQ=255)     ctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGGCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GC‑gggc                                                                                                                                           <  6:247293/147‑3 (MQ=255)           gCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCTGGCTGATCCAt                                                                                                                                    <  5:331675/149‑1 (MQ=255)           gCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATcctc                                                                                                                                    <  1:272898/149‑3 (MQ=255)           gCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCAt                                                                                                                                    >  7:36301/1‑149 (MQ=255)                aTCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTc                                                                                                                                <  1:299647/148‑1 (MQ=255)                aTCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTc                                                                                                                                <  6:387148/148‑1 (MQ=255)                 tCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGtggc                                                                                                                              >  4:400174/1‑146 (MQ=255)                    cTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCa                                                                                                                           >  1:138541/1‑149 (MQ=255)                       gCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCCCGa                                                                                                                        >  6:426519/1‑149 (MQ=255)                        cAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGAt                                                                                                                       <  7:340730/149‑1 (MQ=255)                          ggCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATgg                                                                                                                     <  7:385488/149‑1 (MQ=255)                           gCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGa                                                                                                                    >  4:386077/1‑149 (MQ=255)                               tGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑Tc                                                                                                                <  1:113641/147‑1 (MQ=255)                                      cATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCgcgc                                                                                                       >  8:388529/1‑149 (MQ=255)                                      cATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCgcgc                                                                                                       >  2:32794/1‑149 (MQ=255)                                                 ccGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGcc                                                                                         >  4:387876/1‑149 (MQ=255)                                                    aCTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCATATACATAGTCAGCCAAGATGGA‑‑TCAGCCCGCGCCCA‑‑‑GCCAAGCCGGt                                                                                      >  5:286562/1‑149 (MQ=255)                                                          ttgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACg                                                                                 >  1:346487/1‑148 (MQ=255)                                                                 ttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGAAACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCaa                                                                           >  1:145265/1‑147 (MQ=255)                                                                   tttgttgTATTCGCGGTCTCTCTGTCTTACGGGTGGCGCCATATAAAGCGCGGTTTTCGGGGCATAAACGCCGTTTTGCGGG‑GGCTTACCCTTCGTCCGCCACTATGTA‑‑TACGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAg                                                                        <  8:36301/147‑1 (MQ=255)                                                                    tggttgtCATCGCGGTCTCTCTTTCTTAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCATAAAAGACGTGTTTC‑GCGGGCTTATCCATCGTCAGCAACTATGTA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAggt                                                                      <  2:346487/147‑1 (MQ=255)                                                                       ttgtCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGc                                                                                          >  3:130175/1‑126 (MQ=255)                                                                       ttgtCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAggtggt                                                                   >  2:355553/1‑149 (MQ=255)                                                                         gtCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTg                                                                                     >  1:306359/1‑129 (MQ=255)                                                                          tCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTc                                                                  >  1:175737/1‑147 (MQ=255)                                                                            aTCGCGGTCTCTCTGTCTGAAGCGTTCCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAg                                                               <  2:33807/148‑1 (MQ=255)                                                                               gcgGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACgcag                                                                              >  7:159665/1‑130 (MQ=255)                                                                                     tctctGTCTTAAGCGTTGCGCCAGACTAAGAGCTGGTTGAGGGGCATAAAAGACGTGTTTC‑GCGGGCTTATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccg                                                      <  6:96068/148‑1 (MQ=255)                                                                                      ctctGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCa                                                                                                                           >  2:332998/1‑83 (MQ=255)                                                                                      ctctGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCa                                                                                                                           <  1:332998/83‑1 (MQ=255)                                                                                       tctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagc                                                   >  7:183454/1‑149 (MQ=255)                                                                                           tctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGGTCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGAc                                               <  4:393058/149‑1 (MQ=255)                                                                                           tctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagccga                                                >  1:346655/1‑148 (MQ=255)                                                                                             tgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAg                                             >  3:241078/1‑149 (MQ=255)                                                                                             tgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAg                                             >  8:166003/1‑149 (MQ=255)                                                                                             tgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAg                                             >  6:185388/1‑149 (MQ=255)                                                                                              gAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACgcagc                                                                             <  1:265622/116‑1 (MQ=255)                                                                                              gAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACgcagc                                                                             >  2:265622/1‑116 (MQ=255)                                                                                                          cAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAG‑CGAGAGCGTGTGc                               >  3:336656/1‑149 (MQ=255)                                                                                                                        gATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCTTCAGCCACGATGGATCTCAGCCCGCGCCCTCTCCCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAG‑CGAGAGCGTGTGCCGGTCgcgc                      <  1:355553/149‑1 (MQ=255)                                                                                                                        gATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAG‑CGAGAGCGTGTGCCGGTCGCGCATg                   >  3:184974/1‑147 (MQ=255)                                                                                                                           gACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAG‑CGAGAGCGTGTGCCGGTCGCGCATGCAGAt              >  2:58365/1‑149 (MQ=255)                                                                                                                            ccGGGCATAAACGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCACCGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAG‑CGAGAGCGTGTGCCGGTCGCGCATGCAGAt              <  3:386077/147‑1 (MQ=255)                                                                                                                            aCGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGTA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAG‑CGAGAGCGTGTGCCGGTCGCGCATGCAGAt              <  7:388529/148‑1 (MQ=255)                                                                                                                            aCGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCCGAGAGCGTGTGCCGGTCGCGCATGCAGAt              <  1:32794/149‑1 (MQ=255)                                                                                                                             cGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAG‑CGAGAGCGTGTGCCGGTCGCGCATGCAGATGa            <  1:58365/149‑1 (MQ=255)                                                                                                                               ggCTGAAAAGACGTGTTTC‑GCGGGCTGATCCTGCTGCCGCCCCGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAG‑CGAGAGCGTGTGCCGGTCGCGCATGCAGATGACg          <  4:241078/149‑1 (MQ=255)                                                                                                                                gCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAG‑CGAGAGCGTGTGCCGGTCGCGCATGCAGATGACg          >  3:134342/1‑148 (MQ=255)                                                                                                                                   gAAAAGACGTGTTTC‑GCGGGCTGATACATCCTCAGCCACCATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACCCACCAACAGGTGGTCAAGCAGGCGCCGAGCCCACAG‑CGAGAGCGTGTGCCCGTCGCGCAAGCAGATGACGAAt       >  2:23499/1‑148 (MQ=255)                                                                                                                                    aaaaGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAG‑CGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAAtt      >  5:209926/1‑148 (MQ=255)                                                                                                                                    aaaaGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAG‑CGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTc     >  2:268359/1‑149 (MQ=255)                                                                                                                                      acgACGTGTTTC‑GCGGGCTTATCCTTGTCCTGCCACGTTGTA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAG‑CGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCAc   <  4:336656/147‑1 (MQ=255)                                                                                                                                       tgaggtgtTTC‑GCGGGCTGACCAATCGTCAGCCACGAGGGA‑‑TCTGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAG‑CGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCAcc  <  3:387876/145‑1 (MQ=255)                                                                                                                                           gtgtTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGCCAG‑CGAGAGCGTGTGCCGGTcgc                        >  3:278925/1‑123 (MQ=255)                                                                                                                                             gtTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAG‑CGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTc     <  7:245239/140‑1 (MQ=255)                                                                                                                                             gtTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAG‑CGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTc     >  8:245239/1‑140 (MQ=255)                                                                                                                                              tttC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAAc                                                                          <  6:430802/71‑1 (MQ=255)                                                                                                                                              tttC‑GCGGGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAAc                                                                          >  5:430802/1‑71 (MQ=255)                                                                                                                                                     |                                                                                                                                              GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GC‑GGCTGATCCATCGTCAGCCACGATGGA‑‑TCAGCCCGCGCCCT‑‑‑GCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAG‑CGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACC  >  NC_002947/3845557‑3845838

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