Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,169,196 T→A 25.1% E57V (GAA→GTA)  yqiA ← acyl CoA esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,169,1960TA25.1% 578.5 / 215.6 307E57V (GAA→GTA) yqiAacyl CoA esterase
Reads supporting (aligned to +/- strand):  ref base T (96/133);  new base A (35/42);  total (132/175)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 310 total aligned reads displayed.

GTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGT  >  CP009273/3169048‑3169345
                                                                                                                                                    |                                                                                                                                                     
gTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGACTCACCGCCATGTTc                                                                                                                                                      >  1:1373797/1‑150 (MQ=255)
gTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCTCCATGTTc                                                                                                                                                      >  2:2121991/1‑150 (MQ=255)
  tCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTATTATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTcct                                                                                                                                                    <  2:4193077/150‑2 (MQ=255)
  tCTGACCGAGATAGTCCGTCAGCACTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCa                                                                                                                                                    <  1:3719566/150‑1 (MQ=255)
   cTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAg                                                                                                                                                   <  1:2913621/150‑1 (MQ=255)
    tGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAgg                                                                                                                                                  >  1:4971370/1‑150 (MQ=255)
    tGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTGTACCACCACTGCGGGCAGCATAAAACCTTGTGACAACCAGGTGGCGTACTATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAgg                                                                                                                                                  <  2:1758095/150‑1 (MQ=255)
      aCCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACCATACCCAGCGAATCACCGCCATGTTCCAGGAc                                                                                                                                                >  2:3040473/1‑150 (MQ=255)
        cGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACaa                                                                                                                                              >  1:4491261/1‑150 (MQ=255)
         gagaTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAAt                                                                                                                                             <  1:1057835/150‑1 (MQ=255)
         gagaTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAAt                                                                                                                                             >  2:1518148/1‑150 (MQ=255)
         gagaTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAAt                                                                                                                                             <  2:4503564/150‑1 (MQ=255)
          agaTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATg                                                                                                                                            <  2:3619129/150‑1 (MQ=255)
           gaTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATgg                                                                                                                                           <  2:3252008/150‑1 (MQ=255)
           gaTAGTCCGTCAGCAGGTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATgg                                                                                                                                           <  2:4676947/150‑1 (MQ=255)
             tAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGAt                                                                                                                                         >  1:574190/1‑150 (MQ=255)
             tAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGAt                                                                                                                                         <  2:5311958/150‑1 (MQ=255)
              aGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGAtt                                                                                                                                        <  1:2325079/150‑1 (MQ=255)
                 ccGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCca                                                                                                                                     <  2:3748087/150‑1 (MQ=255)
                 ccGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCca                                                                                                                                     >  2:1581754/1‑150 (MQ=255)
                   gTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCcagc                                                                                                                                   <  2:5104095/150‑1 (MQ=255)
                   gTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCcagc                                                                                                                                   <  1:349200/150‑1 (MQ=255)
                   gTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCcagc                                                                                                                                   >  2:3358418/1‑150 (MQ=255)
                        cagTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGctc                                                                                                                              <  2:1147037/150‑1 (MQ=255)
                              aaaCGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCcgc                                                                                                                        <  1:4467089/150‑1 (MQ=255)
                                aCGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGt                                                                                                                      <  1:2294605/150‑1 (MQ=255)
                                 cGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTc                                                                                                                     >  1:3095749/1‑150 (MQ=255)
                                 cGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTc                                                                                                                     <  2:3682684/150‑1 (MQ=255)
                                    gcgcACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGa                                                                                                                  <  2:4402543/150‑1 (MQ=255)
                                    gcgcACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACAGTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGa                                                                                                                  <  2:2474619/150‑1 (MQ=255)
                                     cgcACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGaa                                                                                                                 <  1:5160662/150‑1 (MQ=255)
                                       cACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAgg                                                                                                               <  2:1957526/150‑1 (MQ=255)
                                        aCCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGa                                                                                                              >  1:2941747/1‑150 (MQ=255)
                                         ccgccgGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGAt                                                                                                             <  2:1176917/150‑1 (MQ=255)
                                           gccgGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATAc                                                                                                           <  1:1790954/150‑1 (MQ=255)
                                            ccgGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAAGACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATAcg                                                                                                          <  2:3807516/150‑1 (MQ=255)
                                             cgGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATAcgg                                                                                                         >  1:2023598/1‑150 (MQ=255)
                                             cgGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATAcgg                                                                                                         >  2:1376992/1‑150 (MQ=255)
                                                gTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATAcggcgg                                                                                                      >  1:2398974/1‑150 (MQ=255)
                                                 tttACCACCACTGCGGGCAGCATAAAACAATGAGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAAAGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATAcggcgga                                                                                                     >  2:734111/1‑149 (MQ=255)
                                                     ccaccaCTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACt                                                                                                 <  1:1519697/150‑1 (MQ=255)
                                                      caccacTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTg                                                                                                >  1:623419/1‑150 (MQ=255)
                                                       accacTGCGGGCAGCATAAAACATTGTGACAACCAAGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACtac                                                                                               <  2:1902132/150‑3 (MQ=255)
                                                        ccacTTCGGGCAGCATAAAACATTGTGACAAACAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATAACCGCCATGTTCCAGGACAATTGATTCCAGCAGATCTGCCGCGTCTGAAGGATACGGCGGCAACTGCg                                                                                              <  1:3509940/150‑1 (MQ=255)
                                                        ccacTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCg                                                                                              <  2:86329/150‑1 (MQ=255)
                                                        ccacTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCg                                                                                              <  1:2646281/150‑1 (MQ=255)
                                                          acTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGa                                                                                            <  1:1167433/150‑1 (MQ=255)
                                                           cTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGaa                                                                                           >  1:883686/1‑150 (MQ=255)
                                                            tGCGGGCAGCATAAAACCTTGTGACATCCCGGTGGCGTAATATCCCCCCAGTGACGAACCGACAAGACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAAt                                                                                          <  2:1258388/150‑1 (MQ=255)
                                                            tGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAAt                                                                                          <  2:2351532/150‑1 (MQ=255)
                                                                  cagcaTAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGatcat                                                                                    >  1:5273294/1‑150 (MQ=255)
                                                                   agcaTAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGatcatc                                                                                   >  1:2750220/1‑150 (MQ=255)
                                                                   agcaTAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGatcatc                                                                                   <  2:3975639/150‑1 (MQ=255)
                                                                    gcaTAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCAtct                                                                                  >  1:933517/1‑150 (MQ=255)
                                                                    gcaTAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCAtct                                                                                  >  1:3206698/1‑150 (MQ=255)
                                                                      aTAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCa                                                                                >  2:2867777/1‑150 (MQ=255)
                                                                         aaaCATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTCCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACg                                                                             >  2:5029442/1‑150 (MQ=255)
                                                                         aaaCATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGAAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACg                                                                             <  1:2200710/150‑1 (MQ=255)
                                                                            cATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCa                                                                          >  2:2226908/1‑150 (MQ=255)
                                                                             aTTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAg                                                                         <  2:61404/150‑1 (MQ=255)
                                                                             aTTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAg                                                                         <  1:2742347/150‑1 (MQ=255)
                                                                                   aCAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGAt                                                                   <  1:1839471/150‑1 (MQ=255)
                                                                                    cAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATg                                                                  <  1:4433455/150‑1 (MQ=255)
                                                                                    cAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATg                                                                  <  1:2997078/150‑1 (MQ=255)
                                                                                     aaCCAGGTGGCGTAATATCCCCCCATTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGTATTCCAGCAGCTATGCCGCGTCGGAAGGATACGGCGGCAAATGCGGAATGATCATATCAACGACAGGGTGATGt                                                                 >  1:2436493/1‑150 (MQ=255)
                                                                                     aaCCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGt                                                                 >  1:5023902/1‑150 (MQ=255)
                                                                                        cAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTcc                                                              <  1:3636977/150‑1 (MQ=255)
                                                                                                gTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCgccagcca                                                      >  2:2757550/1‑150 (MQ=255)
                                                                                                gTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCgccagcca                                                      >  1:1691473/1‑150 (MQ=255)
                                                                                                gTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCgccagcca                                                      >  2:4073394/1‑150 (MQ=255)
                                                                                                 tAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGccagccag                                                     <  1:4020148/150‑1 (MQ=255)
                                                                                                   atatCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGtt                                                   <  2:4548465/150‑1 (MQ=255)
                                                                                                   atatCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGtt                                                   <  2:4205835/150‑1 (MQ=255)
                                                                                                   atatCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGtt                                                   <  1:3365762/150‑1 (MQ=255)
                                                                                                    tatCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGttt                                                  >  2:4292555/1‑150 (MQ=255)
                                                                                                       ccccccAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGtttttt                                               >  1:458646/1‑150 (MQ=255)
                                                                                                          cccAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAAc                                            <  2:723902/150‑1 (MQ=255)
                                                                                                          cccAGCGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGAGTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAAc                                            <  2:802438/150‑1 (MQ=255)
                                                                                                              gTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGc                                        >  1:1430152/1‑150 (MQ=255)
                                                                                                              gTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGc                                        >  1:4716098/1‑150 (MQ=255)
                                                                                                               tGACGAACCGACAATACCCAGCGAATCCCCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCt                                       <  2:2991210/150‑1 (MQ=255)
                                                                                                                gACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTc                                      >  1:2512389/1‑150 (MQ=255)
                                                                                                                gACGAACCGACAAGACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGGTCCGCCAGCCAGTTTTTTAACAAGCTc                                      <  2:1049788/150‑1 (MQ=255)
                                                                                                                 aCGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCg                                     >  1:1253209/1‑150 (MQ=255)
                                                                                                                   gAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCt                                   <  1:1916670/150‑1 (MQ=255)
                                                                                                                   gAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCt                                   <  1:4775891/150‑1 (MQ=255)
                                                                                                                   gAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCt                                   >  1:3065953/1‑150 (MQ=255)
                                                                                                                    aaCCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCtt                                  <  2:4971370/150‑1 (MQ=255)
                                                                                                                      ccGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCtttt                                <  2:2638510/150‑1 (MQ=255)
                                                                                                                      ccGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCtttt                                >  2:2551975/1‑150 (MQ=255)
                                                                                                                       cGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTg                               <  2:3496086/150‑1 (MQ=255)
                                                                                                                       cGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTg                               <  2:1454726/150‑1 (MQ=255)
                                                                                                                            agaCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCagag                          <  2:1285444/148‑1 (MQ=255)
                                                                                                                              aCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTGTTTTAACAAGCTCGCTTTTGCAGAgcg                        >  2:4808755/1‑150 (MQ=255)
                                                                                                                                 cAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCgcgg                     <  2:4761775/150‑1 (MQ=255)
                                                                                                                                        tCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTg              >  1:969090/1‑150 (MQ=255)
                                                                                                                                         cACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGt             <  2:1813902/150‑1 (MQ=255)
                                                                                                                                          aCCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGtt            >  2:586065/1‑150 (MQ=255)
                                                                                                                                            cGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGa          >  2:3223876/1‑150 (MQ=255)
                                                                                                                                               cATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAAc       >  2:4244026/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGT  >  CP009273/3169048‑3169345

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: