Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,572,520 A→T 86.0% Y116N (TAC→AAC)  hsdM ← DNA methyltransferase M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,572,5200AT86.0% 692.1 / 87.8 241Y116N (TAC→AAC) hsdMDNA methyltransferase M
Reads supporting (aligned to +/- strand):  ref base A (11/23);  new base T (104/103);  total (115/126)
Fisher's exact test for biased strand distribution p-value = 6.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 243 total aligned reads displayed.

TTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGAACTGCAAC  >  CP009273/4572375‑4572665
                                                                                                                                                 |                                                                                                                                                 
 tCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAg                                                                                                                                               <  1:4295736/149‑1 (MQ=255)
 tCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAg                                                                                                                                               <  1:4560675/149‑1 (MQ=255)
   agcagATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTc                                                                                                                                             >  1:1875533/1‑149 (MQ=255)
   agcagATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTc                                                                                                                                             >  2:537199/1‑149 (MQ=255)
    gcagATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTcc                                                                                                                                            >  1:3911320/1‑149 (MQ=255)
    gcagATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTcc                                                                                                                                            >  2:4314678/1‑149 (MQ=255)
      agATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTtccagtccag                                                                                                                                          >  2:102148/1‑149 (MQ=255)
       gATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGc                                                                                                                                         <  2:683352/149‑1 (MQ=255)
       gATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGc                                                                                                                                         >  2:2307330/1‑149 (MQ=255)
        aTGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCg                                                                                                                                        >  1:4644834/1‑149 (MQ=255)
          gAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGaa                                                                                                                                      <  1:1431273/149‑1 (MQ=255)
             taatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATcc                                                                                                                                   >  1:63716/1‑149 (MQ=255)
             taatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATcc                                                                                                                                   >  1:3594323/1‑149 (MQ=255)
             taatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATcc                                                                                                                                   >  2:4073562/1‑149 (MQ=255)
               atGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCat                                                                                                                                 >  1:2782040/1‑149 (MQ=255)
                  gTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATAtt                                                                                                                              >  1:2065889/1‑149 (MQ=255)
                   ttttAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTg                                                                                                                             >  1:1701418/1‑149 (MQ=255)
                     ttAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCt                                                                                                                           <  2:3710772/149‑1 (MQ=255)
                      tAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTg                                                                                                                          <  1:3335697/149‑1 (MQ=255)
                      tAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTg                                                                                                                          <  1:3335737/149‑1 (MQ=255)
                       aaTCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGa                                                                                                                         <  2:4718668/149‑1 (MQ=255)
                          cAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCa                                                                                                                      <  1:4045784/149‑1 (MQ=255)
                              ggACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGc                                                                                                                  >  2:1704350/1‑149 (MQ=255)
                                aCGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCgg                                                                                                                <  2:703346/149‑1 (MQ=255)
                                  gcgGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGtt                                                                                                              <  1:533088/149‑1 (MQ=255)
                                  gcgGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGtt                                                                                                              <  1:4445427/149‑1 (MQ=255)
                                   cgGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTa                                                                                                             >  1:3772258/1‑149 (MQ=255)
                                    ggggTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTAt                                                                                                            >  2:2718029/1‑149 (MQ=255)
                                    ggggTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTAt                                                                                                            >  2:326940/1‑149 (MQ=255)
                                     gggTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTAtt                                                                                                           >  1:3869897/1‑149 (MQ=255)
                                      ggTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTAttt                                                                                                          >  2:1071599/1‑149 (MQ=255)
                                      ggTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTAttt                                                                                                          >  2:1852840/1‑149 (MQ=255)
                                             tACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCgg                                                                                                   >  1:1839454/1‑149 (MQ=255)
                                                tGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTc                                                                                                <  1:3038641/149‑1 (MQ=255)
                                                 ggCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCg                                                                                               >  2:1400239/1‑149 (MQ=255)
                                                  gCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCgg                                                                                              <  2:1387345/149‑1 (MQ=255)
                                                  gCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCgg                                                                                              <  1:2467027/149‑1 (MQ=255)
                                                  gCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATACATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCgg                                                                                              <  1:1939691/149‑1 (MQ=255)
                                                    cTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTg                                                                                            >  2:2091622/1‑149 (MQ=255)
                                                    cTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTg                                                                                            <  1:3681852/149‑1 (MQ=255)
                                                      gCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAt                                                                                          <  1:969327/149‑1 (MQ=255)
                                                         ccAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAtggt                                                                                       <  1:2600687/149‑1 (MQ=255)
                                                         ccAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAtggt                                                                                       >  2:619894/1‑149 (MQ=255)
                                                            gACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATggtggt                                                                                    <  1:4210776/149‑1 (MQ=255)
                                                             aCTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTa                                                                                   >  1:4752795/1‑149 (MQ=255)
                                                               ttGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACt                                                                                 >  1:4431308/1‑149 (MQ=255)
                                                               ttGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACt                                                                                 >  2:364128/1‑149 (MQ=255)
                                                                 ggTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTaa                                                                               >  1:1212051/1‑149 (MQ=255)
                                                                  gTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAAc                                                                              <  2:1340441/149‑1 (MQ=255)
                                                                  gTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAAc                                                                              <  1:2893207/149‑1 (MQ=255)
                                                                   tttCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACa                                                                             >  1:695459/1‑149 (MQ=255)
                                                                    ttCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACAt                                                                            >  2:4615721/1‑149 (MQ=255)
                                                                    ttCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACAt                                                                            <  1:2672259/149‑1 (MQ=255)
                                                                    ttCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACAt                                                                            <  1:1183741/149‑1 (MQ=255)
                                                                      cATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTa                                                                          <  2:2034255/149‑1 (MQ=255)
                                                                       aTTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTAt                                                                         >  1:3752998/1‑149 (MQ=255)
                                                                       aTTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTAt                                                                         <  1:2800762/149‑1 (MQ=255)
                                                                        ttCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATg                                                                        >  2:3140613/1‑149 (MQ=255)
                                                                        ttCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATg                                                                        >  2:2007276/1‑149 (MQ=255)
                                                                           gcgTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGaaa                                                                     <  1:2697985/149‑1 (MQ=255)
                                                                              ttcttcTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAAc                                                                  >  2:2409100/1‑149 (MQ=255)
                                                                              ttcttcTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAAc                                                                  <  2:3146412/149‑1 (MQ=255)
                                                                                    tGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGACTCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTg                                                            >  2:223817/1‑149 (MQ=255)
                                                                                    tGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACAGTATGAAAAACAGCCTg                                                            >  2:1208308/1‑149 (MQ=255)
                                                                                     gCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGt                                                           >  1:233194/1‑149 (MQ=255)
                                                                                       aaCAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTAc                                                         >  1:2872897/1‑149 (MQ=255)
                                                                                       aaCAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTAc                                                         <  1:3177580/149‑1 (MQ=255)
                                                                                       aaCAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCGCGGTGAGTGTGGTACTAACATTATTAAAAACTGCCTGTaa                                                         >  1:559083/1‑148 (MQ=255)
                                                                                        aCAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTAcc                                                        >  2:3982894/1‑149 (MQ=255)
                                                                                         cAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCa                                                       <  2:4260455/149‑1 (MQ=255)
                                                                                           gCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGc                                                     <  1:2743845/149‑1 (MQ=255)
                                                                                            cccTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCt                                                    <  1:1902268/149‑1 (MQ=255)
                                                                                             ccTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGGACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGGACCAGCtt                                                   <  2:4111512/149‑1 (MQ=255)
                                                                                               ttCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCtttt                                                 <  2:4747037/149‑1 (MQ=255)
                                                                                                tCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCttttt                                                <  2:2556490/149‑1 (MQ=255)
                                                                                                tCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCttttt                                                <  2:4376801/149‑1 (MQ=255)
                                                                                                  gTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGt                                              >  2:535093/1‑149 (MQ=255)
                                                                                                  gTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGt                                              <  2:3493476/149‑1 (MQ=255)
                                                                                                   tACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTc                                             >  1:1570435/1‑149 (MQ=255)
                                                                                                      atatCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATc                                          >  2:1589213/1‑149 (MQ=255)
                                                                                                      atatCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATc                                          >  2:3922027/1‑149 (MQ=255)
                                                                                                      atatCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATc                                          >  2:1532282/1‑149 (MQ=255)
                                                                                                               aaGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTaa                                 >  1:1553819/1‑149 (MQ=255)
                                                                                                                aGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTaaa                                <  1:1727254/149‑1 (MQ=255)
                                                                                                                  tCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATg                              <  1:4390965/149‑1 (MQ=255)
                                                                                                                    aTCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAACTGCa                            <  1:814170/149‑1 (MQ=255)
                                                                                                                    aTCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCa                            <  1:1049098/149‑1 (MQ=255)
                                                                                                                     tCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCAc                           >  2:2862412/1‑149 (MQ=255)
                                                                                                                     tCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCAGATTGCTGACCAGTGCGGATAGTTGTTTCGGCTCTGTGAAGGTGGTACTACAATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCAc                           >  2:348362/1‑149 (MQ=255)
                                                                                                                     tCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGAAc                           >  2:3257553/1‑149 (MQ=255)
                                                                                                                          cgACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCa                      <  2:3594323/149‑1 (MQ=255)
                                                                                                                          cgACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCa                      <  2:2877539/149‑1 (MQ=255)
                                                                                                                           gACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCAt                     <  2:934884/149‑1 (MQ=255)
                                                                                                                             cTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCAttt                   <  2:2466324/149‑1 (MQ=255)
                                                                                                                              ttACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCAtttt                  >  2:2988504/1‑149 (MQ=255)
                                                                                                                                   gTGCGCGCCGTTGGTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGt             <  2:3380100/149‑1 (MQ=255)
                                                                                                                                      cgcgcCGTTGGTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGa          <  2:3782973/149‑1 (MQ=255)
                                                                                                                                       gcgcCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaa         >  1:3485087/1‑149 (MQ=255)
                                                                                                                                        cgcCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaac        <  1:3830378/149‑1 (MQ=255)
                                                                                                                                              tGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaactgcaac  >  1:1402334/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGAACTGCAAC  >  CP009273/4572375‑4572665

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: