Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,786,016 C→T 7.5% R398H (CGC→CAC)  mhpT ← 3‑(3‑hydroxy‑phenyl)propionate transporter MhpT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,786,0160CT 7.5% 152.0 / 2.5 54R398H (CGC→CAC) mhpT3‑(3‑hydroxy‑phenyl)propionate transporter MhpT
Reads supporting (aligned to +/- strand):  ref base C (26/23);  new base T (2/2);  total (29/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 = 7.46e-01

CATGTACCTTTCCAAGGTATCGGGCCGTAACCGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAACCCAGGCGCCCCA  >  NC_002947/3785870‑3786163
                                                                                                                                                  |                                                                                                                                                    
cATGTACCTTTCCAAGGTATCGGGCCGTAACCGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCgg                                                                                                                                                    <  7:213455/149‑1 (MQ=255)
 aTGTACCTTTCCAAGGTATCGGGCCGTAACCGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCggg                                                                                                                                                   >  7:221030/1‑149 (MQ=255)
 aTGTACCTTTCCAAGGTATCGGGCCGTAACCGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCggg                                                                                                                                                   >  2:262907/1‑149 (MQ=255)
           ccAAGGTATCGGGCCGTAACCGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGTGGGTTCGGTGttt                                                                                                                                         >  8:175782/1‑149 (MQ=255)
           ccAAGGTATCGGGCCGTAACCGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGTGGGTTCGGTGttt                                                                                                                                         >  6:136060/1‑149 (MQ=255)
           ccAAGGTATCGGGCCGTAACCGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGttt                                                                                                                                         >  6:217342/1‑149 (MQ=255)
            cAAGGTATCGGGCCGTAACCGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGTGGGTTCGGTGttt                                                                                                                                         <  5:136060/148‑1 (MQ=255)
            cAAGGTATCGGGCCGTAACCGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGTGGGTTCGGTGttt                                                                                                                                         <  7:175782/148‑1 (MQ=255)
            cAAGGTATCGGGCCGTAACCGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGttt                                                                                                                                         >  5:21844/1‑148 (MQ=255)
            cAAGGTATCGGGCCGTAACCGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGttt                                                                                                                                         <  5:365369/148‑1 (MQ=255)
            cAAGGTATCGGGCCGTAACCGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGttt                                                                                                                                         <  6:305489/148‑1 (MQ=255)
            cAAGGTATCGGGCCGTAACCGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGttt                                                                                                                                         >  3:254863/1‑148 (MQ=255)
            cAAGGTATCGGGCCGTAACCGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTc                                                                                                                                        <  2:341363/149‑1 (MQ=255)
             aaGGTATCGGGCCGTAACCGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCa                                                                                                                                       <  4:254863/149‑1 (MQ=255)
              aGGTATCGGGCCGTAACCGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAg                                                                                                                                      >  3:395899/1‑149 (MQ=255)
                gTATCGGGCCGTAACCGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTa                                                                                                                                    <  4:73048/149‑1 (MQ=255)
                              ccGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTgccgcc                                                                                                                      <  8:221030/149‑1 (MQ=255)
                                         tCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAg                                                                                                           >  5:41071/1‑149 (MQ=255)
                                           gCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCa                                                                                                         <  6:41071/149‑1 (MQ=255)
                                            ccGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCAc                                                                                                        >  6:86386/1‑149 (MQ=255)
                                                   ggAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCt                                                                                                 <  5:217342/149‑1 (MQ=255)
                                                   ggAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCt                                                                                                 >  5:186087/1‑149 (MQ=255)
                                                      gCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGa                                                                                              <  4:73168/149‑1 (MQ=255)
                                                          gaggCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCcg                                                                                          >  8:365469/1‑149 (MQ=255)
                                                                  gcgcAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCAt                                                                                   <  4:395899/148‑1 (MQ=255)
                                                                  gcgcAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCAAGCACCAGGCCTGGAGCCGCGGCCAt                                                                                   >  7:51341/1‑148 (MQ=255)
                                                                    gcAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATcag                                                                                >  1:362361/1‑149 (MQ=255)
                                                                      aaGCAGGACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATcagca                                                                              <  8:319333/149‑1 (MQ=255)
                                                                      aaGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCAcc                                                                                                       <  3:167341/124‑1 (MQ=255)
                                                                      aaGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGACTCACCCGGCTCAGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCAcc                                                                                                       >  4:167341/1‑124 (MQ=255)
                                                                       agcagcACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATcagcag                                                                             >  4:151636/1‑149 (MQ=255)
                                                                              cTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑ccccccc                                                                     <  5:86386/149‑1 (MQ=255)
                                                                               tAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCg                                                                    >  4:193645/1‑149 (MQ=255)
                                                                                 gCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCgca                                                                  >  2:211079/1‑149 (MQ=255)
                                                                                  ccccccAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCgcac                                                                 <  5:339301/149‑1 (MQ=255)
                                                                                     cccAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCgcacctg                                                              >  6:86174/1‑149 (MQ=255)
                                                                                        aaGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCgcacctgcac                                                           <  3:151636/149‑1 (MQ=255)
                                                                                         ggctaGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCgcacctgcacc                                                          <  1:321505/145‑1 (MQ=255)
                                                                                              gTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTAGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCGCACCTGCACctgctg                                                     >  1:90335/1‑149 (MQ=255)
                                                                                                   aTGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCCCGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CACCCCCGCACCTGCACctgctactttg                                                >  5:415117/1‑143 (MQ=255)
                                                                                                    tGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCGCACCTGCACCTGCTGCAAg                                                 <  5:86174/147‑1 (MQ=255)
                                                                                                    tGCAGGGTCCTGGGCTGCAGCCTAACCCGGCTCCGGTCGTTGCGTGCTGGTTCGGTGTTTGATTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGACTGGTGCCGCGGCCATAAGCAG‑CCCCCCCGCACCTGCACATGCTGCAAgg                                                <  8:169578/148‑1 (MQ=255)
                                                                                                       aGGGTCCTGGGCTGCAGCCTCACCCGGCTCAGGGCGTTGCGTGAGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCGCACCAGCACCTGCGGCAAGGATCt                                            >  7:293347/1‑149 (MQ=255)
                                                                                                              tGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGgc                                     >  3:348993/1‑149 (MQ=255)
                                                                                                                ggCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGgccg                                   >  2:43863/1‑149 (MQ=255)
                                                                                                                  cTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGgccgc                                  >  5:269354/1‑149 (MQ=255)
                                                                                                                   tGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGgccgc                                  >  1:323102/1‑147 (MQ=255)
                                                                                                                     cAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCa                                <  6:21844/147‑1 (MQ=255)
                                                                                                                         cTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGAc                          >  3:35736/1‑149 (MQ=255)
                                                                                                                          tCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑TCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGAcc                         >  3:355685/1‑149 (MQ=255)
                                                                                                                          tCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGAcc                         >  2:136755/1‑149 (MQ=255)
                                                                                                                             cccGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGc                      <  4:348993/149‑1 (MQ=255)
                                                                                                                               cGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCa                    <  1:43863/149‑1 (MQ=255)
                                                                                                                                     cGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAAc              <  1:102445/149‑1 (MQ=255)
                                                                                                                                               gtgcgGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCACGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCGCACCTGCACCAGCAGCAAGGAt                                              >  4:377567/1‑107 (MQ=255)
                                                                                                                                                 gcgGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCCACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCCG‑CCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAACCCAGgtgcgcg   >  2:137121/1‑145 (MQ=255)
                                                                                                                                                 gcgGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAACCCAGGCGCCcca  >  3:275197/1‑149 (MQ=255)
                                                                                                                                                  cgGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑TCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAACCCAGGCGCCcca  <  4:355685/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CATGTACCTTTCCAAGGTATCGGGCCGTAACCGGGTCACGGTCGCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAG‑CCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAACCCAGGCGCCCCA  >  NC_002947/3785870‑3786163

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