Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,749,958 C→A 8.1% I147I (ATC→ATA intA → CP4‑57 prophage; integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,749,9580CA 8.1% 764.7 / 26.6 210I147I (ATC→ATAintACP4‑57 prophage; integrase
Reads supporting (aligned to +/- strand):  ref base C (98/95);  new base A (9/8);  total (107/103)
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.97e-01

Only 100 of 212 total aligned reads displayed.

TTCTCAAGAGGCCAAAACCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTA  >  CP009273/2749811‑2750103
                                                                                                                                                   |                                                                                                                                                 
 tctcAAGAGGCCAAAACCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCa                                                                                                                                                  <  1:3021135/148‑1 (MQ=255)
 tctcAAGAGGCCAAAACCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCa                                                                                                                                                  <  2:2516543/148‑1 (MQ=255)
 tctcAAGAGGCCAAAACCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCa                                                                                                                                                  <  1:2474550/148‑1 (MQ=255)
  ctcAAGAGGCCAAAACCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCAc                                                                                                                                                 <  1:551510/148‑1 (MQ=255)
    cAAGAGGCCAAAACCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTg                                                                                                                                               <  2:19473/148‑1 (MQ=255)
    cAAGAGGCCAAAACCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTg                                                                                                                                               <  1:1556629/148‑1 (MQ=255)
    cAAGAGGCCAAAACCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTg                                                                                                                                               <  2:369531/148‑1 (MQ=255)
     aaGAGGCAAAAACCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTga                                                                                                                                              >  2:3132359/1‑148 (MQ=255)
        agGCCAAAACCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGAt                                                                                                                                           >  2:127981/1‑148 (MQ=255)
        agGCCAAAACCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGAt                                                                                                                                           >  1:2706685/1‑148 (MQ=255)
        agGCCAAAACCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATAACTGAGAt                                                                                                                                           >  1:2322494/1‑148 (MQ=255)
        agGCCAAAACCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATAACTGAGAt                                                                                                                                           >  1:662006/1‑148 (MQ=255)
             aaaaCCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAgg                                                                                                                                      >  1:124460/1‑148 (MQ=255)
             aaaaCCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAgg                                                                                                                                      >  2:731238/1‑148 (MQ=255)
                aCCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTc                                                                                                                                   >  2:1738056/1‑148 (MQ=255)
                aCCAATAACTTCTTGTTAGTTGCCGAGCGTTGTTGGAATGTGAAGAAAACCATCGTAACACAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATATTCCCGGCAATCTGTGATATCAGTATCACTGAGATTAAGGCTc                                                                                                                                   >  2:1261802/1‑148 (MQ=255)
                 ccAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCa                                                                                                                                  >  1:2561104/1‑148 (MQ=255)
                   aaTACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATa                                                                                                                                <  1:2368536/148‑1 (MQ=255)
                    aTACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATAc                                                                                                                               >  2:1968223/1‑148 (MQ=255)
                       ccTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATActct                                                                                                                            >  1:470284/1‑148 (MQ=255)
                       ccTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATActct                                                                                                                            >  1:1708035/1‑148 (MQ=255)
                        cttcttGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTg                                                                                                                           >  1:2455936/1‑148 (MQ=255)
                         ttcttGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATAACTGAGATTAAGGCTCATACTCggg                                                                                                                          <  1:864500/148‑1 (MQ=255)
                         ttcttGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATAACTGAGATTAAGGCTCATACTCTgg                                                                                                                          >  2:232465/1‑148 (MQ=255)
                          tcttGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGt                                                                                                                         <  2:1462950/148‑1 (MQ=255)
                             tGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTaa                                                                                                                      <  2:2823190/148‑1 (MQ=255)
                              gTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTaaa                                                                                                                     <  2:1721895/148‑1 (MQ=255)
                                tAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGc                                                                                                                   >  1:1940273/1‑148 (MQ=255)
                                   ttGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATCCTCTGGTTAAAGCAGt                                                                                                                >  1:581089/1‑148 (MQ=255)
                                        gAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGc                                                                                                           >  1:2711563/1‑148 (MQ=255)
                                          gCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCg                                                                                                         >  1:3226417/1‑148 (MQ=255)
                                          gCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCg                                                                                                         >  2:317740/1‑148 (MQ=255)
                                           cGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCgg                                                                                                        >  2:1920416/1‑148 (MQ=255)
                                           cGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCgg                                                                                                        >  2:2777914/1‑148 (MQ=255)
                                              tggtggAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTc                                                                                                     >  1:2296111/1‑148 (MQ=255)
                                                     aTGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCag                                                                                              <  1:1452222/148‑1 (MQ=255)
                                                     aTGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCag                                                                                              <  1:1053084/148‑1 (MQ=255)
                                                         gaagaaAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGt                                                                                          >  2:2053781/1‑148 (MQ=255)
                                                          aagaaAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTg                                                                                         >  2:3121186/1‑148 (MQ=255)
                                                           agaaAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGc                                                                                        >  2:601566/1‑148 (MQ=255)
                                                             aaaaCCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCAt                                                                                      >  1:1146136/1‑148 (MQ=255)
                                                             aaaaCCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCAt                                                                                      <  1:2217786/148‑1 (MQ=255)
                                                             aaaaCCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCAt                                                                                      >  1:1644090/1‑148 (MQ=255)
                                                                 ccAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTaga                                                                                  >  2:1818083/1‑148 (MQ=255)
                                                                  cAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTagag                                                                                 >  1:748558/1‑148 (MQ=255)
                                                                      gTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTg                                                                             <  1:1151239/148‑1 (MQ=255)
                                                                       tAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGt                                                                            <  2:324273/148‑1 (MQ=255)
                                                                        aaCAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGtt                                                                           >  1:459663/1‑148 (MQ=255)
                                                                               gACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCgccgcc                                                                    >  1:1828987/1‑148 (MQ=255)
                                                                                aCTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCt                                                                   >  2:761671/1‑148 (MQ=255)
                                                                                  tATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCttt                                                                 <  1:1702802/148‑1 (MQ=255)
                                                                                   aTGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTg                                                                <  1:1968223/148‑1 (MQ=255)
                                                                                   aTGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTg                                                                <  2:843143/148‑1 (MQ=255)
                                                                                    tGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGt                                                               >  2:3205933/1‑148 (MQ=255)
                                                                                       cgacgaTATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAg                                                            <  2:513220/148‑1 (MQ=255)
                                                                                       cgacgaTATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAg                                                            <  1:730276/148‑1 (MQ=255)
                                                                                       cgacgaTATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAg                                                            <  2:124460/148‑1 (MQ=255)
                                                                                         acgaTATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCg                                                          >  2:2999058/1‑148 (MQ=255)
                                                                                             tatCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTAtt                                                      <  1:784669/148‑1 (MQ=255)
                                                                                              atCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTa                                                     <  1:926149/148‑1 (MQ=255)
                                                                                               tCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTaa                                                    <  2:2353409/148‑1 (MQ=255)
                                                                                               tCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTaa                                                    <  1:3139434/148‑1 (MQ=255)
                                                                                                cTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAAc                                                   >  2:1155589/1‑148 (MQ=255)
                                                                                                    cgctcgctTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAg                                               <  1:3103256/148‑1 (MQ=255)
                                                                                                      ctcgctTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTc                                             <  1:127981/148‑1 (MQ=255)
                                                                                                        cgctTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGGGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCAt                                           >  2:142759/1‑148 (MQ=255)
                                                                                                        cgctTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCAt                                           >  2:894179/1‑148 (MQ=255)
                                                                                                             gagagagaTATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGAttt                                      <  1:29366/148‑1 (MQ=255)
                                                                                                              agagagaTATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTa                                     <  2:2412695/148‑1 (MQ=255)
                                                                                                               gagagaTATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTAt                                    >  2:1893570/1‑148 (MQ=255)
                                                                                                               gagagaTATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTAt                                    >  2:2163355/1‑148 (MQ=255)
                                                                                                                   gaTATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATgcgc                                >  1:1252154/1‑148 (MQ=255)
                                                                                                                    atatTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCa                               >  1:1603393/1‑148 (MQ=255)
                                                                                                                         ttCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAaca                          <  2:1146136/148‑1 (MQ=255)
                                                                                                                             cGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAgg                      <  2:12925/148‑1 (MQ=255)
                                                                                                                              ggCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGc                     <  1:2274001/148‑1 (MQ=255)
                                                                                                                               gCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGcc                    <  1:1686830/148‑1 (MQ=255)
                                                                                                                               gCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGcc                    >  2:2856152/1‑148 (MQ=255)
                                                                                                                                 aaTCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtg                  <  1:2160088/148‑1 (MQ=255)
                                                                                                                                  aTCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGAATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtga                 >  2:1601987/1‑148 (MQ=255)
                                                                                                                                  aTCGGTGATATCAGTATAACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtga                 >  1:3041079/1‑148 (MQ=255)
                                                                                                                                   tCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgat                >  1:3304196/1‑148 (MQ=255)
                                                                                                                                   tCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgat                <  1:284287/148‑1 (MQ=255)
                                                                                                                                   tCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgat                <  1:369283/148‑1 (MQ=255)
                                                                                                                                   tCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgat                >  2:2761032/1‑148 (MQ=255)
                                                                                                                                   tCGGTGATATCAGTATAACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgat                <  1:1421620/148‑1 (MQ=255)
                                                                                                                                    cGGTGATATCAGTATAACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgatt               >  1:3285828/1‑148 (MQ=255)
                                                                                                                                     ggTGTTATCAGAATCTATGAGATTAAGGCTCATACTGTTGTTAATGCTGTTAAGCCGGTTCAGGCCAGAGGTGCTTTAGATACTTTTCGCCGCCTTTGTCTGAGTATTAACGAAGTCATGATTTATGCTCAGAACACAGGCCtgattg              <  2:1392794/148‑1 (MQ=255)
                                                                                                                                     ggTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgattg              <  2:2661822/148‑1 (MQ=255)
                                                                                                                                      gTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgattga             <  1:590992/148‑1 (MQ=255)
                                                                                                                                      gTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgattga             <  2:1036738/148‑1 (MQ=255)
                                                                                                                                       tGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgattgat            >  1:1186836/1‑148 (MQ=255)
                                                                                                                                       tGATATCAGTATAACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgattgat            <  1:2022231/148‑1 (MQ=255)
                                                                                                                                       tGATATCAGTATAACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgattgat            <  1:1784588/148‑1 (MQ=255)
                                                                                                                                        gATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATg           <  1:2053781/148‑1 (MQ=255)
                                                                                                                                        gATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATg           <  2:261690/148‑1 (MQ=255)
                                                                                                                                           atCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTg        <  1:2595571/148‑1 (MQ=255)
                                                                                                                                            tCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGt       >  2:3037519/1‑148 (MQ=255)
                                                                                                                                               gTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTcc    <  1:429436/148‑1 (MQ=255)
                                                                                                                                                 aTCACTGAGATTAAGGCTCATACTCGGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATTCTGTTCCTa  <  2:2233239/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTCTCAAGAGGCCAAAACCAATACCTTCTTGTTAGTTGCCGAGCGTTGGTGGAATGTGAAGAAAACCAGCGTAACAGAGGACTATGCCGACGATATCTGGCGCTCGCTTGAGAGAGATATTTTCCCGGCAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTA  >  CP009273/2749811‑2750103

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: