Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,227,000 G→T 6.4% G120V (GGC→GTC)  umuC → DNA polymerase V, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,227,0000GT 6.4% 926.9 / 23.3 249G120V (GGC→GTC) umuCDNA polymerase V, subunit C
Reads supporting (aligned to +/- strand):  ref base G (109/124);  new base T (8/8);  total (117/132)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 253 total aligned reads displayed.

AGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTG  >  CP009273/1226852‑1227148
                                                                                                                                                    |                                                                                                                                                    
 gcagcaATTATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGc                                                                                                                                                     >  1:2253745/1‑149 (MQ=255)
 gcagcaATTATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGc                                                                                                                                                     >  2:4094240/1‑149 (MQ=255)
    gcaATTATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCaga                                                                                                                                                  >  1:130831/1‑149 (MQ=255)
    gcaATTATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCaga                                                                                                                                                  <  2:3400013/149‑1 (MQ=255)
      aaTTATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGGAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAgaga                                                                                                                                                <  1:2237863/149‑1 (MQ=255)
      aaTTATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAgaga                                                                                                                                                <  1:3124018/149‑1 (MQ=255)
      aaTTATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAgaga                                                                                                                                                <  1:918994/149‑1 (MQ=255)
      aaTTATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAgaga                                                                                                                                                <  1:164776/149‑1 (MQ=255)
        ttATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGaaa                                                                                                                                              <  2:1966855/149‑1 (MQ=255)
         tATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAAt                                                                                                                                             >  2:1100847/1‑149 (MQ=255)
          aTGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAAtt                                                                                                                                            <  1:3333896/149‑1 (MQ=255)
           tGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTc                                                                                                                                           >  2:1478834/1‑149 (MQ=255)
            gAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTcg                                                                                                                                          >  1:4059321/1‑149 (MQ=255)
             aGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTcgc                                                                                                                                         >  1:2596854/1‑149 (MQ=255)
                 ttACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCaa                                                                                                                                     <  1:1040256/149‑1 (MQ=255)
                    cGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACgg                                                                                                                                  <  2:544015/149‑1 (MQ=255)
                    cGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACgg                                                                                                                                  <  2:3337515/149‑1 (MQ=255)
                     gCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGt                                                                                                                                 >  2:407261/1‑149 (MQ=255)
                       aGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGc                                                                                                                               <  2:3009146/149‑1 (MQ=255)
                        gACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCt                                                                                                                              >  1:3153164/1‑149 (MQ=255)
                         aCATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTa                                                                                                                             <  1:1055610/149‑1 (MQ=255)
                           aTGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACa                                                                                                                           >  2:2855175/1‑149 (MQ=255)
                             gAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAAc                                                                                                                         <  2:1728859/149‑1 (MQ=255)
                               gCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGt                                                                                                                       >  1:110131/1‑149 (MQ=255)
                                 aaTCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTAc                                                                                                                     <  2:3080030/149‑1 (MQ=255)
                                 aaTCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTAc                                                                                                                     <  2:1021959/149‑1 (MQ=255)
                                  aTCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTAcc                                                                                                                    <  2:1866151/149‑1 (MQ=255)
                                  aTCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTAcc                                                                                                                    >  1:870998/1‑149 (MQ=255)
                                  aTCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTAcc                                                                                                                    <  1:1450863/149‑1 (MQ=255)
                                     gggTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCAt                                                                                                                 >  2:59513/1‑149 (MQ=255)
                                         gATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTa                                                                                                             <  1:2366820/149‑1 (MQ=255)
                                          aTGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGTCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTAc                                                                                                            <  1:2862946/149‑1 (MQ=255)
                                                aCGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTgg                                                                                                      >  2:3990970/1‑149 (MQ=255)
                                                  gCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGgtg                                                                                                    <  1:2933066/149‑1 (MQ=255)
                                                      gAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACCACGTACCCCTCTTACGGGGGGTGTggg                                                                                                >  2:347029/1‑149 (MQ=255)
                                                      gAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTggg                                                                                                <  2:3734150/149‑1 (MQ=255)
                                                      gAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTggg                                                                                                >  2:346104/1‑149 (MQ=255)
                                                      gAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTggg                                                                                                >  1:300948/1‑149 (MQ=255)
                                                        agagCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGa                                                                                              >  1:2216861/1‑149 (MQ=255)
                                                        agagCTATCGCCCCGCGTCGAGATTTACAGAATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGTGGa                                                                                              >  1:974727/1‑149 (MQ=255)
                                                         gagCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGAt                                                                                             <  2:2459514/149‑1 (MQ=255)
                                                          agCTATCGCCCCGCGTCGGGATTTACAGTATTGATGAGGCATTATGCGAGCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATc                                                                                            >  1:923774/1‑149 (MQ=255)
                                                          agCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATc                                                                                            >  2:3905504/1‑149 (MQ=255)
                                                          agCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATc                                                                                            <  2:913891/149‑1 (MQ=255)
                                                             tATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGcc                                                                                         >  2:1699809/1‑149 (MQ=255)
                                                              aTCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGTCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGccc                                                                                        <  1:2897516/149‑1 (MQ=255)
                                                              aTCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGccc                                                                                        >  1:269653/1‑149 (MQ=255)
                                                               tCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCa                                                                                       <  2:2253754/149‑1 (MQ=255)
                                                               tCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCa                                                                                       >  1:205539/1‑149 (MQ=255)
                                                                cGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAg                                                                                      >  1:1665151/1‑149 (MQ=255)
                                                                   cccGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGAcc                                                                                   <  2:671900/149‑1 (MQ=255)
                                                                    ccGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCa                                                                                  >  1:3817762/1‑149 (MQ=255)
                                                                    ccGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCa                                                                                  >  1:635483/1‑149 (MQ=255)
                                                                       cgTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCaaaa                                                                               <  1:2417033/149‑1 (MQ=255)
                                                                           gagaTTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCt                                                                           >  2:955545/1‑149 (MQ=255)
                                                                           gagaTTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCt                                                                           >  1:1690919/1‑149 (MQ=255)
                                                                             gaTTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTgg                                                                         >  1:364931/1‑149 (MQ=255)
                                                                              aTTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGc                                                                        >  2:415677/1‑149 (MQ=255)
                                                                               tttACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCt                                                                       >  1:1993657/1‑149 (MQ=255)
                                                                                    aGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCt                                                                  <  1:2605847/149‑1 (MQ=255)
                                                                                       aTTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGc                                                               <  2:3327308/149‑1 (MQ=255)
                                                                                        ttGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGcc                                                              <  1:4102362/149‑1 (MQ=255)
                                                                                          gatgaGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCaa                                                            <  2:2738809/149‑1 (MQ=255)
                                                                                          gatgaGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCaa                                                            >  2:3524332/1‑149 (MQ=255)
                                                                                            tgaGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATc                                                          <  1:4094201/149‑1 (MQ=255)
                                                                                            tgaGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATc                                                          <  2:2877735/149‑1 (MQ=255)
                                                                                              aGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCAt                                                        <  2:710603/149‑1 (MQ=255)
                                                                                               ggCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATg                                                       >  1:2588056/1‑149 (MQ=255)
                                                                                                gCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGTCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGc                                                      >  1:1358079/1‑149 (MQ=255)
                                                                                                   ttCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGc                                                   >  2:2695243/1‑149 (MQ=255)
                                                                                                   ttCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGc                                                   >  2:1811837/1‑149 (MQ=255)
                                                                                                       gCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCaaaa                                               >  2:1194571/1‑149 (MQ=255)
                                                                                                          aTCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCaaaaaaa                                            >  1:3200895/1‑149 (MQ=255)
                                                                                                          aTCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCaaaaaaa                                            >  2:881823/1‑149 (MQ=255)
                                                                                                          aTCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCaaaaaaa                                            <  1:3979624/149‑1 (MQ=255)
                                                                                                           tCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGTCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAAt                                           <  2:236754/149‑1 (MQ=255)
                                                                                                              gACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGc                                        >  1:184402/1‑149 (MQ=255)
                                                                                                                   gtgtGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCgg                                   >  2:3931511/1‑149 (MQ=255)
                                                                                                                     gtgCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCa                                 <  1:1727760/149‑1 (MQ=255)
                                                                                                                      tgCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAg                                <  1:1100826/149‑1 (MQ=255)
                                                                                                                        cGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGAc                              <  2:1039338/149‑1 (MQ=255)
                                                                                                                            aTTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACgggt                          >  2:3977802/1‑149 (MQ=255)
                                                                                                                               gTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACgggtggg                       >  2:1408408/1‑149 (MQ=255)
                                                                                                                                  gcgATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACgggtggggtg                    >  1:2287759/1‑149 (MQ=255)
                                                                                                                                     aTCTGACTGATTTTGTCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGgtggtg                 >  2:3659229/1‑149 (MQ=255)
                                                                                                                                     aTCTGACTGATTTTGTCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGgtggtg                 >  2:1616686/1‑149 (MQ=255)
                                                                                                                                     aTCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGgtggtg                 <  2:4035852/149‑1 (MQ=255)
                                                                                                                                       ctgactgaTTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGa               <  2:3767659/149‑1 (MQ=255)
                                                                                                                                        tgactgaTTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGTGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGAt              <  1:1909023/149‑1 (MQ=255)
                                                                                                                                        tgactgaTTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGTGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGAt              <  1:1464739/149‑1 (MQ=255)
                                                                                                                                          actgaTTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGAttt            >  2:1267758/1‑149 (MQ=255)
                                                                                                                                           ctgaTTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTa           <  2:2497225/149‑1 (MQ=255)
                                                                                                                                           ctgaTTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTa           <  1:2458017/149‑1 (MQ=255)
                                                                                                                                            tgaTTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTAt          <  1:3676259/149‑1 (MQ=255)
                                                                                                                                             gaTTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCAGGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATc         <  2:357226/149‑1 (MQ=255)
                                                                                                                                               ttttGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCaa       <  1:3286568/149‑1 (MQ=255)
                                                                                                                                                tttGTCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCaaa      <  2:3423426/149‑1 (MQ=255)
                                                                                                                                                tttGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCACCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCaaa      <  2:2082884/149‑1 (MQ=255)
                                                                                                                                                  tGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGTGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATc    <  1:388261/149‑1 (MQ=255)
                                                                                                                                                    gCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTg  <  2:1320807/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGCAGCAATTATGAGCTTTACGCAGACATGAGCAATCGGGTGATGTCGACGCTGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTG  >  CP009273/1226852‑1227148

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: