Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,470,873 G→T 64.9% R928S (CGT→AGT)  valS ← valyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,470,8730GT64.9% 164.5 / 262.3 217R928S (CGT→AGT) valSvalyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (40/36);  new base T (60/80);  total (101/116)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

CTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGGCCAGCCATCG  >  CP009273/4470724‑4471022
                                                                                                                                                     |                                                                                                                                                     
cttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACg                                                                                                                                                       >  2:1680497/1‑150 (MQ=255)
cttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACg                                                                                                                                                       <  1:1512092/150‑1 (MQ=255)
    tttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCActctct                                                                                                                                                   <  2:3594675/150‑1 (MQ=255)
     ttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCAGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGACTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCATCCAGCTTCTCACGCTCtt                                                                                                                                                  >  2:773600/1‑150 (MQ=255)
       aCCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTc                                                                                                                                                <  1:4115685/150‑1 (MQ=255)
          aTTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTcgcg                                                                                                                                             <  1:2889915/150‑1 (MQ=255)
            ttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGAt                                                                                                                                           <  2:2060463/150‑1 (MQ=255)
                gTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGAcc                                                                                                                                       <  1:1046385/150‑1 (MQ=255)
                 tAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCg                                                                                                                                      <  2:3362332/150‑1 (MQ=255)
                 gaaGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCg                                                                                                                                      <  2:2291292/149‑1 (MQ=255)
                  aaGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGc                                                                                                                                     >  1:2966254/1‑150 (MQ=255)
                    gagaTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCtt                                                                                                                                   <  1:610229/150‑1 (MQ=255)
                     agaTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTc                                                                                                                                  <  1:2646072/150‑1 (MQ=255)
                     agaTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTc                                                                                                                                  >  1:871251/1‑150 (MQ=255)
                        tAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCgg                                                                                                                               <  1:4302985/150‑1 (MQ=255)
                         aaaaatGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGt                                                                                                                              >  1:4264657/1‑150 (MQ=255)
                         aaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGt                                                                                                                              >  1:2421226/1‑150 (MQ=255)
                          aaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTTACGCTCTTTCGCGATGACCGCTTCCGGTg                                                                                                                             >  1:3308739/1‑150 (MQ=255)
                          aaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTg                                                                                                                             >  2:170928/1‑150 (MQ=255)
                          aaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTg                                                                                                                             <  1:2709426/150‑1 (MQ=255)
                          aaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGAGGACCGCTTCCGGTg                                                                                                                             <  2:780265/150‑1 (MQ=255)
                                  cGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCg                                                                                                                     <  2:664356/150‑1 (MQ=255)
                                   ggAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGa                                                                                                                    <  2:382811/150‑1 (MQ=255)
                                   ggAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGa                                                                                                                    <  2:382815/150‑1 (MQ=255)
                                   ggAGCATGCTCCGGCCTTCGTTTACATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGa                                                                                                                    >  2:3938574/1‑150 (MQ=255)
                                   ggAGCATGCTCCGGCCTTCGTTATCATCACTGAGTTTTGATTACAGCGCGGAGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCATTACAGCTTCTAACTCTCTTTCGCGATGAACGATTACGGTGCGAGGGCGa                                                                                                                    >  1:4332765/1‑150 (MQ=255)
                                    gAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGAc                                                                                                                   >  2:532380/1‑150 (MQ=255)
                                    gAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGAc                                                                                                                   <  1:2027899/150‑1 (MQ=255)
                                    gAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGAc                                                                                                                   <  2:3474295/150‑1 (MQ=255)
                                         tGCTCCGGCCTTCTTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCAGCATAGCCTTCCAGCTTCTAACTCTCTTTCGCGATGACCGCTTCCGGTGAGCGGGCGACaaatc                                                                                                              >  2:4236698/1‑148 (MQ=255)
                                           cTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCt                                                                                                            <  2:2132961/150‑1 (MQ=255)
                                            tCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCtt                                                                                                           <  1:4047209/150‑1 (MQ=255)
                                             ccGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGAGCGCGGGCGACAAATCCttc                                                                                                          >  2:1644544/1‑150 (MQ=255)
                                                gCCTTCGTTGTCATCACTGTGTTTTGAGTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTGTCGCTGTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCttcgtt                                                                                                       <  2:4239138/150‑1 (MQ=255)
                                                 ccTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCttcgttc                                                                                                      <  1:2259842/150‑1 (MQ=255)
                                                        tttCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGtt                                                                                               <  2:2895788/150‑1 (MQ=255)
                                                         ttCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGttt                                                                                              >  1:28682/1‑150 (MQ=255)
                                                         ttCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGttt                                                                                              <  2:1685505/150‑1 (MQ=255)
                                                         ttCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTACCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTACGTTAGCCAGttt                                                                                              >  1:2724061/1‑150 (MQ=255)
                                                          tcatcaCTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTg                                                                                             <  2:47325/150‑1 (MQ=255)
                                                           catcaCTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGt                                                                                            <  2:3960970/150‑1 (MQ=255)
                                                            atcaCTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGtt                                                                                           >  2:1240381/1‑150 (MQ=255)
                                                            atcaCTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGtt                                                                                           <  2:4071252/150‑1 (MQ=255)
                                                            atcaCTGTGTTGTGATTACATCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGtt                                                                                           <  2:4303500/150‑1 (MQ=255)
                                                              caCTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTtct                                                                                         >  1:3282773/1‑150 (MQ=255)
                                                               aCTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTtctc                                                                                        >  1:580929/1‑150 (MQ=255)
                                                               aCTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTtctc                                                                                        >  1:3527148/1‑150 (MQ=255)
                                                               aCTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTtctc                                                                                        >  1:3526352/1‑150 (MQ=255)
                                                                cTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCg                                                                                       <  1:3863524/150‑1 (MQ=255)
                                                                cTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCg                                                                                       <  1:191195/150‑1 (MQ=255)
                                                                cTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCg                                                                                       <  1:2138478/150‑1 (MQ=255)
                                                                  gtgtTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGAt                                                                                     <  2:533853/150‑1 (MQ=255)
                                                                    gtTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATAc                                                                                   >  2:4007960/1‑150 (MQ=255)
                                                                    gtTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATAc                                                                                   <  2:3190236/150‑1 (MQ=255)
                                                                        tGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCt                                                                               >  1:2116793/1‑150 (MQ=255)
                                                                        tGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCt                                                                               >  2:2303296/1‑150 (MQ=255)
                                                                          aTTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGa                                                                             >  1:2594913/1‑150 (MQ=255)
                                                                            tACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGAtt                                                                           >  2:2186113/1‑150 (MQ=255)
                                                                                      cgATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCTGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCaa                                                                 >  1:4429261/1‑150 (MQ=255)
                                                                                       gATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAAt                                                                <  1:1001558/150‑1 (MQ=255)
                                                                                         tAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCt                                                              >  2:1132213/1‑150 (MQ=255)
                                                                                           aCAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTc                                                            <  1:1154601/150‑1 (MQ=255)
                                                                                            cAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCg                                                           >  2:2673835/1‑150 (MQ=255)
                                                                                              gCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGcc                                                         >  1:3551418/1‑150 (MQ=255)
                                                                                              gCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACAGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGcc                                                         >  1:2404882/1‑150 (MQ=255)
                                                                                                 tgctgTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCAct                                                      <  1:1379637/150‑1 (MQ=255)
                                                                                                  gctgTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCActt                                                     <  2:1473548/150‑1 (MQ=255)
                                                                                                   ctgTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCActtc                                                    >  2:3995241/1‑150 (MQ=255)
                                                                                                   ctgTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCActtc                                                    <  2:871118/150‑1 (MQ=255)
                                                                                                     gTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCActtctt                                                  <  1:170928/150‑1 (MQ=255)
                                                                                                      ttCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCttt                                                 <  1:1638581/150‑1 (MQ=255)
                                                                                                      ttCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCttt                                                 <  1:3900329/150‑1 (MQ=255)
                                                                                                       tCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTc                                                <  2:3736720/150‑1 (MQ=255)
                                                                                                       tCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTc                                                >  1:1011647/1‑150 (MQ=255)
                                                                                                        cAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCg                                               >  2:2011749/1‑150 (MQ=255)
                                                                                                         aaTCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGc                                              >  1:2379171/1‑150 (MQ=255)
                                                                                                         aaTCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGc                                              >  1:940163/1‑150 (MQ=255)
                                                                                                          aTCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGcc                                             <  1:3676886/150‑1 (MQ=255)
                                                                                                           tCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCa                                            >  1:1126470/1‑150 (MQ=255)
                                                                                                            cAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAg                                           >  2:2229531/1‑150 (MQ=255)
                                                                                                            cAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGATTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATAAGGCTGATATCACCTTCAATCTACGCCACTTCTATCGCCAg                                           >  2:2229601/1‑150 (MQ=255)
                                                                                                             aGTTTCGCTTTCGCTTCCGCATAGCCTTACAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGa                                          <  1:4389419/150‑1 (MQ=255)
                                                                                                                ttcgctttcgctTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGAcgc                                       >  2:25230/1‑150 (MQ=255)
                                                                                                                ttcgctttcgctTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGAcgc                                       >  1:740536/1‑150 (MQ=255)
                                                                                                                 tcgctttcgcttCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGAcgcg                                      >  1:4420293/1‑150 (MQ=255)
                                                                                                                  cgctttcgcttCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACgcgc                                     >  1:3749557/1‑150 (MQ=255)
                                                                                                                   gctttcgcttCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGcc                                    <  2:2116793/150‑1 (MQ=255)
                                                                                                                     tttcgcttCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAg                                  >  2:1043471/1‑150 (MQ=255)
                                                                                                                         gcttCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATCCGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCa                              >  2:3009097/1‑150 (MQ=255)
                                                                                                                         gcttCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCa                              >  2:3008642/1‑150 (MQ=255)
                                                                                                                           ttCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCAtc                            >  1:359663/1‑150 (MQ=255)
                                                                                                                           ttCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCAtc                            >  2:4087410/1‑150 (MQ=255)
                                                                                                                             ccGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCAtctt                          <  1:3994036/150‑1 (MQ=255)
                                                                                                                              cGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCAtcttc                         <  2:3975701/150‑1 (MQ=255)
                                                                                                                                cATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATcttctt                       >  1:4076602/1‑150 (MQ=255)
                                                                                                                                       ttCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTtgat                >  2:2194092/1‑150 (MQ=255)
                                                                                                                                             cTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGgcc          <  2:4278543/150‑1 (MQ=255)
                                                                                                                                                   aCGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGGCCAGCCAt    <  2:3728182/150‑1 (MQ=255)
                                                                                                                                                   aCGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGGCCAGCCAt    >  1:217337/1‑150 (MQ=255)
                                                                                                                                                     gCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGGCCAGCCATCg  >  1:3044712/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGGCCAGCCATCG  >  CP009273/4470724‑4471022

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