Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,953,398 C→A 11.3% G797V (GGC→GTC)  recC ← exonuclease V (RecBCD complex), gamma chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,953,3980CA11.3% 1086.3 / 86.4 300G797V (GGC→GTC) recCexonuclease V (RecBCD complex), gamma chain
Reads supporting (aligned to +/- strand):  ref base C (117/149);  new base A (20/14);  total (137/163)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 307 total aligned reads displayed.

CCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCAC  >  CP009273/2953250‑2953546
                                                                                                                                                    |                                                                                                                                                    
   cACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCctg                                                                                                                                                   <  2:4462554/149‑1 (MQ=255)
    aCCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGActgg                                                                                                                                                  >  1:2777904/1‑149 (MQ=255)
        gATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCctggctgg                                                                                                                                              >  2:2527754/1‑149 (MQ=255)
         aTGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGt                                                                                                                                             >  1:260547/1‑149 (MQ=255)
         aTGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGACTGGCTGGt                                                                                                                                             >  2:10756/1‑149 (MQ=255)
             gCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTATGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGACTGGCTGGTAGtt                                                                                                                                         <  1:1316688/149‑1 (MQ=255)
             gCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGtt                                                                                                                                         >  2:1969509/1‑149 (MQ=255)
             gCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGACTGGCTGGTAGtt                                                                                                                                         >  2:1725940/1‑149 (MQ=255)
                cAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTtg                                                                                                                                      <  1:4676293/149‑1 (MQ=255)
                   aaTCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTgg                                                                                                                                   <  2:2154591/149‑1 (MQ=255)
                   aaTCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTgg                                                                                                                                   <  2:3856637/149‑1 (MQ=255)
                    aTCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGa                                                                                                                                  <  1:2087973/149‑1 (MQ=255)
                    aTCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGa                                                                                                                                  <  1:3978958/149‑1 (MQ=255)
                     tCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGAt                                                                                                                                 <  1:3758785/149‑1 (MQ=255)
                     tCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGACTGGCTGGTAGTTTTGTGGAt                                                                                                                                 >  1:3198035/1‑149 (MQ=255)
                       gTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCa                                                                                                                               <  1:1597287/149‑1 (MQ=255)
                         tGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCaaa                                                                                                                             <  2:607077/149‑1 (MQ=255)
                                tttCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATc                                                                                                                      <  1:2808785/149‑1 (MQ=255)
                                tttCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATc                                                                                                                      <  2:2651634/149‑1 (MQ=255)
                                   ccAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCgg                                                                                                                   <  1:4409049/149‑1 (MQ=255)
                                   ccAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCgg                                                                                                                   >  1:4024837/1‑149 (MQ=255)
                                   ccAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCgg                                                                                                                   >  2:235872/1‑149 (MQ=255)
                                   ccAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCgg                                                                                                                   <  1:2217587/149‑1 (MQ=255)
                                    cAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCggg                                                                                                                  <  2:1281576/149‑1 (MQ=255)
                                     aGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGt                                                                                                                 <  2:5117368/149‑1 (MQ=255)
                                        ggCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATg                                                                                                              >  1:65775/1‑149 (MQ=255)
                                        ggCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATg                                                                                                              <  2:2496034/149‑1 (MQ=255)
                                         gCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATgg                                                                                                             >  1:287905/1‑149 (MQ=255)
                                          cACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGa                                                                                                            >  2:1158743/1‑149 (MQ=255)
                                          cACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGa                                                                                                            <  1:4757267/149‑1 (MQ=255)
                                           aCGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAg                                                                                                           <  1:2917967/149‑1 (MQ=255)
                                           aCGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGACTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAg                                                                                                           <  2:5040576/149‑1 (MQ=255)
                                            cGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAgg                                                                                                          >  1:136778/1‑149 (MQ=255)
                                               tttCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCaa                                                                                                       <  1:228738/149‑1 (MQ=255)
                                               tttCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCaa                                                                                                       <  2:2028699/149‑1 (MQ=255)
                                                ttCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAg                                                                                                      <  1:3909059/149‑1 (MQ=255)
                                                       aaggtaaACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATg                                                                                               >  2:2815648/1‑149 (MQ=255)
                                                       aaggtaaACGGCAGCGGCTGAACAAATTCAGAATGTGCTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATgc                                                                                              <  1:949066/149‑1 (MQ=255)
                                                             aaCGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCttt                                                                                         <  1:4840095/149‑1 (MQ=255)
                                                             aaCGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGACTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCttt                                                                                         >  1:731449/1‑149 (MQ=255)
                                                                ggCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTAc                                                                                      <  2:4349299/149‑1 (MQ=255)
                                                                   aGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCt                                                                                   >  1:798120/1‑149 (MQ=255)
                                                                         tGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTc                                                                             >  1:5136911/1‑149 (MQ=255)
                                                                         tGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTc                                                                             >  2:2220201/1‑149 (MQ=255)
                                                                           aaCAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGc                                                                           <  2:4664694/149‑1 (MQ=255)
                                                                            aCAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCt                                                                          <  1:3694852/149‑1 (MQ=255)
                                                                             cAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCtt                                                                         <  1:4442492/149‑1 (MQ=255)
                                                                                aTTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTtca                                                                      <  1:1526027/149‑1 (MQ=255)
                                                                                 ttCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTtcat                                                                     <  2:2272913/149‑1 (MQ=255)
                                                                                  tCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTtcatc                                                                    >  1:2121238/1‑149 (MQ=255)
                                                                                     gAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATcaca                                                                 >  1:3816513/1‑149 (MQ=255)
                                                                                      aaTGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAg                                                                <  1:2713757/149‑1 (MQ=255)
                                                                                       aTGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGt                                                               >  1:492679/1‑149 (MQ=255)
                                                                                       aTGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGACTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGt                                                               >  1:3067505/1‑149 (MQ=255)
                                                                                        tgtgCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGtt                                                              >  2:2887688/1‑149 (MQ=255)
                                                                                        tgtgCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGtt                                                              >  2:2441640/1‑149 (MQ=255)
                                                                                        tgtgCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGtt                                                              >  2:3193299/1‑149 (MQ=255)
                                                                                          tgCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGa                                                            <  1:3360483/149‑1 (MQ=255)
                                                                                           gCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAg                                                           <  2:5140285/149‑1 (MQ=255)
                                                                                            cTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAgc                                                          <  2:4147524/149‑1 (MQ=255)
                                                                                            cTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAgc                                                          <  2:1876520/149‑1 (MQ=255)
                                                                                            cTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAgc                                                          >  1:1639057/1‑149 (MQ=255)
                                                                                             tttACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAgcg                                                         >  2:3754166/1‑149 (MQ=255)
                                                                                             tttACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAgcg                                                         >  2:3265334/1‑149 (MQ=255)
                                                                                              ttACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAgcgc                                                        <  2:666553/149‑1 (MQ=255)
                                                                                              ttACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGACTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAgcgc                                                        <  2:1377727/149‑1 (MQ=255)
                                                                                               tACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCt                                                       >  1:2319049/1‑149 (MQ=255)
                                                                                                 ccAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTc                                                     >  1:1115527/1‑149 (MQ=255)
                                                                                                 ccAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTc                                                     >  1:5050595/1‑149 (MQ=255)
                                                                                                   aGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCAt                                                   <  2:3362567/149‑1 (MQ=255)
                                                                                                    gCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATc                                                  >  1:1007007/1‑149 (MQ=255)
                                                                                                    gCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATc                                                  >  1:2836213/1‑149 (MQ=255)
                                                                                                    gCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATc                                                  <  1:2985621/149‑1 (MQ=255)
                                                                                                    gCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGACTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATc                                                  >  2:3580042/1‑149 (MQ=255)
                                                                                                     ccTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCg                                                 <  1:2952447/149‑1 (MQ=255)
                                                                                                      ctggctggCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGc                                                >  1:1888269/1‑149 (MQ=255)
                                                                                                      ctggctggCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGACTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGc                                                <  1:4127122/149‑1 (MQ=255)
                                                                                                       tggctggcCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGTGTATGGAGGCAAGTAAGATGCGCTTTTACCCTAGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGcc                                               >  2:2089937/1‑149 (MQ=255)
                                                                                                        ggctggcCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGccc                                              >  1:4524849/1‑149 (MQ=255)
                                                                                                          ctggcCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCgg                                            <  1:1194777/149‑1 (MQ=255)
                                                                                                            ggcCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGGTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTa                                          <  2:867463/149‑1 (MQ=255)
                                                                                                             gcCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAg                                         >  2:1461359/1‑149 (MQ=255)
                                                                                                             gcCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAg                                         <  2:2988455/149‑1 (MQ=255)
                                                                                                                 cAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATaa                                     <  1:4966723/149‑1 (MQ=255)
                                                                                                                  aGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAAt                                    <  1:2527754/149‑1 (MQ=255)
                                                                                                                   ggTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATg                                   <  2:287905/149‑1 (MQ=255)
                                                                                                                       gCCCTTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACtt                               >  2:4882434/1‑149 (MQ=255)
                                                                                                                               cGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGAt                       <  2:4589495/149‑1 (MQ=255)
                                                                                                                                 aGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGt                     >  1:2362946/1‑149 (MQ=255)
                                                                                                                                 aGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGt                     <  2:5173164/149‑1 (MQ=255)
                                                                                                                                   cATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAg                   >  2:3655591/1‑149 (MQ=255)
                                                                                                                                     tAGCTTTGTCGTTCGACTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTc                 >  1:875892/1‑149 (MQ=255)
                                                                                                                                        cTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGAt              >  1:4266997/1‑149 (MQ=255)
                                                                                                                                        cTTTGTCGTTCGACTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGAt              >  1:603268/1‑149 (MQ=255)
                                                                                                                                          ttGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCa            <  1:460008/149‑1 (MQ=255)
                                                                                                                                          ttGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCa            <  2:340795/149‑1 (MQ=255)
                                                                                                                                           tGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAg           <  1:544890/149‑1 (MQ=255)
                                                                                                                                            gTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGt          <  2:4199831/149‑1 (MQ=255)
                                                                                                                                            gTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGt          >  2:4777184/1‑149 (MQ=255)
                                                                                                                                                  cGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGc    <  1:3545736/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCAC  >  CP009273/2953250‑2953546

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: