Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,881,645 C→A 100% D202Y (GAT→TAT)  rnd ← ribonuclease D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,881,6450CA98.9% 1136.3 / 1.9 354D202Y (GAT→TAT) rndribonuclease D
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base A (165/185);  total (166/188)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-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 360 total aligned reads displayed.

TATAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATCCGCCG  >  CP009273/1881498‑1881791
                                                                                                                                                   |                                                                                                                                                  
tataACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAAc                                                                                                                                                  >  2:3137180/1‑150 (MQ=255)
 ataACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACg                                                                                                                                                 >  2:1175035/1‑150 (MQ=255)
  taACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGc                                                                                                                                                >  2:1011545/1‑150 (MQ=255)
   aaCGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGcc                                                                                                                                               >  1:1157163/1‑150 (MQ=255)
    aCGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCa                                                                                                                                              <  2:3903367/150‑1 (MQ=255)
       cgcTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGc                                                                                                                                           >  1:4393597/1‑150 (MQ=255)
         cTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCAt                                                                                                                                         <  1:3681795/150‑1 (MQ=255)
           aCCGACCACACATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCTCCAGTCGGCTAACAGTTGCTGACAGGCCAGTTGGCGTGTGAGTAATTGCCAGGCTTTGGTGATATAACGCCAGGcttgt                                                                                                                                       <  2:1736213/150‑5 (MQ=255)
             cGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTc                                                                                                                                     >  2:1303465/1‑150 (MQ=255)
              gACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTcc                                                                                                                                    <  2:2915897/150‑1 (MQ=255)
                ccACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCgg                                                                                                                                  >  1:887879/1‑150 (MQ=255)
                 cacaAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGgc                                                                                                                                 >  1:233260/1‑150 (MQ=255)
                   caAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGgcgc                                                                                                                               >  1:1013820/1‑150 (MQ=255)
                     aaTGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCaa                                                                                                                             <  2:1605675/150‑1 (MQ=255)
                     aaTGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCaa                                                                                                                             <  2:844816/150‑1 (MQ=255)
                     aaTGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCaa                                                                                                                             <  1:3249854/150‑1 (MQ=255)
                     aaTGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCaa                                                                                                                             >  2:233821/1‑150 (MQ=255)
                              cacgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCt                                                                                                                    <  1:489278/150‑1 (MQ=255)
                              cacgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCt                                                                                                                    >  1:2819613/1‑150 (MQ=255)
                               acgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTg                                                                                                                   <  1:4999852/150‑1 (MQ=255)
                                  cacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTgacg                                                                                                                >  2:4785327/1‑150 (MQ=255)
                                  cacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTgacg                                                                                                                <  2:4123933/150‑1 (MQ=255)
                                   acgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTCACAGTTGCAGACCGGCCAGTTGGCGTGTGCGTAATTGCCAGGCCGTGGTGATATAACGCCCGGCATCTTCCGGCGCACCGACTTCCTgacga                                                                                                               >  2:4005621/1‑150 (MQ=255)
                                       cAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCa                                                                                                           >  2:643201/1‑150 (MQ=255)
                                       cAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCa                                                                                                           >  2:1488603/1‑150 (MQ=255)
                                        aaaGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCAt                                                                                                          <  2:750970/150‑1 (MQ=255)
                                        aaaGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCAt                                                                                                          >  2:4521465/1‑150 (MQ=255)
                                         aaGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCAtt                                                                                                         >  1:2519970/1‑150 (MQ=255)
                                         aaGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCAtt                                                                                                         <  2:1785616/150‑1 (MQ=255)
                                             tCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCa                                                                                                     >  1:1895461/1‑150 (MQ=255)
                                              cACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCAt                                                                                                    >  1:1574647/1‑150 (MQ=255)
                                               aCCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATc                                                                                                   >  1:3086340/1‑150 (MQ=255)
                                                 cGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAg                                                                                                 <  1:2159704/150‑1 (MQ=255)
                                                    cAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAggcg                                                                                              <  1:4108965/150‑1 (MQ=255)
                                                     aGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAggcgg                                                                                             >  2:3595417/1‑150 (MQ=255)
                                                      gATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAggcggc                                                                                            <  2:2989248/150‑1 (MQ=255)
                                                         cgcgCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCAtt                                                                                         <  1:207619/150‑1 (MQ=255)
                                                         cgcaCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCAtt                                                                                         <  2:199729/150‑1 (MQ=255)
                                                             ctctCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATc                                                                                     >  1:1060452/1‑150 (MQ=255)
                                                             ctctCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATc                                                                                     >  1:2591652/1‑150 (MQ=255)
                                                              tctcGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATcc                                                                                    <  2:4849984/150‑1 (MQ=255)
                                                                 cgcgCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAgc                                                                                 <  2:4656594/150‑1 (MQ=255)
                                                                   cgcCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAgcgc                                                                               >  1:457497/1‑150 (MQ=255)
                                                                   cgcCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAgcgc                                                                               <  1:3567373/150‑1 (MQ=255)
                                                                    gcCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGcgcc                                                                              >  2:3664253/1‑150 (MQ=255)
                                                                     ccTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGcgccg                                                                             <  2:2228929/150‑1 (MQ=255)
                                                                     ccTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGcgccg                                                                             <  1:1728551/150‑1 (MQ=255)
                                                                      cTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGcgccgc                                                                            >  2:2219637/1‑150 (MQ=255)
                                                                      cTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGcgccgc                                                                            >  2:1084967/1‑150 (MQ=255)
                                                                           tcAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTa                                                                       >  2:5110654/2‑150 (MQ=255)
                                                                           gcAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTa                                                                       <  2:2653899/150‑1 (MQ=255)
                                                                             aGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTagc                                                                     >  1:1802947/1‑150 (MQ=255)
                                                                             aGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTagc                                                                     >  1:1397659/1‑150 (MQ=255)
                                                                               tCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTagcca                                                                   <  1:4358166/150‑1 (MQ=255)
                                                                                    aGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCgg                                                              >  1:2994479/1‑150 (MQ=255)
                                                                                     gTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGa                                                             <  1:4505495/150‑1 (MQ=255)
                                                                                     gTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGa                                                             <  1:4092697/150‑1 (MQ=255)
                                                                                     gTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGa                                                             <  1:1525755/150‑1 (MQ=255)
                                                                                      tCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAg                                                            >  2:189415/1‑150 (MQ=255)
                                                                                       cGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAgg                                                           >  2:3418462/1‑150 (MQ=255)
                                                                                           tAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTc                                                       >  1:403957/1‑150 (MQ=255)
                                                                                            aaCAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTcc                                                      >  2:4283703/1‑150 (MQ=255)
                                                                                            aaCAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTcc                                                      >  2:963683/1‑150 (MQ=255)
                                                                                                 ttGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTc                                                 <  2:4353185/150‑1 (MQ=255)
                                                                                                   gCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTa                                               <  2:564803/150‑1 (MQ=255)
                                                                                                     aGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTAcc                                             <  1:720827/150‑1 (MQ=255)
                                                                                                     aGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTAcc                                             >  1:2050421/1‑150 (MQ=255)
                                                                                                      gACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCa                                            <  1:237294/150‑1 (MQ=255)
                                                                                                       aCAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCAt                                           <  2:2779400/150‑1 (MQ=255)
                                                                                                       aCAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCAt                                           <  1:1717789/150‑1 (MQ=255)
                                                                                                        cAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATa                                          <  2:4492852/150‑1 (MQ=255)
                                                                                                            ccAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCt                                      >  2:2415980/1‑150 (MQ=255)
                                                                                                              aGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTg                                    >  2:3246109/1‑150 (MQ=255)
                                                                                                                 tGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCg                                 <  2:2092760/150‑1 (MQ=255)
                                                                                                                 tGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCg                                 >  1:2754707/1‑150 (MQ=255)
                                                                                                                  ggCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGTCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCgg                                >  1:386374/1‑150 (MQ=255)
                                                                                                                  ggCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCgg                                >  1:2810013/1‑150 (MQ=255)
                                                                                                                     gtgtGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGa                             >  1:4442414/1‑150 (MQ=255)
                                                                                                                       gtgCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATc                           >  1:4064555/1‑150 (MQ=255)
                                                                                                                           gTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTa                       >  1:577768/1‑150 (MQ=255)
                                                                                                                           gTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTa                       >  1:3359075/1‑150 (MQ=255)
                                                                                                                            tAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTaa                      >  1:1275246/1‑150 (MQ=255)
                                                                                                                                tGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGa                  <  2:3036508/150‑1 (MQ=255)
                                                                                                                                  ccAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATa                <  1:3468503/150‑1 (MQ=255)
                                                                                                                                  ccAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATa                >  2:4940677/1‑150 (MQ=255)
                                                                                                                                    aGGCATTGGTGATATACCGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATAcc              <  2:2832394/150‑1 (MQ=255)
                                                                                                                                    aGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATAcc              <  1:1962456/150‑1 (MQ=255)
                                                                                                                                       cATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGa           >  1:3951392/1‑150 (MQ=255)
                                                                                                                                       cATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGa           <  1:4785327/150‑1 (MQ=255)
                                                                                                                                         ttGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACa         <  2:1962388/150‑1 (MQ=255)
                                                                                                                                          tGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACAt        <  1:656401/150‑1 (MQ=255)
                                                                                                                                           ggTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATc       >  2:3344366/1‑150 (MQ=255)
                                                                                                                                              gATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTATGAGATACCAGACATccgc    <  2:3763902/150‑1 (MQ=255)
                                                                                                                                              gATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATccgc    >  2:2746556/1‑150 (MQ=255)
                                                                                                                                               atatAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATccgcc   >  1:834198/1‑150 (MQ=255)
                                                                                                                                                tataACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATccgccg  >  1:4680135/1‑150 (MQ=255)
                                                                                                                                                tataACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATccgccg  >  1:524908/1‑150 (MQ=255)
                                                                                                                                                tataACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATccgccg  >  2:5055674/1‑150 (MQ=255)
                                                                                                                                                tataACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATccgccg  >  2:2445815/1‑150 (MQ=255)
                                                                                                                                                tataACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATccgccg  >  1:3112201/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TATAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATCCGCCG  >  CP009273/1881498‑1881791

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: