Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,679,324 C→A 100% M1070I (ATG→ATT)  bcsC ← cellulose synthase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,679,3240CA99.4% 1059.1 / ‑2.2 309M1070I (ATG→ATT) bcsCcellulose synthase subunit
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (170/137);  total (170/139)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 322 total aligned reads displayed.

TGGATATCAATTGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGC  >  CP009273/3679178‑3679469
                                                                                                                                                  |                                                                                                                                                 
tgGATATCAATTGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATtaat                                                                                                                                                 >  2:3168341/1‑149 (MQ=255)
tgGATATCAATTGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATtaat                                                                                                                                                 >  1:2675250/1‑149 (MQ=255)
  gATATCAATTGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATtaatca                                                                                                                                               <  1:834081/149‑1 (MQ=255)
  gATATCAATTGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATtaatca                                                                                                                                               >  1:2721059/1‑149 (MQ=255)
   atatCAATTGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAg                                                                                                                                              <  2:2236491/149‑1 (MQ=255)
    tatCAATTGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGc                                                                                                                                             <  1:3160695/149‑1 (MQ=255)
    tatCAATTGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGc                                                                                                                                             <  2:421548/149‑1 (MQ=255)
     atCAATTGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCg                                                                                                                                            <  2:2946852/149‑1 (MQ=255)
       cAATTGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGa                                                                                                                                          <  1:2831194/149‑1 (MQ=255)
       cAATTGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGa                                                                                                                                          >  1:2067452/1‑149 (MQ=255)
           tGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAAc                                                                                                                                      >  1:245183/1‑149 (MQ=255)
           tGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAAc                                                                                                                                      <  1:408656/149‑1 (MQ=255)
             ccGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAAcgc                                                                                                                                    <  2:487275/149‑1 (MQ=255)
              cGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAAcgcg                                                                                                                                   <  1:3495606/149‑1 (MQ=255)
                tGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGc                                                                                                                                 >  1:3797974/1‑149 (MQ=255)
                  ccAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCAt                                                                                                                               >  1:1970037/1‑149 (MQ=255)
                   caacaaACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATg                                                                                                                              <  1:786020/149‑1 (MQ=255)
                   caacaaACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATg                                                                                                                              <  1:4681414/149‑1 (MQ=255)
                     acaaACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGGGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGt                                                                                                                            >  1:3207824/1‑149 (MQ=255)
                       aaaCCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGttt                                                                                                                          >  2:1450707/1‑149 (MQ=255)
                       aaaCCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGttt                                                                                                                          <  2:2497372/149‑1 (MQ=255)
                        aaCCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGtttt                                                                                                                         <  2:438673/149‑1 (MQ=255)
                        aaCCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGtttt                                                                                                                         >  2:2398007/1‑149 (MQ=255)
                           cAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGt                                                                                                                      <  2:91801/149‑1 (MQ=255)
                           cAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCAGCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGt                                                                                                                      >  2:1986389/1‑149 (MQ=255)
                            aGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTg                                                                                                                     >  1:4389980/1‑149 (MQ=255)
                                ggAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTg                                                                                                                 >  1:4252132/1‑149 (MQ=255)
                                ggAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTg                                                                                                                 <  1:291447/149‑1 (MQ=255)
                                ggAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTg                                                                                                                 <  1:2480500/149‑1 (MQ=255)
                                ggAAGTAACACGTCGTTCAAGTAATGACCGCGACGTGTAGCCGAAGCCCTGACTGCTGCCGCCATAGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGACGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTg                                                                                                                 >  1:2347157/1‑149 (MQ=255)
                                 gAAGTAACACTTCGTTCAAGTAATTCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGACGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTAGTTCGGGATCAGATTAATCAGCGGATAACACGTCATGGTTTTGGTGCGTGa                                                                                                                >  2:632416/1‑149 (MQ=255)
                                      aaCACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCg                                                                                                           >  2:2097369/1‑149 (MQ=255)
                                       acacGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGa                                                                                                          >  2:4620144/1‑149 (MQ=255)
                                        cacGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGAc                                                                                                         <  1:3045579/149‑1 (MQ=255)
                                        cacGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGAc                                                                                                         <  1:654636/149‑1 (MQ=255)
                                         acGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGAcc                                                                                                        <  2:2015351/149‑1 (MQ=255)
                                               ttCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGc                                                                                                  <  1:1036719/149‑1 (MQ=255)
                                               ttCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGc                                                                                                  >  2:2814213/1‑149 (MQ=255)
                                                tCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGcc                                                                                                 <  2:263637/149‑1 (MQ=255)
                                                   aGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGa                                                                                              <  2:1179067/149‑1 (MQ=255)
                                                    gTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGAc                                                                                             >  2:3679571/1‑149 (MQ=255)
                                                        tGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCAc                                                                                         <  1:591628/149‑1 (MQ=255)
                                                         gCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCAcc                                                                                        >  1:1185607/1‑149 (MQ=255)
                                                           ccGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCa                                                                                      <  2:1670101/149‑1 (MQ=255)
                                                           ccGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCa                                                                                      <  2:1918292/149‑1 (MQ=255)
                                                           ccGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCa                                                                                      >  1:1320457/1‑149 (MQ=255)
                                                            cgcgCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAg                                                                                     <  1:1616152/149‑1 (MQ=255)
                                                             gcgcCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGc                                                                                    >  1:1398788/1‑149 (MQ=255)
                                                               gcCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTc                                                                                  >  1:3913009/1‑149 (MQ=255)
                                                               gcCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTc                                                                                  <  2:2147600/149‑1 (MQ=255)
                                                                  gtgtAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCa                                                                               <  1:3633489/149‑1 (MQ=255)
                                                                   tgtAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCAc                                                                              >  2:1321210/1‑149 (MQ=255)
                                                                    gtAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACg                                                                             >  2:2868260/1‑149 (MQ=255)
                                                                       gCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGAcc                                                                          >  2:3002202/1‑149 (MQ=255)
                                                                       gCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGAcc                                                                          <  2:124012/149‑1 (MQ=255)
                                                                        ccGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCa                                                                         >  2:3334697/1‑149 (MQ=255)
                                                                          gAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGt                                                                       <  1:2890430/149‑1 (MQ=255)
                                                                          gAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGt                                                                       >  2:1426424/1‑149 (MQ=255)
                                                                          gAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGt                                                                       >  2:2697910/1‑149 (MQ=255)
                                                                           aaGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGtt                                                                      >  1:2250909/1‑149 (MQ=255)
                                                                           aaGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTAGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGtt                                                                      <  1:3882476/149‑1 (MQ=255)
                                                                            aGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGttt                                                                     <  2:346409/149‑1 (MQ=255)
                                                                            aGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGttt                                                                     >  1:933230/1‑149 (MQ=255)
                                                                            aGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGttt                                                                     >  1:1799399/1‑149 (MQ=255)
                                                                             gCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGtttt                                                                    >  2:313006/1‑149 (MQ=255)
                                                                              cccTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTc                                                                   <  2:2902559/149‑1 (MQ=255)
                                                                                  gACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTg                                                               <  1:3254332/149‑1 (MQ=255)
                                                                                     tgctgcCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGTGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTATTCCGTgcgc                                                            >  1:41808/1‑149 (MQ=255)
                                                                                       ctgcCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTc                                                          <  2:3508348/149‑1 (MQ=255)
                                                                                         gccgccATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTccc                                                        <  1:1667195/149‑1 (MQ=255)
                                                                                         gccgccATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTccc                                                        <  2:4165322/149‑1 (MQ=255)
                                                                                          ccgccATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCg                                                       >  1:1534178/1‑149 (MQ=255)
                                                                                           cgccATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGc                                                      <  1:891811/149‑1 (MQ=255)
                                                                                              cATCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCcac                                                   >  2:4625647/1‑149 (MQ=255)
                                                                                                tCGTTGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACAt                                                 >  1:2857801/1‑149 (MQ=255)
                                                                                                   ttGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCAc                                              >  2:4345919/1‑149 (MQ=255)
                                                                                                   ttGGATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCAc                                              >  2:4023841/1‑149 (MQ=255)
                                                                                                      gATTGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCgg                                           <  2:1522278/149‑1 (MQ=255)
                                                                                                        ttGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTa                                         <  1:3410262/149‑1 (MQ=255)
                                                                                                         tGCCGCGCAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTAt                                        <  1:193462/149‑1 (MQ=255)
                                                                                                             gcgcAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCa                                    >  1:2941968/1‑149 (MQ=255)
                                                                                                              cgcAGCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCaa                                   <  1:813501/149‑1 (MQ=255)
                                                                                                                  gCTTCTTCCGGCCAGTCGGTCGTGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGa                               <  1:4198287/149‑1 (MQ=255)
                                                                                                                  gCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGa                               >  1:4501840/1‑149 (MQ=255)
                                                                                                                  gCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGa                               >  1:1081459/1‑149 (MQ=255)
                                                                                                                  gCTTCTTCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACACACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGa                               >  1:3362078/1‑149 (MQ=255)
                                                                                                                   cttcttCCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGAc                              >  2:4499611/1‑149 (MQ=255)
                                                                                                                     tcttcCGGCCTGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAg                            <  2:3717770/149‑1 (MQ=255)
                                                                                                                       ttcCGGCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGt                          >  2:1023597/1‑149 (MQ=255)
                                                                                                                           ggCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTc                      >  1:2496533/2‑149 (MQ=255)
                                                                                                                            gCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTcc                     <  2:3912789/149‑1 (MQ=255)
                                                                                                                            gCCAGTCGGTCGGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTcc                     >  2:1806921/1‑149 (MQ=255)
                                                                                                                                gtcggtcgGGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCg                 >  2:2848206/1‑149 (MQ=255)
                                                                                                                                   ggtcggTATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGAc              <  2:1157998/149‑1 (MQ=255)
                                                                                                                                   ggtcggGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGAc              >  2:544768/1‑149 (MQ=255)
                                                                                                                                    gtcggGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACt             >  1:3997022/1‑149 (MQ=255)
                                                                                                                                    gtcggGATCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACt             >  1:4724586/1‑149 (MQ=255)
                                                                                                                                      cggTATCAGATTAATCAGCGGTTAACGCGGCTTGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGCCCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTgt           <  2:1587581/149‑1 (MQ=255)
                                                                                                                                          aTCAGATTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTgtagta       >  2:2325690/1‑149 (MQ=255)
                                                                                                                                               aTTAATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCGACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGc  >  2:1098059/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGGATATCAATTGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGC  >  CP009273/3679178‑3679469

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