Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,858,044 G→C 8.0% R32P (CGC→CCC)  PP_4273 → MerR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,858,0440GC 8.0% 239.0 / 7.9 74R32P (CGC→CCC) PP_4273MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (36/32);  new base C (3/3);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGAC  >  NC_002947/4857900‑4858192
                                                                                                                                                |                                                                                                                                                    
cATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGc                                                                                                                                                  <  5:192458/149‑1 (MQ=255)
     cgcgCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGCGCCCGGCCaagg                                                                                                                                             <  8:260985/149‑4 (MQ=255)
      gcgcAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGt                                                                                                                                            <  7:91773/149‑1 (MQ=255)
       cgcAGCATAGCAGATGACGTCAACTTCAAGTCAAGGGGCCGAAATGGACTTTGCCGACGTCGCCCACCGAACTGGCGTACCTGCATCTACACTGCGTTACTACGAGAAGAATGGCTTGCGCAAGTCCTTGGCCGCGCCCGGCCaaggac                                                                                                                                           <  6:308877/149‑6 (MQ=255)
         cagcaTAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGCGCCCGGCCaaggactt                                                                                                                                         <  6:324769/149‑8 (MQ=255)
             aTAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGt                                                                                                                                     >  7:181971/1‑149 (MQ=255)
              ttgAAGTTGACTTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGGCa                                                                                                                                       >  6:55569/3‑146 (MQ=255)
              ttgAAGTTGACTTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACa                                                                                                                                       <  5:55569/144‑1 (MQ=255)
                gAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGttt                                                                                                                                   <  4:144576/148‑1 (MQ=255)
                gAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTg                                                                                                                                  <  6:300974/149‑1 (MQ=255)
                    ttGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGc                                                                                                                                                  >  8:209430/1‑129 (MQ=255)
                    ttGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGc                                                                                                                                                  <  7:209430/129‑1 (MQ=255)
                          tCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGAtgt                                                                                                                        >  1:145494/1‑149 (MQ=255)
                          tCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGAtgt                                                                                                                        <  2:145494/149‑1 (MQ=255)
                            aaCTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCa                                                                                                                                                >  4:311984/1‑123 (MQ=255)
                            aaCTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCa                                                                                                                                                <  3:311984/123‑1 (MQ=255)
                            aaCTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGcc                                                                                                                             <  1:98758/142‑1 (MQ=255)
                            aaCTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGcc                                                                                                                             >  2:98758/1‑142 (MQ=255)
                                  aagtcaagGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGAc                                                                                                                 <  6:13102/148‑1 (MQ=255)
                                           gCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGgcgc                                                                                                       <  8:181971/149‑1 (MQ=255)
                                             ccGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCt                                                                                                      >  3:117661/1‑148 (MQ=255)
                                              cGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCCGTTGGCGCTGa                                                                                                    >  2:317108/1‑149 (MQ=255)
                                               gAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGAt                                                                                                   >  8:257856/1‑149 (MQ=255)
                                                    ggACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGccc                                                                                              >  3:194277/1‑149 (MQ=255)
                                                          tGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTgggg                                                                                         >  8:40398/1‑148 (MQ=255)
                                                            ccGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAg                                                                                      >  8:198696/1‑149 (MQ=255)
                                                             cGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAgg                                                                                     <  7:40398/149‑1 (MQ=255)
                                                               aCGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGt                                                                                                             >  8:189341/1‑123 (MQ=255)
                                                                    gCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGt                                                                                                             >  8:212274/1‑118 (MQ=255)
                                                                    gCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGt                                                                                                             <  7:212274/118‑1 (MQ=255)
                                                                    gCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAg                                                                                      >  4:213605/1‑141 (MQ=255)
                                                                    gCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAg                                                                                      <  3:213605/141‑1 (MQ=255)
                                                                      cAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTg                                                                                            <  1:150478/133‑1 (MQ=255)
                                                                      cAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTg                                                                                            >  2:150478/1‑133 (MQ=255)
                                                                        aaCGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGttt                                                                          >  3:268679/1‑149 (MQ=255)
                                                                                  cGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCg                                                                                <  6:82946/133‑1 (MQ=255)
                                                                                  cGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCg                                                                                >  5:82946/1‑133 (MQ=255)
                                                                                       cTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTg                                                           >  5:128326/1‑149 (MQ=255)
                                                                                        tGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTgg                                                                   <  4:95576/140‑1 (MQ=255)
                                                                                        tGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTgg                                                                   >  3:95576/1‑140 (MQ=255)
                                                                                          cATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTgggg                                                        >  3:239132/1‑149 (MQ=255)
                                                                                          cATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTgg                                                          <  4:117661/147‑1 (MQ=255)
                                                                                           aTCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTggggg                                                       <  4:268679/149‑1 (MQ=255)
                                                                                            tCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCg                                                                                                                            <  4:70907/79‑1 (MQ=255)
                                                                                            tCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCg                                                                                                                            >  3:70907/1‑79 (MQ=255)
                                                                                              gACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCg                                                     <  4:239132/148‑1 (MQ=255)
                                                                                                cacTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATg                                                  <  4:194277/149‑1 (MQ=255)
                                                                                                 acTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGc                                                 >  1:283235/1‑149 (MQ=255)
                                                                                                      gTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCtggtg                                            <  6:179963/149‑1 (MQ=255)
                                                                                                             cGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTg                                      >  6:129627/1‑148 (MQ=255)
                                                                                                              gagaAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTg                                      >  5:240409/1‑147 (MQ=255)
                                                                                                              gagaAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCa                                    <  3:119813/149‑1 (MQ=255)
                                                                                                               tgaagaaAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGTTCGCGCTGGGGCAGGCTGCGGGGGTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGAACTGCAg                                   >  5:254697/2‑149 (MQ=255)
                                                                                                               agaagaAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAg                                   <  6:240409/149‑1 (MQ=255)
                                                                                                                gaagaaAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATg                                                  <  4:232713/133‑1 (MQ=255)
                                                                                                                gaagaaAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATg                                                  >  3:232713/1‑133 (MQ=255)
                                                                                                                   gaaAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTc                                                                        >  5:192445/1‑108 (MQ=255)
                                                                                                                   gaaAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTc                                                                        <  6:192445/108‑1 (MQ=255)
                                                                                                                    aaaGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGa                              >  1:320855/1‑149 (MQ=255)
                                                                                                                      aGGCTTGCTCAAGTCCTTGGCCGGGCGGGGCCAGCGTCGTCATTTTGCAGCCTATTTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGATTTGGGGTCGATGCTGGTGGACCTGCAGGTCGAcc                            <  7:198696/149‑1 (MQ=255)
                                                                                                                      aGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGCCCGGTTGGAGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGAcc                            >  2:295269/1‑149 (MQ=255)
                                                                                                                        gCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCgg                          >  8:263868/1‑149 (MQ=255)
                                                                                                                         cTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGc                         <  2:283235/149‑1 (MQ=255)
                                                                                                                            gCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCgg                          >  8:257991/1‑145 (MQ=255)
                                                                                                                              tCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATg                    >  7:129014/1‑149 (MQ=255)
                                                                                                                              tCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATg                    >  8:225717/1‑149 (MQ=255)
                                                                                                                                   acgctGGCCGCGCCCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATTCTGGTGGACCTGCAGGTCGACCAGCAGAAGCAGAt               >  4:238501/5‑149 (MQ=255)
                                                                                                                                    cgctGGCCGCGCCCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACTGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATg              >  2:184902/4‑149 (MQ=255)
                                                                                                                                     gctGGCCGCGCCCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATgg             >  7:150857/3‑149 (MQ=255)
                                                                                                                                     cTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATgg             <  7:257856/149‑1 (MQ=255)
                                                                                                                                         gCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCa          <  5:129627/148‑1 (MQ=255)
                                                                                                                                          ccGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAg        >  4:111865/1‑149 (MQ=255)
                                                                                                                                           cGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAg                                   >  6:131154/1‑121 (MQ=255)
                                                                                                                                           cGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAg                                   <  5:131154/121‑1 (MQ=255)
                                                                                                                                              gcgcGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCgg    <  8:129014/149‑1 (MQ=255)
                                                                                                                                                gcgGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGAc  >  3:121176/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGAC  >  NC_002947/4857900‑4858192

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