Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 294,929 C→A 7.6% E573* (GAA→TAA)  paoC ← PaoABC aldehyde oxidoreductase, Moco‑containing subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273294,9290CA 7.6% 1076.4 / 39.5 284E573* (GAA→TAA) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base C (142/120);  new base A (8/14);  total (150/134)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 287 total aligned reads displayed.

CGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATAC  >  CP009273/294780‑295078
                                                                                                                                                     |                                                                                                                                                     
 gcAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCa                                                                                                                                                      <  1:902276/150‑1 (MQ=255)
     ttttCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGc                                                                                                                                                  >  1:4151226/1‑150 (MQ=255)
        tCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGctct                                                                                                                                               <  2:3650289/150‑1 (MQ=255)
         cGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGctctc                                                                                                                                              >  2:1515919/1‑150 (MQ=255)
          ggATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGctctct                                                                                                                                             >  1:3172539/1‑150 (MQ=255)
           gATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCtt                                                                                                                                            >  1:2307471/1‑150 (MQ=255)
           gATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCtt                                                                                                                                            >  2:3020442/1‑150 (MQ=255)
             ttCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTcc                                                                                                                                          <  1:3023846/150‑1 (MQ=255)
                aGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCt                                                                                                                                       <  2:3484185/150‑1 (MQ=255)
                aGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACATCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCt                                                                                                                                       >  1:1787377/1‑150 (MQ=255)
                 ggATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTg                                                                                                                                      >  1:3528209/1‑150 (MQ=255)
                           tGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTgcc                                                                                                                            >  1:605391/1‑150 (MQ=255)
                           tGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTgcc                                                                                                                            <  1:2374526/150‑1 (MQ=255)
                              agcaCACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTgccgcc                                                                                                                         >  1:739273/1‑150 (MQ=255)
                               gcaCACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCt                                                                                                                        <  1:2036000/150‑1 (MQ=255)
                                 acacacaGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGc                                                                                                                      <  1:4591469/150‑1 (MQ=255)
                                  cacacaGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCg                                                                                                                     >  1:1726603/1‑150 (MQ=255)
                                     acaGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTg                                                                                                                  >  1:2238654/1‑150 (MQ=255)
                                     acaGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTg                                                                                                                  >  2:2413917/1‑150 (MQ=255)
                                          gAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTc                                                                                                             >  1:2847437/1‑150 (MQ=255)
                                            gCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCAt                                                                                                           <  2:3216168/150‑1 (MQ=255)
                                              aTACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGt                                                                                                         <  1:188852/150‑1 (MQ=255)
                                              aTACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGt                                                                                                         >  1:3975875/1‑150 (MQ=255)
                                              aTACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGt                                                                                                         >  1:4429405/1‑150 (MQ=255)
                                                   ccGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGt                                                                                                    >  2:4080966/1‑150 (MQ=255)
                                                   ccGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGt                                                                                                    >  2:1673532/1‑150 (MQ=255)
                                                    cGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTg                                                                                                   >  1:67431/1‑150 (MQ=255)
                                                     ggACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTgg                                                                                                  >  2:4254395/1‑150 (MQ=255)
                                                       aacagAAATTCTCCCGTCGCTCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGgcg                                                                                                <  1:379385/146‑1 (MQ=255)
                                                         ccGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCt                                                                                              >  2:3801291/1‑150 (MQ=255)
                                                           gAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTc                                                                                            <  2:399471/150‑1 (MQ=255)
                                                               ttCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCACAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGt                                                                                        <  2:14668/150‑1 (MQ=255)
                                                                 ctcCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTAc                                                                                      >  1:3770323/1‑150 (MQ=255)
                                                                  tcCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTAcc                                                                                     <  1:4350157/150‑1 (MQ=255)
                                                                   cccGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCg                                                                                    <  1:1775994/150‑1 (MQ=255)
                                                                    ccGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGt                                                                                   >  1:3357871/1‑150 (MQ=255)
                                                                    ccGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGt                                                                                   >  1:4005970/1‑150 (MQ=255)
                                                                       tCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTgg                                                                                <  1:4522779/150‑1 (MQ=255)
                                                                        cgcgCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGt                                                                               <  1:343871/150‑1 (MQ=255)
                                                                        cgcgCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGt                                                                               >  1:1480321/1‑150 (MQ=255)
                                                                           gcTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAAt                                                                            >  2:3147327/1‑150 (MQ=255)
                                                                             tATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCt                                                                          <  1:4324173/150‑1 (MQ=255)
                                                                              aTGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCtt                                                                         <  2:4371443/150‑1 (MQ=255)
                                                                              aTGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCtt                                                                         <  2:2627919/150‑1 (MQ=255)
                                                                               tGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTg                                                                        >  1:3582749/1‑150 (MQ=255)
                                                                               tGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGGTGGTAATCTTg                                                                        >  2:3327908/1‑150 (MQ=255)
                                                                                  aCGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCg                                                                     >  1:467428/1‑150 (MQ=255)
                                                                                     ccGACCTCCACAAACTGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCg                                                                  <  2:1328673/150‑1 (MQ=255)
                                                                                     ccGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCg                                                                  <  1:2943655/150‑1 (MQ=255)
                                                                                       gACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGc                                                                >  1:1271708/1‑150 (MQ=255)
                                                                                       gACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGc                                                                >  2:2878199/1‑150 (MQ=255)
                                                                                       gACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGc                                                                <  2:1826230/150‑1 (MQ=255)
                                                                                         ccTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCaa                                                              <  1:3641619/150‑1 (MQ=255)
                                                                                         ccTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCaa                                                              <  1:2681492/150‑1 (MQ=255)
                                                                                            ccACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACt                                                           <  2:3431648/150‑1 (MQ=255)
                                                                                                aaaaTGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGa                                                       >  2:4536512/1‑150 (MQ=255)
                                                                                                   aTGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTg                                                    >  1:1013720/1‑150 (MQ=255)
                                                                                                    tGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGc                                                   >  2:1717372/1‑150 (MQ=255)
                                                                                                      cccGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTc                                                 <  1:2299573/150‑1 (MQ=255)
                                                                                                       ccGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCa                                                <  2:4440969/150‑1 (MQ=255)
                                                                                                        cGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAg                                               >  2:4211712/1‑150 (MQ=255)
                                                                                                         ggCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAgg                                              >  2:483248/1‑150 (MQ=255)
                                                                                                         ggCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAgg                                              <  2:830218/150‑1 (MQ=255)
                                                                                                         ggCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAgg                                              >  1:3264233/1‑150 (MQ=255)
                                                                                                          gCAATGGTCGACTGCTGGTACTCTTTGCTAATTGTTCCGAATTCAATGCTCTCTTCCGCTGTCATTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTTCTCAGGa                                             <  1:1336389/150‑1 (MQ=255)
                                                                                                          gCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGa                                             <  1:2364947/150‑1 (MQ=255)
                                                                                                          gCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGATCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGa                                             >  1:52051/1‑150 (MQ=255)
                                                                                                           cAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGAt                                            >  2:780313/1‑150 (MQ=255)
                                                                                                            aaaGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATc                                           >  2:2266482/1‑150 (MQ=255)
                                                                                                            aaaGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATc                                           >  1:3013320/1‑150 (MQ=255)
                                                                                                            aaaGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATc                                           >  2:563956/1‑150 (MQ=255)
                                                                                                            aaaGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATc                                           >  1:2555884/1‑150 (MQ=255)
                                                                                                            aaaGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATc                                           >  2:3066878/1‑150 (MQ=255)
                                                                                                                   gACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCg                                    >  1:2048947/1‑150 (MQ=255)
                                                                                                                         tgGTACTCTTTGCTCAGTGGTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGcc                              <  1:4234432/150‑1 (MQ=255)
                                                                                                                          ggTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCg                             <  1:821257/150‑1 (MQ=255)
                                                                                                                          ggTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCg                             >  1:2782703/1‑150 (MQ=255)
                                                                                                                             aCTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAgg                          <  1:2452303/150‑1 (MQ=255)
                                                                                                                              ctctTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGc                         >  1:2669330/1‑150 (MQ=255)
                                                                                                                               tctTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCa                        >  2:539232/1‑150 (MQ=255)
                                                                                                                                ctTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCaa                       >  1:1032094/1‑150 (MQ=255)
                                                                                                                                 tttGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAAt                      >  2:4349829/1‑150 (MQ=255)
                                                                                                                                 tttGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAAt                      >  2:2695051/1‑150 (MQ=255)
                                                                                                                                  ttGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATc                     <  2:1868446/150‑1 (MQ=255)
                                                                                                                                     cTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCAtt                  >  1:4495645/1‑150 (MQ=255)
                                                                                                                                      tCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCAttt                 <  2:1743670/150‑1 (MQ=255)
                                                                                                                                        aGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTcg               <  2:604651/150‑1 (MQ=255)
                                                                                                                                          tgtTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTcgcg             <  1:1177034/150‑1 (MQ=255)
                                                                                                                                          tgtTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTcgcg             <  2:3602828/150‑1 (MQ=255)
                                                                                                                                            ttCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCAGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGaa           >  2:1181017/1‑150 (MQ=255)
                                                                                                                                             tCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAg          <  1:3733356/150‑1 (MQ=255)
                                                                                                                                              ccGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGc         >  1:3685374/1‑150 (MQ=255)
                                                                                                                                               cGAATTCAATGCTCTCTTCCGCTGTCAGTCTGTCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCt        <  2:4491566/150‑1 (MQ=255)
                                                                                                                                                gAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCtt       <  2:2546242/150‑1 (MQ=255)
                                                                                                                                                 aaTTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTc      >  2:2182923/1‑150 (MQ=255)
                                                                                                                                                 aaTTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTc      >  2:3976859/1‑150 (MQ=255)
                                                                                                                                                  aTTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCa     >  1:936393/1‑150 (MQ=255)
                                                                                                                                                   ttCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCAt    <  1:2399062/150‑1 (MQ=255)
                                                                                                                                                    tCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATa   <  2:2542320/150‑1 (MQ=255)
                                                                                                                                                     cAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATAc  >  1:3483685/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATAC  >  CP009273/294780‑295078

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: