Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,070,058 C→A 55.0% A90A (GCC→GCA yeeS → CP4‑44 prophage; putative DNA repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,070,0580CA55.0% 27.2 / inf 351A90A (GCC→GCAyeeSCP4‑44 prophage; putative DNA repair protein
Reads supporting (aligned to +/- strand):  ref base C (96/62);  new base A (111/81);  total (208/143)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 352 total aligned reads displayed.

TGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGTGGCAACCAGG  >  CP009273/2069910‑2070203
                                                                                                                                                    |                                                                                                                                                 
    gggACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTg                                                                                                                                               >  1:4834698/1‑149 (MQ=255)
     ggACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGc                                                                                                                                              <  2:1658393/149‑1 (MQ=255)
     ggACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGc                                                                                                                                              <  2:3205172/149‑1 (MQ=255)
        cTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGt                                                                                                                                           <  1:2786899/149‑1 (MQ=255)
        cTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGt                                                                                                                                           <  2:403564/149‑1 (MQ=255)
          ggAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTgg                                                                                                                                         >  2:1593333/1‑149 (MQ=255)
           gAGCGTGAAGACTTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGt                                                                                                                                        <  2:3509673/149‑1 (MQ=255)
             gCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTAc                                                                                                                                      <  1:1653302/149‑1 (MQ=255)
                 gaagaaTTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGc                                                                                                                                  >  2:2781144/1‑149 (MQ=255)
                  aagaaTTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGca                                                                                                                                 >  1:4475656/1‑149 (MQ=255)
                  aagaaTTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGca                                                                                                                                 >  2:2012677/1‑149 (MQ=255)
                    gaaTTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGcaca                                                                                                                               >  1:5202782/1‑149 (MQ=255)
                     aaTTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGcacac                                                                                                                              >  2:1090905/1‑149 (MQ=255)
                      aTTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGCAGTGAGTAATCGCGCCCTGTACCACACTGCCGCTGCCGTGGTACTGGcacaca                                                                                                                             >  2:3075328/1‑149 (MQ=255)
                      aTTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGcacaca                                                                                                                             >  1:2608086/1‑149 (MQ=255)
                         ccGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATc                                                                                                                          >  2:2927218/1‑149 (MQ=255)
                         ccGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATc                                                                                                                          >  1:1099741/1‑149 (MQ=255)
                          cGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCa                                                                                                                         >  1:624379/1‑149 (MQ=255)
                           gggTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTTCCACAATGCCGCTGCCGTGGTACTGGCCCACAAGCAc                                                                                                                        >  1:1905588/1‑149 (MQ=255)
                           gggTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCAc                                                                                                                        <  1:391917/149‑1 (MQ=255)
                                 tGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTcc                                                                                                                  <  1:2550948/149‑1 (MQ=255)
                                 tGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTcc                                                                                                                  <  2:8081/149‑1 (MQ=255)
                                    aTCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGt                                                                                                               >  2:2514666/1‑149 (MQ=255)
                                    aTCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGt                                                                                                               >  1:2886103/1‑149 (MQ=255)
                                    aTCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGt                                                                                                               <  1:2960807/149‑1 (MQ=255)
                                     tCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTg                                                                                                              <  2:4326270/149‑1 (MQ=255)
                                     tCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTg                                                                                                              >  1:4784369/1‑149 (MQ=255)
                                       tGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGaa                                                                                                            <  1:3762617/149‑1 (MQ=255)
                                         aacaacCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGACGt                                                                                                          >  2:2643602/1‑149 (MQ=255)
                                         aacaacCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACAGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACACTCACCCGTCCGGTGAAGt                                                                                                          >  2:2991604/1‑149 (MQ=255)
                                          acccccAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTc                                                                                                         <  1:3597017/145‑1 (MQ=255)
                                          acaacCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTc                                                                                                         >  2:3868506/1‑149 (MQ=255)
                                          acaacCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTc                                                                                                         <  2:1810875/149‑1 (MQ=255)
                                              ccAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCacac                                                                                                     >  2:4787516/1‑149 (MQ=255)
                                                 gAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCa                                                                                                  >  1:1283098/1‑149 (MQ=255)
                                                  aaTCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAg                                                                                                 >  1:1410133/1‑149 (MQ=255)
                                                   aTCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGt                                                                                                >  2:2509571/1‑149 (MQ=255)
                                                    tCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTa                                                                                               >  1:2267514/1‑149 (MQ=255)
                                                     cAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTaa                                                                                              >  2:5231635/1‑149 (MQ=255)
                                                     cAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTaa                                                                                              <  1:307195/149‑1 (MQ=255)
                                                      aGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAg                                                                                             >  2:2365200/1‑149 (MQ=255)
                                                       gCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAgg                                                                                            >  2:4560505/1‑149 (MQ=255)
                                                          gATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAg                                                                                         <  2:2133274/149‑1 (MQ=255)
                                                          gATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAg                                                                                         >  1:1511522/1‑149 (MQ=255)
                                                           aTTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGa                                                                                        <  2:519528/149‑1 (MQ=255)
                                                              gCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCg                                                                                     <  2:2697961/149‑1 (MQ=255)
                                                                 ggTGAAACCCTCTTCACCGGCACCATCAACCGCACGGACGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCt                                                                                  >  2:1047664/1‑149 (MQ=255)
                                                                 ggTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCt                                                                                  <  1:3251667/149‑1 (MQ=255)
                                                                  gTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCtt                                                                                 <  1:2010577/149‑1 (MQ=255)
                                                                    gAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTAt                                                                               >  1:4670683/1‑149 (MQ=255)
                                                                     aaaCCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATc                                                                              <  2:5075067/149‑1 (MQ=255)
                                                                     aaaCCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATc                                                                              >  2:1960157/1‑149 (MQ=255)
                                                                     aaaCCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATc                                                                              >  2:1933991/1‑149 (MQ=255)
                                                                       aCCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCAc                                                                            >  2:3890075/1‑149 (MQ=255)
                                                                       aCCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCAc                                                                            >  1:1738261/1‑149 (MQ=255)
                                                                        cccTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCAcc                                                                           >  1:5043421/1‑149 (MQ=255)
                                                                           tctTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGaa                                                                        >  1:861357/1‑149 (MQ=255)
                                                                           tctTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGaa                                                                        <  2:1826407/149‑1 (MQ=255)
                                                                                aCCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGACGTCACACCCAGTAAGGCAGACCGGCTAATCACCTAACctct                                                                   >  2:3709250/1‑149 (MQ=255)
                                                                                aCCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACgaaa                                                                   >  2:149771/1‑146 (MQ=255)
                                                                                aCCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCt                                                                   <  2:93548/149‑1 (MQ=255)
                                                                                 ccGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTg                                                                  <  1:797031/149‑1 (MQ=255)
                                                                                 ccGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTg                                                                  >  2:4891074/1‑149 (MQ=255)
                                                                                 ccGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTg                                                                  >  1:4924056/1‑149 (MQ=255)
                                                                                    gCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTa                                                               >  1:3999175/1‑149 (MQ=255)
                                                                                      aCCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACa                                                             >  1:1466861/1‑149 (MQ=255)
                                                                                            aaCCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACt                                                       >  2:4494377/1‑149 (MQ=255)
                                                                                                gCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGc                                                   <  2:624978/149‑1 (MQ=255)
                                                                                                    ggAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCtgg                                               >  1:4078009/1‑149 (MQ=255)
                                                                                                     gAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCtggt                                              <  2:4681496/149‑1 (MQ=255)
                                                                                                       aGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCtggtgg                                            >  2:719557/1‑149 (MQ=255)
                                                                                                       aGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCtggtgg                                            >  1:3147201/1‑149 (MQ=255)
                                                                                                       aGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCtggtgg                                            >  1:2696413/1‑149 (MQ=255)
                                                                                                           cATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGatat                                        >  2:4394122/1‑149 (MQ=255)
                                                                                                            aTCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATc                                       <  2:1904619/149‑1 (MQ=255)
                                                                                                            aTCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATc                                       >  2:3834315/1‑149 (MQ=255)
                                                                                                              ccccGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCg                                     <  2:4998980/149‑1 (MQ=255)
                                                                                                               cccGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCgg                                    <  2:1375498/149‑1 (MQ=255)
                                                                                                                ccGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCggg                                   >  1:2466881/1‑149 (MQ=255)
                                                                                                                ccGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCggg                                   <  2:324932/149‑1 (MQ=255)
                                                                                                                   ggAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGc                                <  2:2954657/149‑1 (MQ=255)
                                                                                                                      aGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTGAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCgg                             >  1:3472235/1‑149 (MQ=255)
                                                                                                                       gTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGa                            >  1:4294947/1‑149 (MQ=255)
                                                                                                                           ttAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCAt                        >  1:4972743/1‑149 (MQ=255)
                                                                                                                           ttAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCAt                        >  2:1564463/1‑149 (MQ=255)
                                                                                                                             aaaCGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCt                      >  2:1732860/1‑149 (MQ=255)
                                                                                                                              aaCGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTg                     >  1:1802728/1‑149 (MQ=255)
                                                                                                                                  cgcCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAg                 <  2:2195930/149‑1 (MQ=255)
                                                                                                                                  cgcCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAg                 <  1:100762/149‑1 (MQ=255)
                                                                                                                                   gcCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGt                <  2:3361835/149‑1 (MQ=255)
                                                                                                                                     ccTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCg              <  1:525327/149‑1 (MQ=255)
                                                                                                                                      cTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCgg             <  1:3237197/149‑1 (MQ=255)
                                                                                                                                       tGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGt            >  2:342109/1‑149 (MQ=255)
                                                                                                                                       tGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGt            <  1:2393940/149‑1 (MQ=255)
                                                                                                                                          aCCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGTGGc         >  1:3418083/1‑149 (MQ=255)
                                                                                                                                          aCCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGTGGc         <  2:2306311/149‑1 (MQ=255)
                                                                                                                                          aCCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGTGGc         <  1:169168/149‑1 (MQ=255)
                                                                                                                                            cacaATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGTGGCaa       >  2:282533/1‑149 (MQ=255)
                                                                                                                                              caATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGTGGCAAcc     >  2:1266658/1‑149 (MQ=255)
                                                                                                                                               aaTGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGTGGCAACCa    <  2:4096279/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGTGGCAACCAGG  >  CP009273/2069910‑2070203

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: