Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,726,289 A→T 9.1% P982P (CCA→CCT)  lhr → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,726,2890AT 9.1% 1162.3 / 86.4 353P982P (CCA→CCT) lhrputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base A (170/150);  new base T (18/14);  total (188/165)
Fisher's exact test for biased strand distribution p-value = 8.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

Only 100 of 364 total aligned reads displayed.

AGAAGAGCGAAAGGTGAAGAAGGATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCAGTTACCCGCTTCGCGACCTCTTTCTGCGCTATCT  >  CP009273/1726140‑1726438
                                                                                                                                                     |                                                                                                                                                     
agaagaGCGAAAGGTGAAGAAGGATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCa                                                                                                                                                       >  1:3060843/1‑150 (MQ=255)
agaagaGCGAAAGGTGAAGAAGGATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCa                                                                                                                                                       >  1:4815817/1‑150 (MQ=255)
agaagaGCGAAAGGTGAAGAAGGATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCa                                                                                                                                                       >  2:5134742/1‑150 (MQ=255)
 gaagaGCGAAAGGTGAAGAAGGATTGTTCGACCTGCTGCGCGAACTGGGGCCAATCACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCAGGAAAATCTTCTTGCAGTAAAACGAAACTTCCCAg                                                                                                                                                      >  2:5032493/1‑150 (MQ=255)
  aagaGCGAAAGGTGAAGAAGGATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGc                                                                                                                                                     >  1:4349621/1‑150 (MQ=255)
     agCGAAAGGTGAAGAAGGATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCTGCgat                                                                                                                                                  <  1:997184/150‑1 (MQ=255)
     agCGAAAGGTGAAGAAGGATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCgat                                                                                                                                                  <  2:1536757/150‑1 (MQ=255)
        gAAAGGTGAAGAAGGATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCgatgat                                                                                                                                               <  2:2846795/150‑1 (MQ=255)
         aaaGGTGAAGAAGGATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCTGCGATGAtt                                                                                                                                              >  1:4380781/1‑150 (MQ=255)
            ggTGAAGAAGGATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGc                                                                                                                                           >  2:2945332/1‑150 (MQ=255)
             gTGAAGAAGGATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAATCGAGTCTTCCCAGCGATGAGTAGCg                                                                                                                                          >  2:2681510/1‑150 (MQ=255)
                aagaagGATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGAc                                                                                                                                       <  1:3614406/150‑1 (MQ=255)
                 agaagGATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCTGCGATGATTAGCGGACa                                                                                                                                      >  1:2116616/1‑150 (MQ=255)
                   aagGATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCTGCGATGATTAGCGGACAgg                                                                                                                                    <  2:2462174/150‑1 (MQ=255)
                     ggATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAg                                                                                                                                  >  1:1656345/1‑150 (MQ=255)
                     ggATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAg                                                                                                                                  >  1:3132801/1‑150 (MQ=255)
                      gATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGc                                                                                                                                 >  1:4781556/1‑150 (MQ=255)
                      gATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGc                                                                                                                                 >  2:715673/1‑150 (MQ=255)
                        ttGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACCGGCAGCCGTGAAGAGGTTGCGGCGTATCTGGAAATTCTTCTTGCAGTAAAACGAATCTTCCCAGCGGTGAGTAGCGGACAGGAGCGt                                                                                                                               >  2:4487878/1‑150 (MQ=255)
                        ttGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGt                                                                                                                               <  2:5237746/150‑1 (MQ=255)
                         tGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTc                                                                                                                              <  2:5148737/150‑1 (MQ=255)
                          gTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCt                                                                                                                             >  2:2970561/1‑150 (MQ=255)
                          gTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCt                                                                                                                             >  1:2938120/1‑150 (MQ=255)
                           ttCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCtt                                                                                                                            <  2:1124093/150‑1 (MQ=255)
                           ttCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCtt                                                                                                                            <  2:4232249/150‑1 (MQ=255)
                            tCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTATAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTg                                                                                                                           <  1:1040950/150‑1 (MQ=255)
                            tCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTg                                                                                                                           <  1:3338646/150‑1 (MQ=255)
                             cGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGc                                                                                                                          >  1:3429798/1‑150 (MQ=255)
                             cGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGc                                                                                                                          <  1:1576582/150‑1 (MQ=255)
                                ccTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTg                                                                                                                       <  2:1758710/150‑1 (MQ=255)
                                  tgctgcGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTa                                                                                                                     >  1:3480718/1‑150 (MQ=255)
                                     tgcGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATgg                                                                                                                  >  1:1457430/1‑150 (MQ=255)
                                     tgcGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATgg                                                                                                                  <  1:5154400/150‑1 (MQ=255)
                                      gcgcGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGga                                                                                                                 <  1:5226682/150‑1 (MQ=255)
                                       cgcgAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGgat                                                                                                                <  2:1222151/150‑1 (MQ=255)
                                         cgAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGgatga                                                                                                              >  2:3670771/1‑150 (MQ=255)
                                          gAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGgatgat                                                                                                             >  2:4201948/1‑150 (MQ=255)
                                          gAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGgatgat                                                                                                             >  2:1948643/1‑150 (MQ=255)
                                           aaCTGGGGCCAATGACCGTTGTAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGatgatg                                                                                                            <  1:4507972/150‑1 (MQ=255)
                                             cTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATgcc                                                                                                          <  2:3100921/150‑1 (MQ=255)
                                             cTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATgcc                                                                                                          >  1:573534/1‑150 (MQ=255)
                                              tGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATgccg                                                                                                         <  1:4976632/150‑1 (MQ=255)
                                                gggCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATgccgcc                                                                                                       >  1:4317881/1‑150 (MQ=255)
                                                 ggCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCa                                                                                                      <  2:4583673/150‑1 (MQ=255)
                                                 ggCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCa                                                                                                      <  1:5127396/150‑1 (MQ=255)
                                                    cAATGACCGTTGAAGACCTGGCGCAACGTCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGc                                                                                                   <  2:3161826/150‑1 (MQ=255)
                                                     aaTGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCt                                                                                                  <  1:246198/150‑1 (MQ=255)
                                                      aTGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTg                                                                                                 >  1:2360701/1‑150 (MQ=255)
                                                       tGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGc                                                                                                <  1:2325942/150‑1 (MQ=255)
                                                            gTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGAt                                                                                           >  1:84831/1‑150 (MQ=255)
                                                              tGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGc                                                                                         <  2:3849372/150‑1 (MQ=255)
                                                               gAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGcc                                                                                        <  1:3309708/150‑1 (MQ=255)
                                                                  gACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTc                                                                                     <  2:3301430/150‑1 (MQ=255)
                                                                    ccTGGCGCAACGGCATACCGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCgg                                                                                   <  2:4018978/150‑1 (MQ=255)
                                                                      tGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCg                                                                                 >  1:2261967/1‑150 (MQ=255)
                                                                         cgcAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTAc                                                                              <  1:734043/150‑1 (MQ=255)
                                                                         cgcAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTAc                                                                              <  1:3369113/150‑1 (MQ=255)
                                                                               ggCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTAc                                                                        >  2:3312308/1‑150 (MQ=255)
                                                                                   tACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCaga                                                                    >  2:4467775/1‑150 (MQ=255)
                                                                                     cAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGt                                                                  <  2:1016160/150‑1 (MQ=255)
                                                                                     cAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGt                                                                  >  2:3456573/1‑150 (MQ=255)
                                                                                     cAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGt                                                                  <  1:4164368/150‑1 (MQ=255)
                                                                                       ggCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCa                                                                <  2:1411623/150‑1 (MQ=255)
                                                                                                 aaGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAgaga                                                      <  1:1672904/150‑1 (MQ=255)
                                                                                                  agagGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGat                                                     <  2:3772251/150‑1 (MQ=255)
                                                                                                    agGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCTGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGattt                                                   <  1:3758428/150‑1 (MQ=255)
                                                                                                    agGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGattt                                                   <  1:4037290/150‑1 (MQ=255)
                                                                                                     ggTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGattta                                                  <  1:1394026/150‑1 (MQ=255)
                                                                                                      gTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGatttat                                                 >  1:786859/1‑150 (MQ=255)
                                                                                                       gtGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGatttatt                                                <  1:3912399/149‑1 (MQ=255)
                                                                                                        tGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGatttattt                                               >  2:1796222/1‑150 (MQ=255)
                                                                                                             cgtATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATa                                          >  1:1060389/1‑150 (MQ=255)
                                                                                                              gtATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCATCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATag                                         <  2:1355417/150‑1 (MQ=255)
                                                                                                              gtATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATag                                         <  1:1815868/150‑1 (MQ=255)
                                                                                                              gtATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATag                                         >  1:3595920/1‑150 (MQ=255)
                                                                                                              gtATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATag                                         <  1:2672360/150‑1 (MQ=255)
                                                                                                                aTCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCATGCGATGATTAGAGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGCCGAACGACTACCAGAGTCATTGCCAGCGATTTATGTACATagag                                       >  2:2253467/1‑150 (MQ=255)
                                                                                                                 tCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCTGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGt                                      <  1:1991849/150‑1 (MQ=255)
                                                                                                                  cTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTc                                     <  2:1594457/150‑1 (MQ=255)
                                                                                                                  cTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTc                                     >  1:1208108/1‑150 (MQ=255)
                                                                                                                  cTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTc                                     >  1:212788/1‑150 (MQ=255)
                                                                                                                   tGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCa                                    >  2:2236671/1‑150 (MQ=255)
                                                                                                                   tGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCa                                    <  1:6683/150‑1 (MQ=255)
                                                                                                                    ggAAAATCTTCTTGCAGTAAAACGAATCTTCCCTGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCAg                                   <  2:4853656/150‑1 (MQ=255)
                                                                                                                       aaaTCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCAGTTa                                >  2:4849376/1‑150 (MQ=255)
                                                                                                                        aaTCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCAGTTAc                               <  1:389263/150‑1 (MQ=255)
                                                                                                                        aaTCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATGTACATAGAGTCAGTTAc                               >  1:2638662/1‑150 (MQ=255)
                                                                                                                         aTCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCAGTTAcc                              >  1:3644984/1‑150 (MQ=255)
                                                                                                                          tcttctTGCAGTAAAACGAATCTTCCCTGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCAGTTAccc                             <  1:4599127/150‑1 (MQ=255)
                                                                                                                           cttcttGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCAGTTACCCg                            >  1:38209/1‑150 (MQ=255)
                                                                                                                             tcttGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCAGTTACCCGCt                          >  2:654295/1‑150 (MQ=255)
                                                                                                                              cttGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCAGTTACCCGCtt                         <  1:2036326/150‑1 (MQ=255)
                                                                                                                              cttGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCAGTTACCCGCtt                         >  1:1768589/1‑150 (MQ=255)
                                                                                                                               ttGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCAGTTACCCGCTTc                        <  1:4584166/150‑1 (MQ=255)
                                                                                                                                tGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCAGTTACCCGCTTcg                       <  2:2562503/150‑1 (MQ=255)
                                                                                                                                 gCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCAGTTACCCGCTTcgc                      >  2:3386012/1‑150 (MQ=255)
                                                                                                                                          cGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCAGTTACCCGCTTCGCGACCTCttt             >  2:4698945/1‑150 (MQ=255)
                                                                                                                                               cTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCAGTTACCCGCTTCGCGACCTCTTTCTgcg        >  1:1872035/1‑150 (MQ=255)
                                                                                                                                                 tCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCAGTTACCCGCTTCGCGACCTCTTTCTGCGCt      >  2:4547171/1‑150 (MQ=255)
                                                                                                                                                   ccAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCAGTTACCCGCTTCGCGACCTCTTTCTGCGCTAt    <  2:829135/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGAAGAGCGAAAGGTGAAGAAGGATTGTTCGACCTGCTGCGCGAACTGGGGCCAATGACCGTTGAAGACCTGGCGCAACGGCATACAGGCAGCAGTGAAGAGGTTGCGTCGTATCTGGAAAATCTTCTTGCAGTAAAACGAATCTTCCCAGCGATGATTAGCGGACAGGAGCGTCTTGCCTGTATGGATGATGCCGCCAGGCTGCGTGATGCCCTCGGCGTACGACTACCAGAGTCATTGCCAGAGATTTATTTACATAGAGTCAGTTACCCGCTTCGCGACCTCTTTCTGCGCTATCT  >  CP009273/1726140‑1726438

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