Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,446,543 G→A 9.8% pseudogene (1938/2646 nt) yfcU ← pseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,446,5430GA 9.8% 947.8 / 66.6 287pseudogene (1938/2646 nt)yfcUpseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (118/140);  new base A (10/18);  total (128/159)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 296 total aligned reads displayed.

CGGCAACGCCATCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGT  >  CP009273/2446395‑2446691
                                                                                                                                                    |                                                                                                                                                    
cGGCAACGCCATCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAgg                                                                                                                                                      >  2:2520955/1‑149 (MQ=255)
cGGCAACGCCATCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAgg                                                                                                                                                      >  1:2537748/1‑149 (MQ=255)
 ggCAACGCCATCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGt                                                                                                                                                     <  2:2490858/149‑1 (MQ=255)
  gCAACGCCATCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTc                                                                                                                                                    >  1:3125901/1‑149 (MQ=255)
  gCAACGCCATCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTc                                                                                                                                                    >  2:1714695/1‑149 (MQ=255)
     aCGCCATCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAAc                                                                                                                                                 >  2:2269958/1‑149 (MQ=255)
      cGCCATCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAAcc                                                                                                                                                <  2:494110/149‑1 (MQ=255)
      cGCCATCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAAcc                                                                                                                                                <  2:4664036/149‑1 (MQ=255)
        ccATCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTg                                                                                                                                              <  2:4031096/149‑1 (MQ=255)
          aTCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTgcg                                                                                                                                            >  2:843938/1‑149 (MQ=255)
           tCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTgcgc                                                                                                                                           >  1:2019512/1‑149 (MQ=255)
            cGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGcc                                                                                                                                          <  2:775139/149‑1 (MQ=255)
              gCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAg                                                                                                                                        >  2:928788/1‑149 (MQ=255)
              gCATCAATCAGCAGGCGTGTCCCGCCCATAATCTGGGTACGGTGAAGTGACCCGACGTGGGTCGTCAGCGAGGCGCCGCCCTGTAACGACAAGCCCGCCGAGGTGTATTGCCCTTCATCGTAGTTCGCACTGAGATCAACCTGAGCCAg                                                                                                                                        >  1:3215748/1‑149 (MQ=255)
                atcaatcaGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCg                                                                                                                                      >  1:3393259/1‑149 (MQ=255)
                  caatcaGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGaa                                                                                                                                    >  2:2980846/1‑149 (MQ=255)
                      cagcagGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACcat                                                                                                                                <  1:3597964/149‑1 (MQ=255)
                        gcagGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACcatca                                                                                                                              >  2:2949368/1‑149 (MQ=255)
                           ggCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATgg                                                                                                                           >  2:1250700/1‑149 (MQ=255)
                             cGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCt                                                                                                                         <  1:2687906/149‑1 (MQ=255)
                             cGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCt                                                                                                                         <  1:4352144/149‑1 (MQ=255)
                              gtgtCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTg                                                                                                                        >  2:3641665/1‑149 (MQ=255)
                                   ccGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTgtagta                                                                                                                   <  1:1548515/149‑1 (MQ=255)
                                   ccGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTgtagta                                                                                                                   <  2:3806349/149‑1 (MQ=255)
                                      cccATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGcc                                                                                                                >  2:1931597/1‑149 (MQ=255)
                                       ccATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCa                                                                                                               <  2:5189831/149‑1 (MQ=255)
                                         atatTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATc                                                                                                             >  1:697440/1‑149 (MQ=255)
                                          tatTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCa                                                                                                            <  2:3019047/149‑1 (MQ=255)
                                          tatTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGATCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCa                                                                                                            <  2:4790037/149‑1 (MQ=255)
                                           atTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCaa                                                                                                           <  1:3469610/149‑1 (MQ=255)
                                           atTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCaa                                                                                                           <  2:3848389/149‑1 (MQ=255)
                                              cTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGc                                                                                                        <  2:3278465/149‑1 (MQ=255)
                                               tGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCt                                                                                                       >  1:1275086/1‑149 (MQ=255)
                                               tGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGATCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCt                                                                                                       >  1:5033622/1‑149 (MQ=255)
                                                  gTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGgt                                                                                                    <  2:1750282/149‑1 (MQ=255)
                                                    aCGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGgtgt                                                                                                  <  1:1594030/149‑1 (MQ=255)
                                                        tGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGttt                                                                                              <  1:2287442/149‑1 (MQ=255)
                                                          aaGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGATCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGt                                                                                            <  2:1493237/149‑1 (MQ=255)
                                                              gcgcCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACt                                                                                        <  1:4115895/149‑1 (MQ=255)
                                                               cgccgcCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTg                                                                                       <  2:2815979/149‑1 (MQ=255)
                                                                gccgccGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTgg                                                                                      <  1:4165454/149‑1 (MQ=255)
                                                                 ccgccgTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGt                                                                                     >  2:1628358/1‑149 (MQ=255)
                                                                 ccgccgTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGATCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGt                                                                                     >  2:475313/1‑149 (MQ=255)
                                                                   gccgTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGc                                                                                   <  1:3186761/149‑1 (MQ=255)
                                                                       tGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGAt                                                                               <  1:784032/149‑1 (MQ=255)
                                                                         ggTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGt                                                                             <  2:4150002/149‑1 (MQ=255)
                                                                          gtagtaAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGACGATGtt                                                                            >  2:3103825/1‑149 (MQ=255)
                                                                           tagtaAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTc                                                                           <  1:2442220/149‑1 (MQ=255)
                                                                              taAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAAc                                                                        <  1:588205/149‑1 (MQ=255)
                                                                               aaGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACt                                                                       <  2:1554462/149‑1 (MQ=255)
                                                                                 gCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGa                                                                     >  1:1115373/1‑149 (MQ=255)
                                                                                    ttGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAg                                                                  >  2:1115786/1‑149 (MQ=255)
                                                                                    ttGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAg                                                                  <  2:2589814/149‑1 (MQ=255)
                                                                                    ttGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAg                                                                  >  2:1732116/1‑149 (MQ=255)
                                                                                      gcgcCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTg                                                                <  1:4735675/149‑1 (MQ=255)
                                                                                                       cAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGacac                                               >  2:3557170/1‑149 (MQ=255)
                                                                                                       cAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGacac                                               >  1:704859/1‑149 (MQ=255)
                                                                                                         aGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACgg                                             <  2:2115697/149‑1 (MQ=255)
                                                                                                          gCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGc                                            <  2:3125925/149‑1 (MQ=255)
                                                                                                           cccGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGATCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCt                                           <  1:2779456/149‑1 (MQ=255)
                                                                                                            ccGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTa                                          >  2:1016018/1‑149 (MQ=255)
                                                                                                               cAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTaaaa                                       <  1:1949259/149‑1 (MQ=255)
                                                                                                                 gagGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATa                                     <  1:1847135/149‑1 (MQ=255)
                                                                                                                    gtgtATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAAcc                                  >  1:1000496/1‑149 (MQ=255)
                                                                                                                     tgtATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCg                                 >  1:3574055/1‑149 (MQ=255)
                                                                                                                      gtATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGa                                >  2:553555/1‑149 (MQ=255)
                                                                                                                         ttGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCt                             <  1:4840987/149‑1 (MQ=255)
                                                                                                                         ttGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCt                             >  2:792567/1‑149 (MQ=255)
                                                                                                                            cccTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGgc                          >  1:3936051/1‑149 (MQ=255)
                                                                                                                            cccTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGgc                          >  2:269207/1‑149 (MQ=255)
                                                                                                                            cccTTCATGGTAGTTCGCACTGAGATCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGgc                          <  2:4950906/149‑1 (MQ=255)
                                                                                                                             ccTTCATGGTAGTTCGCACTGAGATCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGgct                         <  1:2091130/149‑1 (MQ=255)
                                                                                                                             ccTTCATGGTAGTTCGCACTGAGATCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGgct                         <  1:2090680/149‑1 (MQ=255)
                                                                                                                              cTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGgctg                        <  1:746312/149‑1 (MQ=255)
                                                                                                                                tCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGgctgct                      <  1:4017581/149‑1 (MQ=255)
                                                                                                                                  aTGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTgt                    <  1:1938060/149‑1 (MQ=255)
                                                                                                                                  aTGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTgt                    >  1:3646449/1‑149 (MQ=255)
                                                                                                                                  aTGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTgt                    >  1:3966695/1‑149 (MQ=255)
                                                                                                                                   tGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTgtc                   <  1:2844970/149‑1 (MQ=255)
                                                                                                                                   tGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTgtc                   <  1:555962/149‑1 (MQ=255)
                                                                                                                                   tGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTgtc                   >  1:834707/1‑149 (MQ=255)
                                                                                                                                    ggTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTgtcc                  <  1:786467/149‑1 (MQ=255)
                                                                                                                                       aGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTgtccgtc               >  2:4124435/1‑149 (MQ=255)
                                                                                                                                        gTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTgtccgtcc              <  2:1202601/149‑1 (MQ=255)
                                                                                                                                         ttCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTccc             <  1:3738098/149‑1 (MQ=255)
                                                                                                                                           cGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCAc           <  1:1913131/149‑1 (MQ=255)
                                                                                                                                           cGCACTGAGATCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCAc           <  1:4928973/149‑1 (MQ=255)
                                                                                                                                            gCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACt          >  2:1038122/1‑149 (MQ=255)
                                                                                                                                              aCTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGc        >  1:2870768/1‑149 (MQ=255)
                                                                                                                                               cTGAGGTCAACCTGCGCCAGCTAACCATCATGGCTGTAGTAGCCATAAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGcc       <  1:3499502/149‑1 (MQ=255)
                                                                                                                                               cTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGcc       >  1:4772803/1‑149 (MQ=255)
                                                                                                                                               cTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGcc       <  1:3143640/149‑1 (MQ=255)
                                                                                                                                                tGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCa      >  2:1316956/1‑149 (MQ=255)
                                                                                                                                                tGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCa      >  2:4355312/1‑149 (MQ=255)
                                                                                                                                                 gAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCAt     >  1:4696431/1‑149 (MQ=255)
                                                                                                                                                  aGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATa    >  1:2742497/1‑149 (MQ=255)
                                                                                                                                                   ggTCAACCTGCGCCAGCGAACCATCATGGCTGTATTAGCCATCAACGCTGGTGTGTTTTTCACTTTTGCCTATGTTCATCTGATTGTGTGTCGCTTCATCGACACGTCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAg   <  1:4880760/149‑1 (MQ=255)
                                                                                                                                                   ggTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAg   >  1:865503/1‑149 (MQ=255)
                                                                                                                                                    gTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGt  >  2:5139463/1‑149 (MQ=255)
                                                                                                                                                    gTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGt  <  1:4611185/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGGCAACGCCATCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGT  >  CP009273/2446395‑2446691

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: