Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,149 A→G 17.5% intergenic (+103/‑210) cyoE → / → alaC protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,1490AG17.5% 113.9 / 16.0 53intergenic (+103/‑210)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (21/21);  new base G (4/5);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.40e-01

GACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAG  >  NC_002947/955002‑955297
                                                                                                                                                          |                                                                                                                                                    
gACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACAAg                                                                                                                                                     >  6:103376/1‑149 (MQ=255)
gACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCCGCCTCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                                      <  1:26638/148‑1 (MQ=255)
       aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTg                                                                                                                                              <  1:44815/149‑1 (MQ=255)
        gACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGt                                                                                                                                             <  4:18114/149‑1 (MQ=255)
        gACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGt                                                                                                                                             <  4:293625/149‑1 (MQ=255)
                cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                     <  3:52502/149‑1 (MQ=38)
                      tgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAg                                                                                                                                           >  3:94/1‑137 (MQ=255)
                      tgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAg                                                                                                                                           <  4:94/137‑1 (MQ=255)
                          gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTTAAACCCCTACTTAAAGGTTGCTGTATGGGCGGGTTTTCCcgc                                                                                                                           >  3:110820/1‑149 (MQ=255)
                          gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTCCTACAAGGTTACTGTAGGAGCGGGTTTACCcgc                                                                                                                           >  4:70236/1‑149 (MQ=255)
                          gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCTCCGACAACGTTCCTGTAGGAGCGGGGTTTCCcgc                                                                                                                           >  1:8926/1‑149 (MQ=39)
                              aCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCcgcg                                                                                                                          >  4:244503/1‑146 (MQ=255)
                               ccTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                    <  1:60231/135‑1 (MQ=37)
                               ccTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                    >  2:60231/1‑135 (MQ=37)
                                cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTATTTTTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGt                                                                                                                                   >  7:45013/1‑140 (MQ=17)
                                cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTATTTTTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGt                                                                                                                                   <  8:45013/140‑1 (MQ=17)
                                    gcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                                   >  8:225635/1‑147 (MQ=255)
                                         cTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCACCTACAAGGTTACTGTTTGAGCGGGTTTTACCGCGAACAGGCCAg                                                                                                                >  1:293710/1‑145 (MQ=255)
                                              aCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCTCCAAAAAGGTTGCGGTATGAGGGGGGTTAACCGCGAATAGGCCAGCACAGaa                                                                                                         >  3:109057/1‑148 (MQ=255)
                                                          ttCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                    >  4:69233/1‑108 (MQ=21)
                                                          ttCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                    <  3:69233/108‑1 (MQ=21)
                                                             actaaaaCCCCGGCAACTGGGCGGGGTTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGTGTAAACCCGCTCCTAAAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGAt                                                                                        <  4:110820/146‑1 (MQ=37)
                                                              cGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTACTACAAAGTTACTGTACGATCCGGTTTAACCCCGAAGACGCCAGCAAAGAAAAATAAAAAAAACTGaat                                                                                       >  8:66980/1‑147 (MQ=255)
                                                                       cGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaaa                                                                              <  5:103376/149‑1 (MQ=255)
                                                                       cGGCCACTCTGCAGGGGTTTTTTTTTGGGCGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCCCGGTAAAACCCTCTCATACAATTTACCTTTAGTAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaaa                                                                              <  3:65110/149‑1 (MQ=37)
                                                                         gCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTTCCTACAAGGTCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAg                                                                                                           <  1:241658/119‑1 (MQ=25)
                                                                           cACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                    >  3:31696/1‑91 (MQ=17)
                                                                           cACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                    <  4:31696/91‑1 (MQ=17)
                                                                            aCTGTGCCGGGGTTTTTTTTTGGGGTTAT‑‑‑‑‑TTTTCTTGTGCGGGCCTCTTCGCGGTTAAACCCGCTCTTACATGTTTCCTTTTGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGAt                                                                         <  2:8926/149‑1 (MQ=38)
                                                                            aCTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGAt                                                                         <  3:244503/149‑1 (MQ=255)
                                                                             cTGGGCGGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATc                                                                        <  2:196010/149‑1 (MQ=255)
                                                                                       gttttttttttGGGGTTTT‑‑‑‑‑TTTTCCTGGGCCGGCGGCTGCGCGGTTAACCCCGCTCAACCACGGTTCCTTTTGGCGGGGTTTAACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGtt                                                               <  2:280024/147‑1 (MQ=255)
                                                                                         tttttttttGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCTGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCt                                                             <  7:225635/149‑1 (MQ=255)
                                                                                                 tGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGTAAAACCCGCTCCTACAAGTTTCCTTTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATaaaa                                                      <  3:127192/148‑1 (MQ=38)
                                                                                                   gggtTAT‑‑‑‑‑TTTTCTTGTTCCGGCCTCTTCTCGTTTAAACCCTATATTACAAGTTTCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAg                                                   <  5:116719/149‑1 (MQ=37)
                                                                                                                      gtgCCGGCCTCTTCGCGGTTAAACCCGCTCCAACAATGTTCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGAtt                                                                                       <  5:64885/99‑1 (MQ=25)
                                                                                                                       tgCCGGCCTCTTCGCGGGTAAACCCGATCCTACAAGGTTCCTGTAGTAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGAt                                                                         <  1:252110/112‑1 (MQ=255)
                                                                                                                                ctTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTTCCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCt                           >  3:242568/1‑149 (MQ=255)
                                                                                                                                    gcgGTTAAACCCGCTCCTAAAAGGTTCCTGTAGTAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGtt                        <  1:120212/148‑1 (MQ=255)
                                                                                                                                    gcgGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCt                           <  2:74669/145‑1 (MQ=255)
                                                                                                                                    gcgGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCt                           >  1:74669/1‑145 (MQ=255)
                                                                                                                                       ggTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCgg                    >  6:131063/1‑149 (MQ=255)
                                                                                                                                           aaCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCAc                 >  8:179592/1‑148 (MQ=255)
                                                                                                                                            aCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTa               >  5:41882/1‑149 (MQ=255)
                                                                                                                                              ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCta                                                            >  2:237095/1‑102 (MQ=255)
                                                                                                                                              ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCta                                                            <  1:237095/102‑1 (MQ=255)
                                                                                                                                              ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGaaa                                               <  2:112793/115‑1 (MQ=255)
                                                                                                                                              ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGaaa                                               >  1:112793/1‑115 (MQ=255)
                                                                                                                                              ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc             >  1:117283/1‑149 (MQ=255)
                                                                                                                                                 cTCCTACAAGTTTCCTTTAGTAGGGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc          <  7:276016/149‑1 (MQ=255)
                                                                                                                                                     cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGa       >  2:37588/1‑149 (MQ=255)
                                                                                                                                                      tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt      >  2:230955/1‑149 (MQ=255)
                                                                                                                                                      ctgaaGTTTCCTGTAGGAGCTGGTTTACACGCGAAGAGGCCAGCACAGATAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGa       <  8:207243/145‑1 (MQ=255)
                                                                                                                                                       aCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt      >  1:54031/1‑148 (MQ=255)
                                                                                                                                                        cAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt      <  7:179592/147‑1 (MQ=255)
                                                                                                                                                         aaGGTTCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCa   <  3:70236/149‑1 (MQ=255)
                                                                                                                                                         aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAATATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt     >  5:230987/1‑147 (MQ=255)
                                                                                                                                                         aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATATGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTACTGAtag    >  7:207243/1‑146 (MQ=255)
                                                                                                                                                          aGGTTCCTGTAGTTGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg  <  7:293662/149‑1 (MQ=255)
                                                                                                                                                          |                                                                                                                                                    
GACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAG  >  NC_002947/955002‑955297

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: