Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 271,650 G→T 10.0% P205P (CCG→CCT)  mmuP → CP4‑6 prophage; putative S‑methylmethionine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273271,6500GT10.0% 683.5 / 35.4 200P205P (CCG→CCT) mmuPCP4‑6 prophage; putative S‑methylmethionine transporter
Reads supporting (aligned to +/- strand):  ref base G (99/81);  new base T (6/14);  total (105/95)
Fisher's exact test for biased strand distribution p-value = 5.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 201 total aligned reads displayed.

TCGCTGGTCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTATCGGC  >  CP009273/271504‑271797
                                                                                                                                                  |                                                                                                                                                   
 cGCTGGTCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGAt                                                                                                                                                   <  2:1061459/148‑1 (MQ=255)
 cGCTGGTCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGAt                                                                                                                                                   >  1:2358249/1‑148 (MQ=255)
 cGCTGGTCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGAt                                                                                                                                                   >  1:2928338/1‑148 (MQ=255)
      gTCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGAGTTtga                                                                                                                                              >  1:847873/1‑148 (MQ=255)
        cAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTtgatg                                                                                                                                            <  2:3218608/148‑1 (MQ=255)
         aaaGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTtgatga                                                                                                                                           <  2:684287/148‑1 (MQ=255)
          aaGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGAc                                                                                                                                          <  1:3484993/148‑1 (MQ=255)
          aaGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGAc                                                                                                                                          >  1:2719944/1‑148 (MQ=255)
           aGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGCTTTTGATGACt                                                                                                                                         >  1:334569/1‑148 (MQ=255)
           aGTGGTCACTATCATCGCCTTTATCATAATAGGTTGTGCGGAGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGACCTCGCCGGGGCTGAGTAATAACACGGCAGAATGCTGGTTCACGCACGGTGGCTTAACGATTTtgatgaat                                                                                                                                         >  2:1362901/1‑146 (MQ=255)
             tGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTAt                                                                                                                                       <  2:2161432/148‑1 (MQ=255)
                 cACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTAtggtg                                                                                                                                   >  2:712860/1‑148 (MQ=255)
                 cACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTAtggtg                                                                                                                                   >  2:3291692/1‑148 (MQ=255)
                  aCTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTAtggtgg                                                                                                                                  <  2:2074317/148‑1 (MQ=255)
                    tATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCa                                                                                                                                <  1:3254378/148‑1 (MQ=255)
                     atcatcGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAg                                                                                                                               <  2:369289/148‑1 (MQ=255)
                     atcatcGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAg                                                                                                                               >  2:3278751/1‑148 (MQ=255)
                      tcatcGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTTATGACTATGGTGGCAGt                                                                                                                              >  2:272957/1‑148 (MQ=255)
                         tcGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGTTGACTATGGTGGCAGTGaa                                                                                                                           >  1:1030439/1‑148 (MQ=255)
                         tcGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGaa                                                                                                                           >  2:826362/1‑148 (MQ=255)
                          cGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAAc                                                                                                                          <  1:530783/148‑1 (MQ=255)
                             cTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACttt                                                                                                                       >  2:2227101/1‑148 (MQ=255)
                               ttATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGc                                                                                                                     >  1:1294620/1‑148 (MQ=255)
                               ttATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGc                                                                                                                     <  1:745272/148‑1 (MQ=255)
                                 atcatcCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCtt                                                                                                                   <  1:3326075/148‑1 (MQ=255)
                                  tcatcCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCttt                                                                                                                  >  2:253270/1‑148 (MQ=255)
                                  tcatcCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCttt                                                                                                                  >  2:1218083/1‑148 (MQ=255)
                                   catcCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCtttt                                                                                                                 <  2:2793680/148‑1 (MQ=255)
                                     tcCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTc                                                                                                               >  2:2664345/1‑148 (MQ=255)
                                      ccTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCg                                                                                                              >  1:2265915/1‑148 (MQ=255)
                                            tGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTAcc                                                                                                        >  2:1410223/1‑148 (MQ=255)
                                            tGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTAcc                                                                                                        >  2:1410078/1‑148 (MQ=255)
                                             ggggCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCg                                                                                                       <  1:2774124/148‑1 (MQ=255)
                                               ggcggcGATTTTAGGCTTTATTACGATGCAGGATGGCTCGCCCGCGCCCGGGCTGAGTAAAATAACCGCAGAAGGCTGGTTACCGCACCGTGGCTTAACGAATTTGATTACTATTGTTGCAGTGAACTTTTATTTTTAGGGTAACGAg                                                                                                     >  1:1637425/1‑148 (MQ=255)
                                                  ggcgATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCtt                                                                                                  <  1:569991/148‑1 (MQ=255)
                                                  ggcgATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCtt                                                                                                  <  1:3374902/148‑1 (MQ=255)
                                                       ttttCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCgg                                                                                             >  2:350121/1‑148 (MQ=255)
                                                       ttttCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCgg                                                                                             >  2:1533457/1‑148 (MQ=255)
                                                        tttCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGc                                                                                            >  2:2472251/1‑148 (MQ=255)
                                                            ggCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTg                                                                                        <  2:1925037/148‑1 (MQ=255)
                                                            ggCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTg                                                                                        <  2:3284248/148‑1 (MQ=255)
                                                                ttATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTgccgc                                                                                    <  2:1292061/148‑1 (MQ=255)
                                                                ttATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTgccgc                                                                                    >  2:2698010/1‑148 (MQ=255)
                                                                 tATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTgccgcc                                                                                   >  2:1090989/1‑148 (MQ=255)
                                                                  aTTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGccgccg                                                                                  <  1:1876043/148‑1 (MQ=255)
                                                                    tCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGt                                                                                <  1:2419892/148‑1 (MQ=255)
                                                                    tCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGt                                                                                >  2:1813104/1‑148 (MQ=255)
                                                                       gATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGaa                                                                             >  1:963256/1‑148 (MQ=255)
                                                                           cAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACgg                                                                         >  1:2432216/1‑148 (MQ=255)
                                                                            aGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGa                                                                        >  1:3045733/1‑148 (MQ=255)
                                                                              gATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGaaa                                                                      >  2:54975/1‑148 (MQ=255)
                                                                               aTGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGaaaa                                                                     <  1:770854/148‑1 (MQ=255)
                                                                               aTGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGaaaa                                                                     <  1:597308/148‑1 (MQ=255)
                                                                                   cTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGTCACGGAAAACCcg                                                                 >  2:1203784/1‑148 (MQ=255)
                                                                                   cTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCcg                                                                 <  1:956615/148‑1 (MQ=255)
                                                                                   cTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCGGACCTATTGTGATGACTATGGTGGCAGTGAACTTTGCTTTGTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCcg                                                                 >  1:2602172/1‑148 (MQ=255)
                                                                                      gCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCgcgc                                                              >  1:2472935/1‑148 (MQ=255)
                                                                                        ccGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCaa                                                            <  1:3162761/148‑1 (MQ=255)
                                                                                         cgcgCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCaaa                                                           <  2:1458316/148‑1 (MQ=255)
                                                                                         cgcgCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCaaa                                                           >  1:2857756/1‑148 (MQ=255)
                                                                                          gcgcCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAg                                                          <  2:1803055/148‑1 (MQ=255)
                                                                                           cgcCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGt                                                         <  2:761250/148‑1 (MQ=255)
                                                                                           cgcCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGt                                                         <  2:3331362/148‑1 (MQ=255)
                                                                                              cGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTAt                                                      <  2:2938212/148‑1 (MQ=255)
                                                                                                 ggCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCTATTTTGATGTCTATGGTGGCAGTGAACTTTGCTTGTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCACCGTTATccc                                                   >  2:2490185/1‑148 (MQ=255)
                                                                                                   cTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCgg                                                 <  2:2079685/148‑1 (MQ=255)
                                                                                                   cTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCgg                                                 >  2:886782/1‑148 (MQ=255)
                                                                                                    tGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGt                                                >  2:1726264/1‑148 (MQ=255)
                                                                                                         aaTATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGa                                           <  1:1370284/148‑1 (MQ=255)
                                                                                                          atatCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGAt                                          <  2:2895921/148‑1 (MQ=255)
                                                                                                          atatCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGAt                                          <  2:2386315/148‑1 (MQ=255)
                                                                                                            atCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTc                                        >  1:1491536/1‑148 (MQ=255)
                                                                                                              cACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGTGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGt                                      >  1:2979668/1‑148 (MQ=255)
                                                                                                               aCGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGta                                     >  1:2730033/1‑148 (MQ=255)
                                                                                                               aCGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGta                                     >  1:1384867/1‑148 (MQ=255)
                                                                                                                  gCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGtacta                                  >  1:2347509/1‑148 (MQ=255)
                                                                                                                    aGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTAcc                                <  2:1703102/148‑1 (MQ=255)
                                                                                                                      aaGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCAt                              <  2:3146924/148‑1 (MQ=255)
                                                                                                                      aaGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCAt                              >  2:2583009/1‑148 (MQ=255)
                                                                                                                      aaGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCAt                              >  1:1075865/1‑148 (MQ=255)
                                                                                                                       aGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATc                             >  2:1679767/1‑148 (MQ=255)
                                                                                                                        ggCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATcg                            <  1:1647262/148‑1 (MQ=255)
                                                                                                                           tGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATcgcgc                         >  1:2172503/1‑148 (MQ=255)
                                                                                                                            ggTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATcgcgcg                        >  1:237069/1‑148 (MQ=255)
                                                                                                                             gTTCCCGCACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGa                       >  2:394656/1‑148 (MQ=255)
                                                                                                                             gTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGa                       >  2:1980326/1‑148 (MQ=255)
                                                                                                                               tCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACt                     >  1:1721210/1‑148 (MQ=255)
                                                                                                                                cccGCACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTg                    <  2:660211/148‑1 (MQ=255)
                                                                                                                                    cACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGatta                <  1:332441/148‑1 (MQ=255)
                                                                                                                                    cACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGatta                <  1:1541903/148‑1 (MQ=255)
                                                                                                                                     aCGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGattat               <  2:2240872/148‑1 (MQ=255)
                                                                                                                                      cGGTGGCTTACCGATTTTGATGTCTATGGTGGGAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGattatt              <  2:1633829/148‑1 (MQ=255)
                                                                                                                                        gTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTAtttt            <  1:2572476/148‑1 (MQ=255)
                                                                                                                                          ggCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCt          <  2:2038052/148‑1 (MQ=255)
                                                                                                                                            cTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCttt        >  1:742393/1‑148 (MQ=255)
                                                                                                                                             ttACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTa       <  2:161246/148‑1 (MQ=255)
                                                                                                                                              tACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTAt      <  1:1045765/148‑1 (MQ=255)
                                                                                                                                               aCCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTATc     <  1:2023509/148‑1 (MQ=255)
                                                                                                                                               aCCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTATc     <  1:432561/148‑1 (MQ=255)
                                                                                                                                                ccGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTATCg    <  1:3413217/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCGCTGGTCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTATCGGC  >  CP009273/271504‑271797

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: ‑