Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 325,161 C→G 5.1% P149A (CCA→GCA)  oprQ → outer‑membrane porin D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947325,1610CG 5.1% 1238.2 / 37.6 363P149A (CCA→GCA) oprQouter‑membrane porin D
Reads supporting (aligned to +/- strand):  ref base C (191/153);  new base G (6/13);  total (197/166)
Fisher's exact test for biased strand distribution p-value = 5.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 372 total aligned reads displayed.

GTGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGC  >  NC_002947/325012‑325310
                                                                                                                                                     |                                                                                                                                                     
   gcggcgTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAg                                                                                                                                                    >  1:5093891/1‑150 (MQ=255)
    cggcgTCGACTTCTTCAAGCCGCACGAAACCACCAAGGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGt                                                                                                                                                   <  1:34318/150‑1 (MQ=255)
     ggcgTCGACTTCTTCAAGCCGAACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGtg                                                                                                                                                  <  2:5792851/150‑1 (MQ=255)
        gTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGGCAGtgctg                                                                                                                                               <  1:4388389/150‑1 (MQ=255)
          cGATTTCTTCAACCCGACCTAAACCACCAATCCCGATGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGGCAGtgctggc                                                                                                                                             <  1:130675/150‑3 (MQ=255)
          cGACTTCTTCAAGCCGCACTAAACCACCAAGGCCGACGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCa                                                                                                                                             <  1:1684162/150‑1 (MQ=255)
          cGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCa                                                                                                                                             <  1:5915993/150‑1 (MQ=255)
            aCTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGGCAGtgctggcat                                                                                                                                           <  2:4323580/150‑5 (MQ=255)
            aCTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGt                                                                                                                                           >  2:3507007/1‑150 (MQ=255)
              ttcttcAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTac                                                                                                                                         >  2:5906262/1‑150 (MQ=255)
              ttcttcAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTac                                                                                                                                         >  2:4932139/1‑150 (MQ=255)
              ttcttcAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTac                                                                                                                                         >  2:6378668/1‑150 (MQ=255)
               tcttcAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGGCAGtgctggcatctg                                                                                                                                        <  1:1424053/150‑8 (MQ=255)
                cttcAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTacga                                                                                                                                       <  1:3133786/150‑1 (MQ=255)
                 ttcAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTacgac                                                                                                                                      <  1:3957111/150‑1 (MQ=255)
                 ttcAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTacgac                                                                                                                                      >  2:1138091/1‑150 (MQ=255)
                 ttcAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTacgac                                                                                                                                      >  2:5361241/1‑150 (MQ=255)
                  tcAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTacgacg                                                                                                                                     >  2:4308101/1‑150 (MQ=255)
                   cAAGCCGAACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTGcgacga                                                                                                                                    <  1:3474712/150‑1 (MQ=255)
                        cGCACGAAACCACCAATGCCGAAGGCAACGCCAGCCCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCAc                                                                                                                               <  1:218054/150‑1 (MQ=255)
                         gCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACg                                                                                                                              >  2:81461/1‑150 (MQ=255)
                               acccacaaTTGCCGATGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTctgct                                                                                                                        <  1:4600786/148‑1 (MQ=255)
                                  caccaATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGcc                                                                                                                     >  1:5202305/1‑150 (MQ=255)
                                  caccaATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCGGCTGTGAAGGTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGcc                                                                                                                     <  2:3876056/150‑1 (MQ=255)
                                   accaATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCa                                                                                                                    >  1:3228346/1‑150 (MQ=255)
                                    ccaATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAg                                                                                                                   >  1:4361052/1‑150 (MQ=255)
                                      aaTGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGaa                                                                                                                 <  2:4144994/150‑1 (MQ=255)
                                       aTGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGaaa                                                                                                                <  1:2251318/150‑1 (MQ=255)
                                       aTGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGaaa                                                                                                                <  1:2610958/150‑1 (MQ=255)
                                         gCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGc                                                                                                              <  1:3290278/150‑1 (MQ=255)
                                           cGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCtt                                                                                                            <  2:2982364/150‑1 (MQ=255)
                                            gAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTc                                                                                                           >  2:649725/1‑150 (MQ=255)
                                            gAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTc                                                                                                           <  1:3052714/150‑1 (MQ=255)
                                               ggCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCAcc                                                                                                        <  1:6534259/150‑1 (MQ=255)
                                                   aCGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCa                                                                                                    <  2:235741/150‑1 (MQ=255)
                                                         gCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATg                                                                                                                                  <  1:505435/114‑1 (MQ=255)
                                                         gCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATg                                                                                                                                  >  2:505435/1‑114 (MQ=255)
                                                          cTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGAt                                                                                             <  1:567601/150‑1 (MQ=255)
                                                                gcACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTc                                                                                       <  2:3766636/150‑1 (MQ=255)
                                                                gcACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTc                                                                                       >  1:4324400/1‑150 (MQ=255)
                                                                gcACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTc                                                                                       <  2:1050225/150‑1 (MQ=255)
                                                                    aaCCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCaaa                                                                                   <  1:6270381/150‑1 (MQ=255)
                                                                    aaCCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCaaa                                                                                   <  2:717583/150‑1 (MQ=255)
                                                                    aaCCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCaaa                                                                                   <  1:6752676/150‑1 (MQ=255)
                                                                              gtggtggtgCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAgg                                                                         >  2:1192410/1‑150 (MQ=255)
                                                                               tggtggtgCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAggg                                                                        <  1:3563326/150‑1 (MQ=255)
                                                                                   ggtgCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTg                                                                    >  1:4213958/1‑150 (MQ=255)
                                                                                    gtgCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTgg                                                                   >  2:3999716/1‑150 (MQ=255)
                                                                                      gctgctGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGGCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGaa                                                                 <  2:4539473/150‑1 (MQ=255)
                                                                                      gctgctGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGaa                                                                 <  1:1061348/150‑1 (MQ=255)
                                                                                      gctgctGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGaa                                                                 <  2:2434338/150‑1 (MQ=255)
                                                                                      gctgctGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGaa                                                                 <  1:6145404/150‑1 (MQ=255)
                                                                                        tgctgTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACt                                                               >  1:781797/1‑150 (MQ=255)
                                                                                         gctgTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTg                                                              <  2:5256194/150‑1 (MQ=255)
                                                                                         gctgTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTg                                                              >  2:2766253/1‑150 (MQ=255)
                                                                                          ctgTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGa                                                             <  2:5529090/150‑1 (MQ=255)
                                                                                          ctgTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGa                                                             <  1:239991/150‑1 (MQ=255)
                                                                                           tgtgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGaa                                                            >  1:961638/1‑150 (MQ=255)
                                                                                           tgtgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGaa                                                            >  1:4734465/1‑150 (MQ=255)
                                                                                            gtgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAAc                                                           >  1:2936998/1‑150 (MQ=255)
                                                                                            gtgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAAc                                                           <  1:1374967/150‑1 (MQ=255)
                                                                                                 gTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACgccgg                                                      >  1:6473471/1‑150 (MQ=255)
                                                                                                  ttCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACgccggc                                                     <  2:6283595/150‑1 (MQ=255)
                                                                                                       cATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCtt                                                >  1:1094547/1‑150 (MQ=255)
                                                                                                        aTTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTc                                               >  2:1399994/1‑150 (MQ=255)
                                                                                                        aTTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACAGAACGCCGGCCGCTTc                                               >  2:6661225/1‑150 (MQ=255)
                                                                                                          ttCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCAc                                             <  1:2835007/150‑1 (MQ=255)
                                                                                                          tctggtcacCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCAc                                             >  1:3713249/7‑150 (MQ=255)
                                                                                                          tctggtcacCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCAc                                             >  1:385328/7‑150 (MQ=255)
                                                                                                           tCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCAcc                                            <  1:6703657/150‑1 (MQ=255)
                                                                                                            ccAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCAccc                                           <  2:2645948/150‑1 (MQ=255)
                                                                                                             cAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCa                                          >  1:1595682/1‑150 (MQ=255)
                                                                                                             cAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCa                                          >  1:1832097/1‑150 (MQ=255)
                                                                                                             cAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCa                                          >  1:2681772/1‑150 (MQ=255)
                                                                                                             cAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCa                                          >  2:5764966/1‑150 (MQ=255)
                                                                                                             cAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCa                                          >  2:1351261/1‑150 (MQ=255)
                                                                                                             cAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCa                                          <  1:5620758/150‑1 (MQ=255)
                                                                                                              aaCACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAg                                         <  2:2491047/150‑1 (MQ=255)
                                                                                                              aaCACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAg                                         >  1:2320924/1‑150 (MQ=255)
                                                                                                                   cGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAgc                                    <  2:4385622/150‑1 (MQ=255)
                                                                                                                   cGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAgc                                    >  2:2258480/1‑150 (MQ=255)
                                                                                                                         gtaGTACGGTGACCAGATGCCAGCCCTGCCAGTGCGGCAGTACGACGATGCACGTCTGCTGCCAGAACGCTGCACCGGCACCCTGATCACCTCCACAGAGATCCAGGGGCTGGAACTGAACGCCGGGCGCTTCACGCGGGACGCGCGCag                              >  2:187247/3‑149 (MQ=255)
                                                                                                                         gAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCaa                              >  2:489159/1‑150 (MQ=255)
                                                                                                                         gAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCaa                              >  2:366526/1‑150 (MQ=255)
                                                                                                                             tACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAgagc                          >  1:5316511/1‑150 (MQ=255)
                                                                                                                                  tGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAACGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTACACCCAGGAAGCGCGCAAgagcgccga                     <  2:3545027/150‑1 (MQ=255)
                                                                                                                                  tGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAgagcgccga                     <  1:6046104/150‑1 (MQ=255)
                                                                                                                                     ccAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCaaa                                                                                   >  2:3156979/1‑85 (MQ=255)
                                                                                                                                      cAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAgagcgccgagcgc                 >  1:489138/1‑150 (MQ=255)
                                                                                                                                      cAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAgagcgccgagcgc                 >  1:6629307/1‑150 (MQ=255)
                                                                                                                                        gATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGagcgccgagcgccg               >  1:1997910/1‑150 (MQ=255)
                                                                                                                                           gCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGa            >  2:1443513/1‑150 (MQ=255)
                                                                                                                                           gCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGa            >  1:84470/1‑150 (MQ=255)
                                                                                                                                            ccAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAt           <  2:3253500/150‑1 (MQ=255)
                                                                                                                                              aGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAtgg         >  2:6225210/1‑150 (MQ=255)
                                                                                                                                              aGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAtgg         >  2:2022301/1‑150 (MQ=255)
                                                                                                                                               gCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAtggt        >  2:555388/1‑150 (MQ=255)
                                                                                                                                                    gCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAtggtggtgg   <  2:1675233/150‑1 (MQ=255)
                                                                                                                                                     ccAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGc  <  2:1275855/150‑1 (MQ=255)
                                                                                                                                                     ccAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGc  >  2:6041780/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGC  >  NC_002947/325012‑325310

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