Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,208,664 G→T 9.5% R107L (CGC→CTC)  mlrA → transcriptional activator of csgD and csgBA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,208,6640GT 9.5% 850.2 / 44.1 231R107L (CGC→CTC) mlrAtranscriptional activator of csgD and csgBA
Reads supporting (aligned to +/- strand):  ref base G (125/84);  new base T (13/9);  total (138/93)
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.92e-01

Only 100 of 238 total aligned reads displayed.

CGACAACGGCGTGCAGGTCAGCAAAGTTAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGACGGCGTACTGATCAACTACATCGCCATTTGTCTGGCTTCG  >  CP009273/2208515‑2208812
                                                                                                                                                     |                                                                                                                                                    
cGACAACGGCGTGCAGGTCAGCAAAGTTAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAgcg                                                                                                                                                      >  2:1690877/1‑150 (MQ=255)
    aaCGGCGTGCAGGTCAGCAAAGTTAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACAAAGCCTGCGAACGTGGATCAAAGAGCGCGGt                                                                                                                                                  >  1:1356560/1‑150 (MQ=255)
      cGGCGTGCAGGTCAGCAAAGTTAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCa                                                                                                                                                >  1:3553940/1‑150 (MQ=255)
      cGGCGTGCAGGTCAGCAAAGTTAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCa                                                                                                                                                <  1:173578/150‑1 (MQ=255)
      cGGCGTGCAGGTCAGCAAAGTTAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCa                                                                                                                                                <  2:2304539/150‑1 (MQ=255)
       ggCGTGCAGGTCAGCAAAGTTAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAg                                                                                                                                               >  1:2560978/1‑150 (MQ=255)
           tGCAGGTCAGCAAAGTTAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGAtt                                                                                                                                           <  1:389033/150‑1 (MQ=255)
           tGCAGGTCAGCAAAGTTAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGAtt                                                                                                                                           >  2:1208121/1‑150 (MQ=255)
           tGCAGGTCAGCAAAGTTAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGAtt                                                                                                                                           >  1:1233346/1‑150 (MQ=255)
           tGCAGGTCAGCAAAGTTAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGAtt                                                                                                                                           <  1:1268809/150‑1 (MQ=255)
            gCAGGTCAGCAAAGTTAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTa                                                                                                                                          <  2:3441737/150‑1 (MQ=255)
               ggTCAGCAAAGTTAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTAccc                                                                                                                                       >  2:256845/1‑150 (MQ=255)
                gTCAGCAAAGTTAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTAcccc                                                                                                                                      >  2:1511280/1‑150 (MQ=255)
                 tCAGCAAAGTTAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCg                                                                                                                                     <  2:81819/150‑1 (MQ=255)
                 tCAGCAAAGTTAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCg                                                                                                                                     >  2:2203259/1‑150 (MQ=255)
                        aGTTAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCTCGGTCAGGATTACCCCGCCCAGac                                                                                                                              >  1:2670746/1‑150 (MQ=255)
                           tAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACt                                                                                                                           >  2:3158103/1‑150 (MQ=255)
                           tAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACt                                                                                                                           >  2:3158098/1‑150 (MQ=255)
                           tAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACt                                                                                                                           >  2:2760504/1‑150 (MQ=255)
                            aaaaTGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTc                                                                                                                          <  1:1951705/150‑1 (MQ=255)
                            aaaaTGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTc                                                                                                                          >  1:2198366/1‑150 (MQ=255)
                                tgctgcTCAGTAATGAAAATGGTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCTCGGTCAGGATTACCCCGCCCAGACACTcacca                                                                                                                      >  2:3905807/1‑150 (MQ=255)
                                  ctgctCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACcaca                                                                                                                    >  2:1532843/1‑150 (MQ=255)
                                   tgctCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCTCGGTCAGGATTACCCCGCCCAGACACTCACcacac                                                                                                                   >  1:1164501/1‑150 (MQ=255)
                                    gctCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACcacaca                                                                                                                  >  2:2632923/1‑150 (MQ=255)
                                     ctCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCTCGGTCAGGATTACCCCGCCCAGACACTCACCACACAt                                                                                                                 <  1:469866/150‑1 (MQ=255)
                                      tCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCTCGGTCAGGATTACCCCGCCCAGACACTCACCACACATc                                                                                                                <  2:396343/150‑1 (MQ=255)
                                         gTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGt                                                                                                             >  2:3105642/1‑150 (MQ=255)
                                         gTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGt                                                                                                             >  1:474120/1‑150 (MQ=255)
                                          tAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCTGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCTGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGtt                                                                                                            >  1:1625354/1‑150 (MQ=255)
                                             tGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTAt                                                                                                         <  2:489326/150‑1 (MQ=255)
                                                aaaTGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTcc                                                                                                      <  2:2117140/150‑1 (MQ=255)
                                                aaaTGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTcc                                                                                                      >  2:2557240/1‑150 (MQ=255)
                                                   tGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCTCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCctct                                                                                                   >  2:3737783/1‑150 (MQ=255)
                                                   tGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCctct                                                                                                   <  2:2074803/150‑1 (MQ=255)
                                                    gTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTg                                                                                                  >  1:1023129/1‑150 (MQ=255)
                                                      tGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTgcg                                                                                                <  1:2160945/150‑1 (MQ=255)
                                                      tGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTgcg                                                                                                >  1:4108550/1‑150 (MQ=255)
                                                         tgtgCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCTCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCg                                                                                             >  2:3203168/1‑150 (MQ=255)
                                                          gtgCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGa                                                                                            <  2:1437995/150‑1 (MQ=255)
                                                          gtgCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGa                                                                                            <  1:1817775/150‑1 (MQ=255)
                                                          gtgCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGa                                                                                            >  1:854067/1‑150 (MQ=255)
                                                            gCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCTCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACg                                                                                          <  2:246122/150‑1 (MQ=255)
                                                             cAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACgg                                                                                         >  1:3094362/1‑150 (MQ=255)
                                                              aGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGc                                                                                        >  1:1527222/1‑150 (MQ=255)
                                                              aGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGc                                                                                        >  1:3587486/1‑150 (MQ=255)
                                                               gAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCt                                                                                       >  1:533900/1‑150 (MQ=255)
                                                               gAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCt                                                                                       >  1:1202378/1‑150 (MQ=255)
                                                               gAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCt                                                                                       >  2:1985779/1‑150 (MQ=255)
                                                               gAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCt                                                                                       >  1:2105298/1‑150 (MQ=255)
                                                                aaCGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCtt                                                                                      >  2:585801/1‑150 (MQ=255)
                                                                 aCGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTc                                                                                     <  1:1872475/150‑1 (MQ=255)
                                                                  cGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCa                                                                                    >  1:2332158/1‑150 (MQ=255)
                                                                  cGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCa                                                                                    <  2:3479656/150‑1 (MQ=255)
                                                                     cTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTg                                                                                 >  2:2908736/1‑150 (MQ=255)
                                                                     cTGGCGCGATCAGCAAGAAACAGTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTg                                                                                 <  1:2028885/150‑1 (MQ=255)
                                                                      tGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGc                                                                                >  1:4007119/1‑150 (MQ=255)
                                                                             aTCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGAGCCaacaac                                                                         >  2:1319457/1‑150 (MQ=255)
                                                                               cAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCg                                                                       >  2:726710/1‑150 (MQ=255)
                                                                                aGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGa                                                                      >  1:2731863/1‑150 (MQ=255)
                                                                                aGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGa                                                                      >  2:2658950/1‑150 (MQ=255)
                                                                                 gCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGAc                                                                     >  1:3384523/1‑150 (MQ=255)
                                                                                   aaGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGActc                                                                   >  2:1618475/1‑150 (MQ=255)
                                                                                    aGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCTCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGActca                                                                  <  1:2405067/150‑2 (MQ=255)
                                                                                      aaaCATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTcc                                                                >  1:1529659/1‑150 (MQ=255)
                                                                                       aaCATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCa                                                               <  2:4060306/150‑1 (MQ=255)
                                                                                        aCATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAg                                                              >  1:3060258/1‑150 (MQ=255)
                                                                                         cATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAgg                                                             <  1:802480/150‑1 (MQ=255)
                                                                                               tGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCgctgc                                                       >  1:1042124/1‑150 (MQ=255)
                                                                                               tGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCgctgc                                                       >  1:3367539/1‑150 (MQ=255)
                                                                                                gACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCgctgct                                                      >  2:2279367/1‑150 (MQ=255)
                                                                                                 aCTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCTCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGATTCAGTGCCAACAACCGACTCTCCAGGCGctgctg                                                     <  2:1114382/150‑1 (MQ=255)
                                                                                                    tACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCg                                                  >  2:3555340/1‑150 (MQ=255)
                                                                                                      ccTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGAt                                                >  2:4177102/1‑150 (MQ=255)
                                                                                                       cTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATc                                               >  2:113989/1‑150 (MQ=255)
                                                                                                       cTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATc                                               <  1:427606/150‑1 (MQ=255)
                                                                                                        tGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATcc                                              <  1:2558597/150‑1 (MQ=255)
                                                                                                          cAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTc                                            >  2:172095/1‑150 (MQ=255)
                                                                                                          cAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTc                                            <  1:158087/150‑1 (MQ=255)
                                                                                                          cAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTc                                            <  1:1751987/150‑1 (MQ=255)
                                                                                                          cAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTc                                            <  2:1253317/150‑1 (MQ=255)
                                                                                                            aaGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGa                                          >  2:4132297/1‑150 (MQ=255)
                                                                                                               cGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGACgg                                       >  2:281772/1‑150 (MQ=255)
                                                                                                                  cAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGACGGCGt                                    >  2:920494/1‑150 (MQ=255)
                                                                                                                   aaTCTACATAGCCTGCGAACGTGGATCAAAGAGCTCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGACGGCGTa                                   <  1:3257416/150‑1 (MQ=255)
                                                                                                                      cTACATAGCCTGCGAACGTGGATCACCGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACACCCGACTCTCCAGGCGCTGCTGGCGATCCTCGACGGCGTACTg                                <  2:1681272/150‑1 (MQ=255)
                                                                                                                        aCATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGACGGCGTACTGAt                              >  2:31188/1‑150 (MQ=255)
                                                                                                                            aGCCTGCGAACGTGTATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCCCACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGACGGCGTACTGATCAAc                          <  2:758725/150‑1 (MQ=255)
                                                                                                                              ccTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGACGGCGTACTGATCAACTa                        >  2:1792351/1‑150 (MQ=255)
                                                                                                                               cTGCGAACGTGGATCAAAGAGCTCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGACGGCGTACTGATCAACTAc                       >  2:1868513/1‑150 (MQ=255)
                                                                                                                                  cGAACGTGGATCAAAGCGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGACGGCGTACTGATCAACTACATc                    <  1:3924970/150‑1 (MQ=255)
                                                                                                                                   gAACGTGGATCAAAGAGCTCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGACGGCGTACTGATCAACTACATCg                   >  2:387080/1‑150 (MQ=255)
                                                                                                                                      cGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGACGGCGTACTGATCAACTACATCGCCa                <  1:1034690/150‑1 (MQ=255)
                                                                                                                                         ggATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGACGGCGTACTGATCAACTACATCGCCAttt             >  1:3752708/1‑150 (MQ=255)
                                                                                                                                             cAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGACGNCGTACTGATCAACTACATCGCCATTTGTCt         <  1:3858519/150‑1 (MQ=255)
                                                                                                                                              aaaGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGACGGCGTACTGATCAACTACATCGCCATTTGTCTg        >  2:3150504/1‑150 (MQ=255)
                                                                                                                                               aaGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGACGGCGTACTGATCAACTACATCGCCATTTGTCTgg       <  1:2990790/150‑1 (MQ=255)
                                                                                                                                                 gagCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACGCTCCAGGCGCTGCTGGCGATCCTCGACGGCGTACTGATCAACTACATCGCCAGTTGTCTGGCt     >  2:142858/1‑150 (MQ=255)
                                                                                                                                                  agCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGACGGCGTACTGATCAACTACATCGCCATTTGTCTGGCtt    <  1:2060998/150‑1 (MQ=255)
                                                                                                                                                   gcgcGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGACGGCGTACTGATCAACTACATCGCCATTTGTCTGGCTTc   <  2:3559241/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGACAACGGCGTGCAGGTCAGCAAAGTTAAAATGCTGCTCAGTAATGAAAATGTTGATGTGCAGAACGGCTGGCGCGATCAGCAAGAAACATTACTGACTTACCTGCAAAGCGGCAATCTACATAGCCTGCGAACGTGGATCAAAGAGCGCGGTCAGGATTACCCCGCCCAGACACTCACCACACATCTGTTTATTCCTCTGCGCCGACGGCTTCAGTGCCAACAACCGACTCTCCAGGCGCTGCTGGCGATCCTCGACGGCGTACTGATCAACTACATCGCCATTTGTCTGGCTTCG  >  CP009273/2208515‑2208812

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: