Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,153 T→C 12.5% intergenic (+107/‑206) cyoE → / → alaC protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,1530TC12.5% 129.7 / 9.8 57intergenic (+107/‑206)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (20/29);  new base C (2/5);  total (22/35)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-02

GCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGCC  >  NC_002947/955006‑955300
                                                                                                                                                        |                                                                                                                                                    
gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                     <  2:359155/149‑1 (MQ=255)
gCCAGACCGCTGCGGATGTGCGGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                     <  8:247704/149‑1 (MQ=255)
 ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTcc                                                                                                                                                    >  7:178431/1‑149 (MQ=255)
 ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCTTCAAGGTTcc                                                                                                                                                    >  4:208409/1‑149 (MQ=255)
  ctatttgCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                                   <  1:309103/143‑1 (MQ=255)
  cAGACCGCTGCGGCTGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGTAAACCCCGGCCACTGGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                                   <  6:361979/149‑1 (MQ=255)
      ccGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAg                                                                                                                                               <  8:51803/149‑1 (MQ=255)
       cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                                   >  2:309103/1‑144 (MQ=255)
            cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGTTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                         <  8:109271/149‑1 (MQ=255)
            cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                         <  1:337936/149‑1 (MQ=38)
            cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                         <  5:156814/149‑1 (MQ=38)
            cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCg                                                                                                                                          <  1:9285/148‑1 (MQ=38)
            cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGc                                                                                                                                           <  7:161293/147‑1 (MQ=255)
             gtaTGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACTAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                        <  1:233788/147‑1 (MQ=37)
             ggATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                        <  4:61479/149‑1 (MQ=37)
              gATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGt                                                                                                                                       <  7:139948/149‑1 (MQ=255)
                    cTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAccc                                                                                                                                 <  4:14946/149‑1 (MQ=255)
                     tgatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTCCAACAGGAACCATGTAGGATCGGGTTTACCcg                                                                                                                                >  2:41298/1‑149 (MQ=37)
                      gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGACAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAGTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCcgc                                                                                                                               <  5:263995/149‑1 (MQ=255)
                       atgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCcgc                                                                                                                               <  1:178304/148‑1 (MQ=255)
                           cttACCCACGCTGAAAGTGCTGCCTTTTTCACTAAAACCCGGGCCAGTTTGGGGGGTTTTTTTTTTGGGGTTATTTTTCGTGTGCCGGCCTCTTGGGGGTAAAACCCGCTCCTACAAGGTTCCTTTTGTTGGGGTTTTACCCGCGAaga                                                                                                                          <  6:43139/147‑1 (MQ=255)
                                    gcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAAACGCTACTAAAATGTTAATTTAGGAGCGGGTTTACCCGCGAAGAGGCCAGc                                                                                                                   >  5:72804/1‑148 (MQ=37)
                                      tGAAACTGCTGCCTGTTTCACGAAAACCCCGCCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCAATTCGCGGGTTAACCCCCTCCTACAAGGTTACAGTAGGAGCGCGTTTAACCGCGAAACGGCCAGca                                                                                                                  >  1:38584/1‑146 (MQ=255)
                                          aCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAGGTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGAACCTACAAGGTTACTTTATGAGCTGCTTTACAAGCCAAAACGCCAGCAAAGaa                                                                                                             >  1:103185/1‑148 (MQ=25)
                                           ctgctgCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCGGCTTCAAGGTTTCTTTTGGAGCGGGGTTTCCCGGGAAGAGGGCCGCACAGaaaaa                                                                                                          >  3:243710/1‑149 (MQ=38)
                                              ctgcCTGTTTCACCAAAACCCCGGCCACTGGGCGGGGTTTTTTTTTTTGGGGTTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGTAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaa                                                                                                          <  6:72804/147‑1 (MQ=25)
                                                  cTGTTTCACGAAACCCCGGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCGGGCGTCTTGGCGTGTAAACCCGCTCCAACATGGTTCGTTTTGGTGCGGGTTTACCCGCGAATAGGCCAGCACAGaaaaataaaaat                                                                                                   <  2:38584/149‑1 (MQ=255)
                                                      aaaaaaaaaaCCCCGGCACAGTGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTTAACCCGCTGCTACAAGGTTGCGGTTGGGGCGGGGTTAACCGGGAAGAGGCCCGCCCAGAAAAATAAAAATACCt                                                                                               >  1:81758/7‑149 (MQ=25)
                                                       tCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTTAACCCCCGGCGACCAGGGTGGGGTTGGAGGGGGGTTTCCCGGGAAGCGGGCCGGCCAGAAAAATAAAAATTCCt                                                                                               >  7:62319/1‑149 (MQ=255)
                                                       tCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCt                                                                                               <  8:178431/148‑1 (MQ=255)
                                                       tCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                         >  3:235071/1‑106 (MQ=255)
                                                       tCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                         <  4:235071/106‑1 (MQ=255)
                                                                 cccGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCACCGACAAGGTTGCTGTATGAGCGGGTTTACCCGCGAAGAGGCCCGCACAGAAAAATTAAAAAAACTGATTTTGaaa                                                                                     >  8:64559/1‑148 (MQ=37)
                                                                    ggCCCCTGGGCGGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTGCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAAt                                                                                 <  6:300256/149‑1 (MQ=255)
                                                                     gCCACTGTGCCGGGGTTTTTTTTTGGGGTAATTTTTCCTGTGCCGGCCTCTTCGCGGTAAAACCCGCTCCTACAAGTTTCCTGTAGTAGCGTGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATa                                                                                <  7:303171/149‑1 (MQ=38)
                                                                         cTGTGCCGGGTTTTTTTTTTTGGGGTATTTTTCCTTTGCCGGCCTCTTCTCGGTTAAACCCTCTCCTACAATGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATc                                                                            <  1:34599/149‑1 (MQ=37)
                                                                            cctctctggTTTTTTTTTGGGGGTATTTTTCCTTTGCCGGCCTCTTCTCGGGTAAACCCGCTCATACAAGTATCCTGTAGTAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGAt                                                                             <  3:30829/139‑1 (MQ=25)
                                                                             ggCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGCCCTCTTCGCGGTAAAACCCGCTCCAACAAGTTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCt                                                                           <  3:208409/146‑1 (MQ=38)
                                                                                 ggggTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCACTTCGCGGGTAAACACGAGAACTACAAGGTTGGTGTATGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaata                                                                                                        >  1:99251/1‑115 (MQ=25)
                                                                                        ttttttttGGGGGTATTTTTCCTGTCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCt                                                                <  8:62319/148‑1 (MQ=255)
                                                                                                            tcttgggCCGGCCTCTTCGCTTTTAAACCCTCTCATAAAATTTTCTTTTATTAGCGTTTTTACCCGCTAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATc                                             <  5:226772/142‑1 (MQ=255)
                                                                                                                gtgCCGGCCTCTTCGCTGTTAAACCCGCTCCTACATTTTTCCTGTAGTTGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATc                                             <  1:117341/144‑1 (MQ=255)
                                                                                                                        ctctTCGCGGGAAAACCCGCTCCTACAAGGTTCCTTTTGGAGGGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATc                                <  3:57329/149‑1 (MQ=39)
                                                                                                                           ttCGCGGTTAAAACCGCTCCAAAAAGGTTCCTTTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAACAGTTCTAATAAAACAGGAAATTATcc                                            <  1:189013/135‑1 (MQ=38)
                                                                                                                           ttCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTTCAGGCATGCTg                             >  7:105174/1‑149 (MQ=255)
                                                                                                                           ttCGCGGGTAAAACCCCGCCGACAAGGTTGCTGTTGGAGCGGGGTTTCCCGCGAAGAGGCCAGCACAGAAAAATAAAAATAACGGATTTTAAAAAAAAAGATCTGAAAAGTTGTAATAAAAAAGGAAATTAtgc                                            >  2:189013/1‑132 (MQ=255)
                                                                                                                                    aaaCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCAc                    >  1:186656/1‑149 (MQ=255)
                                                                                                                                       cccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTa                  >  7:199091/1‑148 (MQ=255)
                                                                                                                                        ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                                       >  7:96856/1‑110 (MQ=255)
                                                                                                                                        ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                                       <  8:96856/110‑1 (MQ=255)
                                                                                                                                          gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct               >  2:100438/1‑148 (MQ=38)
                                                                                                                                               cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGa          >  3:300419/1‑149 (MQ=255)
                                                                                                                                                  cAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCtt                                          <  2:40785/114‑1 (MQ=255)
                                                                                                                                                  cAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCtt                                          >  1:40785/1‑114 (MQ=255)
                                                                                                                                                    cggTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg     <  2:81758/148‑1 (MQ=255)
                                                                                                                                                      gTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgc   <  8:199091/149‑1 (MQ=255)
                                                                                                                                                      gTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgc   >  3:345224/1‑149 (MQ=255)
                                                                                                                                                        tCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttct            <  6:233570/138‑1 (MQ=255)
                                                                                                                                                        tCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttct            >  5:233570/1‑138 (MQ=255)
                                                                                                                                                        tCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgcc  <  4:243710/148‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
GCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGCC  >  NC_002947/955006‑955300

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