Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,533,133 C→G 6.8% H382D (CAC→GAC) ‡ secA → protein translocase subunit SecA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,533,1330CG 6.8% 183.5 / 3.8 59H382D (CAC→GAC) ‡secAprotein translocase subunit SecA
Reads supporting (aligned to +/- strand):  ref base C (25/30);  new base G (2/2);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

TTGCTGACACACGTCTACGCCGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTC  >  NC_002947/1532989‑1533281
                                                                                                                                                |                                                                                                                                                     
ttGCTGACACACGTCTACGCCGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCAc                                                                                                                                                     >  7:210370/1‑147 (MQ=255)
  gCTGACACACGTCTACGCCGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCa                                                                                                                                                   >  4:228473/1‑147 (MQ=255)
  gCTGACACACGTCTACGCCGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCa                                                                                                                                                   <  3:228473/147‑1 (MQ=255)
   cTGACACACGTCTACGCCGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAgg                                                                                                                                                 <  8:210370/148‑1 (MQ=255)
   cTGACACACGTCTACGCCGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAgg                                                                                                                                                 <  5:260869/148‑1 (MQ=255)
    tGACACACGTCTACGCCGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGc                                                                                                                                                >  8:286388/1‑148 (MQ=255)
     gACACACGTCTACGCCGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGACAAGGCCGGCCAGGTTCTGCTGATCGACGACCCCACCGGCCGTCCCATGCCGGGCCGTCGCCTGACCGAGGGCCTGCCCCAGGCCa                                                                                                                                              >  1:346759/1‑149 (MQ=255)
      acacacGTCTACGCCGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTTCGATGGCCTGGTGCAGGccct                                                                                                                                             >  6:132679/1‑147 (MQ=255)
       cacacGTCTACGCCGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCTTGGTGCAGGCCCTc                                                                                                                                            <  6:177393/149‑1 (MQ=255)
        acacGTCTACGCCGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCAt                                                                                                                                             <  4:110319/147‑1 (MQ=255)
        acacGTCTACGCCGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCAt                                                                                                                                             >  3:110319/1‑147 (MQ=255)
        acacGTCTACGCCGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATc                                                                                                                                            >  5:347299/1‑148 (MQ=255)
          acGTCTACGCCGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCg                                                                                                                                           >  1:282053/1‑147 (MQ=255)
          acGTCTACGCCGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCg                                                                                                                                           >  5:143201/1‑147 (MQ=255)
             tCTACGCCGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTTCGATGGCCTGGTGCAGGCCCTCggacag                                                                                                                                      <  5:132679/149‑6 (MQ=255)
               tACGCCGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCGATCGAAgaaa                                                                                                                                     >  6:216329/1‑145 (MQ=255)
                  ggcGGTCTTCGTGCACACAAGCTGTTCCACCGCAAGGTCGAGTAGATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACAGACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCaa                                                                                                                                    <  2:321105/144‑1 (MQ=255)
                  gCCGGTCTTCGTGCACACAAGCTGTTCAACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACCCACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCAACGAAGCCaaaa                                                                                                                                  >  1:321105/1‑147 (MQ=255)
                    cGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAagag                                                                                                                                >  7:195625/1‑148 (MQ=255)
                    cGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAgaga                                                                                                                               <  6:143201/149‑1 (MQ=255)
                            gTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAg                                                                                                                                                  >  3:129031/1‑122 (MQ=255)
                            gTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAg                                                                                                                                                  <  4:129031/122‑1 (MQ=255)
                            gTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAAc                                                                                                                             <  2:186551/143‑1 (MQ=255)
                            gTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAAc                                                                                                                             >  1:186551/1‑143 (MQ=255)
                               cacacaAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATc                                                                                                                                            <  7:302201/125‑1 (MQ=255)
                               cacacaAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATc                                                                                                                                            >  8:302201/1‑125 (MQ=255)
                                acacaAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATc                                                                                                                    >  2:154876/1‑148 (MQ=255)
                                    aaGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGc                                                                                                               <  8:87871/149‑1 (MQ=255)
                                    aaGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGc                                                                                                               >  1:321485/1‑149 (MQ=255)
                                          ttCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCgaga                                                                                                          <  6:158671/148‑1 (MQ=255)
                                            ccACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGc                                                                                                        <  2:105275/148‑1 (MQ=255)
                                             cACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCa                                                                                                      <  7:286388/149‑1 (MQ=255)
                                                           gTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCTTCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGGGTCTGG‑Ctctcg                                                                                       <  2:135588/149‑5 (MQ=255)
                                                            tACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGGCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGa                                                                                       <  6:119223/148‑1 (MQ=255)
                                                                   tCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTc                                                                               >  1:239886/1‑149 (MQ=255)
                                                                      aGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTc                                                                               >  2:79297/1‑146 (MQ=255)
                                                                        gACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGaa                                                                          <  1:195926/149‑1 (MQ=255)
                                                                                ggTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTcc                                                                  >  8:251016/1‑149 (MQ=255)
                                                                                  ttCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGc                                                                <  2:282053/149‑1 (MQ=255)
                                                                                        tgATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGaccacc                                                                                  >  2:360080/1‑125 (MQ=255)
                                                                                        tgATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGaccacc                                                                                  <  1:360080/125‑1 (MQ=255)
                                                                                         gATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTaca                                                         >  1:176532/1‑149 (MQ=255)
                                                                                           tCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTacac                                                        <  8:195625/148‑1 (MQ=255)
                                                                                              acgaACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAg                                                    <  6:347299/149‑1 (MQ=255)
                                                                                                   cacacCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAg                                                    >  7:145347/1‑144 (MQ=255)
                                                                                                   cacacCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAg                                                    <  8:145347/144‑1 (MQ=255)
                                                                                                    acacCGGCCGTACCATGCCGGGCCGTCGGCTGTCCGAGGGCCGGCACCAGGCCATCGTAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTc                                               <  1:79297/149‑1 (MQ=255)
                                                                                                     cacCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCg                                             <  7:251016/149‑1 (MQ=255)
                                                                                                       cgggCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGc                                           <  2:346759/147‑1 (MQ=255)
                                                                                                         ggCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCAt                                         <  7:75432/149‑1 (MQ=255)
                                                                                                               aCCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCg                                             >  7:115714/1‑139 (MQ=255)
                                                                                                               aCCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCg                                             <  8:115714/139‑1 (MQ=255)
                                                                                                                      cGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTAccgcc                            <  2:176532/149‑1 (MQ=255)
                                                                                                                      cGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTAccgcc                            >  8:270734/1‑149 (MQ=255)
                                                                                                                               gCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGaa                   <  2:321485/149‑1 (MQ=255)
                                                                                                                                     tcGATGGCCTGGTGCAGGCCATCGAAGCCAAAGAGACCCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTc             >  1:200261/2‑149 (MQ=255)
                                                                                                                                      cGAGGGCCTGCACCAGGCCATCGTAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTc             <  5:7661/148‑1 (MQ=255)
                                                                                                                                       gAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGTCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCt                                                  <  5:163131/110‑1 (MQ=255)
                                                                                                                                       gAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGCCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCt                                                  >  6:163131/1‑110 (MQ=255)
                                                                                                                                                cACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAg    <  1:154876/147‑1 (MQ=255)
                                                                                                                                                cACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGt   <  7:270734/148‑1 (MQ=255)
                                                                                                                                                cACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTc  >  7:30361/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TTGCTGACACACGTCTACGCCGGTCTTCGTGCACACAAGCTGTTCCACCGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGG‑CTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTC  >  NC_002947/1532989‑1533281

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