Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,970,031 C→A 10.9% S25* (TCA→TAA)  galR → galactose‑inducible d‑galactose regulon transcriptional repressor; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,970,0310CA10.9% 842.5 / 48.1 234S25* (TCA→TAA) galRgalactose‑inducible d‑galactose regulon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (102/106);  new base A (10/16);  total (112/122)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 236 total aligned reads displayed.

CACGCCACCCCTTGAACCAACGGGCGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTC  >  CP009273/2969882‑2970179
                                                                                                                                                     |                                                                                                                                                    
cacGCCACCCCTTGAACCAACGGGCGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTc                                                                                                                                                      >  1:3196848/1‑150 (MQ=255)
  cGCCACCCCTTGAACCAACGGGCGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCAc                                                                                                                                                    >  2:4268695/1‑150 (MQ=255)
  cGCCACCCCTTGAACCAACGGGCGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCAc                                                                                                                                                    <  1:1258812/150‑1 (MQ=255)
   gCCACCCCTTGAACCAACGGGCGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTAAcc                                                                                                                                                   <  1:1211835/150‑1 (MQ=255)
   gCCACCCCTTGAACCAACGGGCGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTAAcc                                                                                                                                                   <  2:1192637/150‑1 (MQ=255)
     cACCCCTTGAACCAACGGGCGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCa                                                                                                                                                 <  1:603206/150‑1 (MQ=255)
     cACCCCTTGAACCAACGGGCGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCa                                                                                                                                                 <  2:4468663/150‑1 (MQ=255)
      aCCCCTTGAACCCACGGGCGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCaa                                                                                                                                                <  2:1226767/150‑1 (MQ=255)
      aCCCCTTGAACCAACGGGCGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCaa                                                                                                                                                <  2:1226822/150‑1 (MQ=255)
      aCCCCTTGAACCAACGGGCGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCaa                                                                                                                                                <  2:2521494/150‑1 (MQ=255)
           ttGAACCAACGGGCGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCATCCAAAGCCa                                                                                                                                           <  2:417906/150‑1 (MQ=255)
           ttGAACCAACGGGCGTTATACCTAAAACAGAAAGAATGTAAGCGTTTACCCAATAAGGTATTTTCATGGAGACCATAAAGTATGTAGCCCGACTGGCAGGCGTTTCAGTCCCCACCGTTTCCCGCGTCAATAATAAATCACCCAAAGCCa                                                                                                                                           >  2:3613761/1‑150 (MQ=255)
              aaCCAACGGGCGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCg                                                                                                                                        >  2:1818133/1‑150 (MQ=255)
                 cAACGGGCGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAg                                                                                                                                     <  2:4295944/150‑1 (MQ=255)
                 cAACGGGCGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAg                                                                                                                                     >  2:4564967/1‑150 (MQ=255)
                    cGGGCGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCtt                                                                                                                                  >  1:3340670/1‑150 (MQ=255)
                      ggCGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTcc                                                                                                                                <  2:182006/150‑1 (MQ=255)
                        cGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCg                                                                                                                              >  2:3355446/1‑150 (MQ=255)
                         gTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCgg                                                                                                                             >  1:1321529/1‑150 (MQ=255)
                         gTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCgg                                                                                                                             >  1:1282993/1‑150 (MQ=255)
                         gTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCgg                                                                                                                             >  1:1282980/1‑150 (MQ=255)
                           tttCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCggct                                                                                                                           <  2:2185881/150‑1 (MQ=255)
                            ttCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCggctg                                                                                                                          <  1:2026696/150‑1 (MQ=255)
                            ttCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCggctg                                                                                                                          <  1:1054072/150‑1 (MQ=255)
                               cGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCggctggct                                                                                                                       <  1:4166349/150‑1 (MQ=255)
                                gTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGgctggctg                                                                                                                      <  2:1856914/150‑1 (MQ=255)
                                  aaCACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCtgtg                                                                                                                    >  2:1639007/1‑150 (MQ=255)
                                  aaCACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCtgtg                                                                                                                    <  1:1116246/150‑1 (MQ=255)
                                   acacTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGc                                                                                                                   >  1:3398592/1‑150 (MQ=255)
                                    cacTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCa                                                                                                                  <  1:1182185/150‑1 (MQ=255)
                                     acTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTAACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCAt                                                                                                                 >  1:756768/1‑150 (MQ=255)
                                      cTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATa                                                                                                                >  1:1004656/1‑150 (MQ=255)
                                        gAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGt                                                                                                              <  2:4088559/150‑1 (MQ=255)
                                            gAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCaa                                                                                                          >  2:1649276/1‑150 (MQ=255)
                                             aaTGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAAt                                                                                                         <  1:2210137/150‑1 (MQ=255)
                                             aaTGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAAt                                                                                                         <  1:2210142/150‑1 (MQ=255)
                                             aaTGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAAt                                                                                                         <  1:2210149/150‑1 (MQ=255)
                                             aaTGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAAt                                                                                                         >  1:1751326/1‑150 (MQ=255)
                                              aTGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATg                                                                                                        >  1:2950132/1‑150 (MQ=255)
                                              aTGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATg                                                                                                        <  1:2649635/150‑1 (MQ=255)
                                               tGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATgg                                                                                                       <  2:1550656/150‑1 (MQ=255)
                                                gTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTGTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGa                                                                                                      >  1:2583539/1‑150 (MQ=255)
                                                 tAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAg                                                                                                     >  2:819763/1‑150 (MQ=255)
                                                  aaGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGt                                                                                                    >  1:295684/1‑150 (MQ=255)
                                                  aaGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGt                                                                                                    <  1:438281/150‑1 (MQ=255)
                                                    gCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGtct                                                                                                  <  1:141182/150‑1 (MQ=255)
                                                      gTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTctct                                                                                                <  2:4069082/150‑1 (MQ=255)
                                                      gTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTctct                                                                                                >  1:3800331/1‑150 (MQ=255)
                                                       tttACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCtt                                                                                               <  2:3143190/150‑1 (MQ=255)
                                                        ttACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTa                                                                                              <  2:2179536/150‑1 (MQ=255)
                                                         tACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAg                                                                                             <  2:3516744/150‑1 (MQ=255)
                                                         tACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAg                                                                                             >  1:2013162/1‑150 (MQ=255)
                                                               cTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCa                                                                                       >  1:3969932/1‑150 (MQ=255)
                                                                 aaGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCAcc                                                                                     >  2:1392850/1‑150 (MQ=255)
                                                                  aGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCAccc                                                                                    >  2:1017253/1‑150 (MQ=255)
                                                                  aGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCAccc                                                                                    <  1:1361223/150‑1 (MQ=255)
                                                                        tttCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGaacgc                                                                              >  2:2320709/1‑150 (MQ=255)
                                                                        tttCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGaacgc                                                                              <  1:3516263/150‑1 (MQ=255)
                                                                         ttCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGaacgcc                                                                             >  2:4357560/1‑150 (MQ=255)
                                                                           cATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGaacgccaa                                                                           >  2:1919682/1‑150 (MQ=255)
                                                                           cATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGaacgccaa                                                                           >  2:1919707/1‑150 (MQ=255)
                                                                            aTGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTAACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGaacgccaac                                                                          >  1:2578230/1‑150 (MQ=255)
                                                                             tGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGaacgccaacg                                                                         >  2:3468285/1‑150 (MQ=255)
                                                                              ggCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGaacgccaacgc                                                                        >  1:844064/1‑150 (MQ=255)
                                                                                 gACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCg                                                                     >  2:2677058/1‑150 (MQ=255)
                                                                                 gACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCg                                                                     >  2:3309800/1‑150 (MQ=255)
                                                                                    cATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTAACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTgc                                                                  <  2:2939801/150‑1 (MQ=255)
                                                                                        aaGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTg                                                              <  1:3302251/150‑1 (MQ=255)
                                                                                             tGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCgca                                                         <  2:389179/150‑1 (MQ=255)
                                                                                                 gCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGcagcag                                                     >  2:3408673/1‑150 (MQ=255)
                                                                                                    cGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGAcc                                                  >  1:4625041/1‑150 (MQ=255)
                                                                                                     gACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCa                                                 >  1:1420536/1‑150 (MQ=255)
                                                                                                        tGGCAGGCGTTTCAGTCGCCACCGTGTCCCGCGTCATTAATAATTCACCCAACGCCAGCGAAGCTTCCCGGCGGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTg                                              <  2:2237761/150‑1 (MQ=255)
                                                                                                         ggCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGa                                             <  1:3779490/150‑1 (MQ=255)
                                                                                                          gCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGaa                                            <  2:3637352/150‑1 (MQ=255)
                                                                                                            aGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAAc                                          >  2:440925/1‑150 (MQ=255)
                                                                                                                gTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAAcggtc                                      >  1:1659047/1‑150 (MQ=255)
                                                                                                                gTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAAcggtc                                      >  2:4598676/1‑150 (MQ=255)
                                                                                                                gTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTAACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAAcggtc                                      <  2:1093978/150‑1 (MQ=255)
                                                                                                                   tCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAAcggtcggt                                   <  1:3607680/150‑1 (MQ=255)
                                                                                                                   tCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAAcggtcggt                                   <  2:4168581/150‑1 (MQ=255)
                                                                                                                     aGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCt                                 <  1:4045678/150‑1 (MQ=255)
                                                                                                                        cGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGt                              <  1:1781419/150‑1 (MQ=255)
                                                                                                                        cGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGt                              <  1:4564967/150‑1 (MQ=255)
                                                                                                                          ccACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCg                            >  2:914771/1‑150 (MQ=255)
                                                                                                                          ccACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCATGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCg                            <  1:1458449/150‑1 (MQ=255)
                                                                                                                            aCCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGtt                          <  1:4268695/150‑1 (MQ=255)
                                                                                                                            aCCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGtt                          >  1:2828554/1‑150 (MQ=255)
                                                                                                                            aCCGTTTCCCGCGTCATTAATAATTAACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGtt                          >  2:3092536/1‑150 (MQ=255)
                                                                                                                               gTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGt                       <  2:2869960/150‑1 (MQ=255)
                                                                                                                                    ccGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGt                  <  1:2733834/150‑1 (MQ=255)
                                                                                                                                        gTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTcc              <  1:1945337/150‑1 (MQ=255)
                                                                                                                                        gTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTcc              <  1:2401744/150‑1 (MQ=255)
                                                                                                                                        gTCATTAATAATTCACCCAAAGCAAGCGAAGCTTCCCGGCTGGCTGTGAATAGTGAAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTcc              <  1:1213941/150‑1 (MQ=255)
                                                                                                                                          cATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGa            >  1:3477895/1‑150 (MQ=255)
                                                                                                                                           tttAATAATTCACCGAAAGCCAGCGACGCTTCCCGGCTGGATGTGCATAGTGAATTGTAGTCTTTTAGCTATCACCCGAAGGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAATAGTTCGGTCTGGTCGTTGGTGTTGTTTCCGAt           <  1:3613761/149‑1 (MQ=255)
                                                                                                                                           aTTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGAt           <  2:3312016/150‑1 (MQ=255)
                                                                                                                                               ataatTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGt       >  2:3041696/1‑150 (MQ=255)
                                                                                                                                                 aatTAACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGttt     >  1:4575045/1‑150 (MQ=255)
                                                                                                                                                    tCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTc  >  2:1345472/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CACGCCACCCCTTGAACCAACGGGCGTTTTCCGTAACACTGAAAGAATGTAAGCGTTTACCCACTAAGGTATTTTCATGGCGACCATAAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTC  >  CP009273/2969882‑2970179

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: