Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,572,520 A→T 95.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,5200AT95.0% 499.5 / 12.6 159Y116N (TAC→AAC) hsdMDNA methyltransferase M
Reads supporting (aligned to +/- strand):  ref base A (4/4);  new base T (81/70);  total (85/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 160 total aligned reads displayed.

AGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGAACT  >  CP009273/4572378‑4572660
                                                                                                                                              |                                                                                                                                            
 gcagATGAATAATGGTTTTAATCAGCGGACTCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTc                                                                                                                                        >  1:2804322/1‑148 (MQ=255)
    gATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTtccagtccag                                                                                                                                     >  2:2024147/1‑148 (MQ=255)
        aataatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGaa                                                                                                                                 >  2:2562195/1‑148 (MQ=255)
        aataatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGaa                                                                                                                                 >  2:180807/1‑148 (MQ=255)
        aataatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGaa                                                                                                                                 >  1:2906134/1‑148 (MQ=255)
         ataatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAAt                                                                                                                                <  2:2814503/148‑1 (MQ=255)
          taatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATc                                                                                                                               >  2:2730810/1‑148 (MQ=255)
          taatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATc                                                                                                                               >  2:2239692/1‑148 (MQ=255)
          taatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATc                                                                                                                               <  1:2680173/148‑1 (MQ=255)
          taatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATc                                                                                                                               >  2:2733245/1‑148 (MQ=255)
              ggTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCata                                                                                                                           <  2:589419/148‑1 (MQ=255)
                ttttAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATAtt                                                                                                                         >  1:2517876/1‑148 (MQ=255)
                      tCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGAc                                                                                                                   >  1:2083386/1‑148 (MQ=255)
                      tCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGAc                                                                                                                   <  1:1902872/148‑1 (MQ=255)
                      tCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGAc                                                                                                                   <  2:235645/148‑1 (MQ=255)
                      tCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGAc                                                                                                                   <  1:2639435/148‑1 (MQ=255)
                       cAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGAcc                                                                                                                  >  1:972326/1‑148 (MQ=255)
                          cGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGt                                                                                                               >  1:1013790/1‑148 (MQ=255)
                                  gggTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTAt                                                                                                       >  1:361033/1‑148 (MQ=255)
                                  gggTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTAt                                                                                                       <  1:1814635/148‑1 (MQ=255)
                                   ggTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTAtt                                                                                                      >  2:2413291/1‑148 (MQ=255)
                                      gAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGt                                                                                                   <  2:3227589/148‑1 (MQ=255)
                                        aGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGttt                                                                                                 <  2:2874223/148‑1 (MQ=255)
                                          tACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCg                                                                                               >  1:1659656/1‑148 (MQ=255)
                                           aCTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCgg                                                                                              <  1:1010409/148‑1 (MQ=255)
                                             tGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCt                                                                                            <  1:686591/148‑1 (MQ=255)
                                             tGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCt                                                                                            >  2:1588507/1‑148 (MQ=255)
                                                ccTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCgg                                                                                         >  1:3134585/1‑148 (MQ=255)
                                                 cTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGt                                                                                        <  1:1283302/148‑1 (MQ=255)
                                                 cTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGt                                                                                        <  1:1283250/148‑1 (MQ=255)
                                                    cACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGag                                                                                     >  2:2304468/1‑147 (MQ=255)
                                                      ccAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAtgg                                                                                   <  2:1659656/148‑1 (MQ=255)
                                                      ccAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAtgg                                                                                   <  1:2413078/148‑1 (MQ=255)
                                                      ccAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTGACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAtgg                                                                                   <  2:1121685/148‑1 (MQ=255)
                                                       cAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACCTCTCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAtggt                                                                                  <  2:3339349/148‑1 (MQ=255)
                                                          aCTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATggtggt                                                                               >  2:115883/1‑148 (MQ=255)
                                                          aCTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATggtggt                                                                               >  2:114882/1‑148 (MQ=255)
                                                           cTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTa                                                                              <  2:3022892/148‑1 (MQ=255)
                                                             tGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACt                                                                            <  1:180893/148‑1 (MQ=255)
                                                               gTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTaa                                                                          >  2:1810058/1‑148 (MQ=255)
                                                                tttCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGATTCGGCTCGGTGATGGTGGTACTAAc                                                                         <  2:2605697/148‑1 (MQ=255)
                                                                 ttCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACa                                                                        <  2:902090/148‑1 (MQ=255)
                                                                  tCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACAt                                                                       <  2:2538976/148‑1 (MQ=255)
                                                                   cATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACAtt                                                                      <  2:2100234/148‑1 (MQ=255)
                                                                   cATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACAtt                                                                      >  1:997304/1‑148 (MQ=255)
                                                                    aTTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTa                                                                     >  1:152606/1‑148 (MQ=255)
                                                                         cgTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGaaa                                                                >  1:3222228/1‑148 (MQ=255)
                                                                         cgTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGaaa                                                                >  1:3222190/1‑148 (MQ=255)
                                                                           ttcttcTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGaaaaa                                                              <  1:753287/148‑1 (MQ=255)
                                                                            tcttctGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAAc                                                             >  1:2781492/1‑148 (MQ=255)
                                                                               tctGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGc                                                          <  2:2382718/148‑1 (MQ=255)
                                                                                ctGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGcc                                                         <  2:2435382/148‑1 (MQ=255)
                                                                                   cAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGt                                                      >  2:2490127/1‑148 (MQ=255)
                                                                                        gCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAg                                                 <  2:426257/148‑1 (MQ=255)
                                                                                         cccTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGc                                                >  2:730768/1‑148 (MQ=255)
                                                                                         cccTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGc                                                <  1:3102445/148‑1 (MQ=255)
                                                                                          ccTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCt                                               >  1:1070818/1‑148 (MQ=255)
                                                                                           cTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCtt                                              <  1:3400430/148‑1 (MQ=255)
                                                                                            ttCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGGTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCttt                                             >  1:640044/1‑148 (MQ=255)
                                                                                            ttCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCttt                                             >  1:2182792/1‑148 (MQ=255)
                                                                                             tCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCtttt                                            <  1:745614/148‑1 (MQ=255)
                                                                                              cGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCttttt                                           <  2:3403110/148‑1 (MQ=255)
                                                                                                tACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGt                                         >  1:523560/1‑148 (MQ=255)
                                                                                                tACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGt                                         >  1:1624114/1‑148 (MQ=255)
                                                                                                 aCATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTc                                        >  1:3479557/1‑148 (MQ=255)
                                                                                                 aCATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTc                                        >  1:250199/1‑148 (MQ=255)
                                                                                                 aCATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTc                                        >  1:1716098/1‑148 (MQ=255)
                                                                                                  cATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCa                                       >  1:1143199/1‑148 (MQ=255)
                                                                                                    tatCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATc                                     <  2:275450/148‑1 (MQ=255)
                                                                                                    tatCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTTCCTGTACCAGCTTTTTTTCATc                                     >  1:117331/1‑148 (MQ=255)
                                                                                                    tatCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGGTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATc                                     <  2:3067781/148‑1 (MQ=255)
                                                                                                     atCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCTGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCt                                    <  1:2367722/148‑1 (MQ=255)
                                                                                                     atCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCt                                    <  1:735139/148‑1 (MQ=255)
                                                                                                     atCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCt                                    >  1:1951160/1‑148 (MQ=255)
                                                                                                     atCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCt                                    >  1:2942665/1‑148 (MQ=255)
                                                                                                       cGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTc                                  >  1:3524543/1‑148 (MQ=255)
                                                                                                            aaGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTa                             >  1:2932158/1‑148 (MQ=255)
                                                                                                             aGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTaa                            >  2:737176/1‑148 (MQ=255)
                                                                                                              gTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTaaa                           >  2:2241628/1‑148 (MQ=255)
                                                                                                                cATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATg                         >  1:1358118/1‑148 (MQ=255)
                                                                                                                   cgcgcgACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCAc                      >  2:1468651/1‑148 (MQ=255)
                                                                                                                    gcgcgACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACg                     <  1:299971/148‑1 (MQ=255)
                                                                                                                     cgcgACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGa                    >  1:1448452/1‑148 (MQ=255)
                                                                                                                        gACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCa                 <  1:2550346/148‑1 (MQ=255)
                                                                                                                              ccGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTc           >  1:3514611/1‑148 (MQ=255)
                                                                                                                               cGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCg          >  2:1656674/1‑148 (MQ=255)
                                                                                                                                gTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCgg         <  2:387134/148‑1 (MQ=255)
                                                                                                                                gTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCgg         <  1:2733245/148‑1 (MQ=255)
                                                                                                                                gTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCgg         <  1:2733136/148‑1 (MQ=255)
                                                                                                                                 tGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGt        >  2:2589365/1‑148 (MQ=255)
                                                                                                                                 tGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGt        >  2:2588234/1‑148 (MQ=255)
                                                                                                                                  gcgcgcCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTa       <  2:3107910/148‑1 (MQ=255)
                                                                                                                                   cgcgcCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAg      <  2:1013790/148‑1 (MQ=255)
                                                                                                                                    gcgcCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGa     <  1:811838/148‑1 (MQ=255)
                                                                                                                                    gcgcCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGa     <  2:2691213/148‑1 (MQ=255)
                                                                                                                                     cgcCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaa    >  2:2191238/1‑148 (MQ=255)
                                                                                                                                     cgcCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaa    >  2:2981322/1‑148 (MQ=255)
                                                                                                                                      gcCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaac   <  1:2024147/148‑1 (MQ=255)
                                                                                                                                       ccGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaact  >  2:1380987/1‑148 (MQ=255)
                                                                                                                                       ccGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCATTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaact  >  2:3283731/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                            
AGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGAACT  >  CP009273/4572378‑4572660

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: