Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,882,760 T→A 6.4% V224E (GTG→GAG)  tnaA → tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,882,7600TA 6.4% 1076.6 / 46.9 311V224E (GTG→GAG) tnaAtryptophanase/L‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (142/149);  new base A (7/13);  total (149/162)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 326 total aligned reads displayed.

GGGATTAGAACGCGGTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGATGCGATGG  >  CP009273/3882611‑3882909
                                                                                                                                                     |                                                                                                                                                     
 ggATTAGAACGCGGTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCggtg                                                                                                                                                      <  1:4084292/150‑1 (MQ=255)
 ggATTAGAACGCGGTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCggtg                                                                                                                                                      <  1:2936345/150‑1 (MQ=255)
   aTTAGAACGCGGTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCggtggt                                                                                                                                                    >  1:226206/1‑150 (MQ=255)
     tAGAACGCGGTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTaa                                                                                                                                                  >  1:3724047/1‑150 (MQ=255)
     tAGAACGCGGTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTaa                                                                                                                                                  >  1:1789257/1‑150 (MQ=255)
           gcgGAATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACt                                                                                                                                            <  1:207408/150‑1 (MQ=255)
            cgGTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTc                                                                                                                                           >  1:3518357/1‑150 (MQ=255)
             ggTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTcc                                                                                                                                          >  1:3143483/1‑150 (MQ=255)
             ggTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTcc                                                                                                                                          <  2:3750869/150‑1 (MQ=255)
               tATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCcgc                                                                                                                                        >  1:2439526/1‑150 (MQ=255)
                 ttGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCcgcgc                                                                                                                                      >  2:602928/1‑150 (MQ=255)
                  tGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCcgcgcg                                                                                                                                     >  1:1134375/1‑150 (MQ=255)
                   gaagaaGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCgcgcgc                                                                                                                                    >  1:3081046/1‑150 (MQ=255)
                     agaagTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCtt                                                                                                                                  <  1:2275764/150‑1 (MQ=255)
                     agaagTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGAGGTAATGGACTCCGCGCGCtt                                                                                                                                  <  2:4637848/150‑1 (MQ=255)
                      gaagTTGGTCCGATTAACGTGCAGTATATCGTAGCAACGATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAATCTTAAATGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCTCttt                                                                                                                                 <  2:4980198/150‑1 (MQ=255)
                       aagTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTg                                                                                                                                >  1:2196993/1‑150 (MQ=255)
                         gTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCt                                                                                                                              >  1:51639/1‑150 (MQ=255)
                           tGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGa                                                                                                                            >  2:5155076/1‑150 (MQ=255)
                           tGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGa                                                                                                                            <  1:58092/150‑1 (MQ=255)
                            ggTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGaa                                                                                                                           >  1:4511433/1‑150 (MQ=255)
                               ccGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAAc                                                                                                                        <  2:1815360/150‑1 (MQ=255)
                                  aaTAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGcc                                                                                                                     >  1:3278504/1‑150 (MQ=255)
                                    tAACGTGCCGTATATCGTTGCAACCATGACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTa                                                                                                                   <  2:4443811/150‑1 (MQ=255)
                                    tAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTa                                                                                                                   >  1:5123120/1‑150 (MQ=255)
                                     aaCGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTAt                                                                                                                  >  1:1816308/1‑150 (MQ=255)
                                     aaCGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATAGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTAt                                                                                                                  <  1:862184/150‑1 (MQ=255)
                                      aCGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTAtt                                                                                                                 <  1:389378/150‑1 (MQ=255)
                                        gTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTtc                                                                                                               >  1:1986643/1‑150 (MQ=255)
                                          gCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTtcat                                                                                                             <  1:2606668/150‑1 (MQ=255)
                                          gCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTtcat                                                                                                             <  2:4967148/150‑1 (MQ=255)
                                             gTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCaa                                                                                                          >  2:2043136/1‑150 (MQ=255)
                                               atatCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAagc                                                                                                        <  1:4357767/150‑1 (MQ=255)
                                                  tCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAagcagc                                                                                                     <  1:1776914/150‑1 (MQ=255)
                                                    gTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGt                                                                                                   >  1:4775871/1‑150 (MQ=255)
                                                         aaCCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGc                                                                                              >  2:4266342/1‑150 (MQ=255)
                                                         aaCCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGc                                                                                              >  1:4043567/1‑150 (MQ=255)
                                                            cATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGa                                                                                           <  2:2546829/150‑1 (MQ=255)
                                                            cATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGa                                                                                           <  2:4120594/150‑1 (MQ=255)
                                                            cATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTGTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGa                                                                                           <  2:1512223/150‑1 (MQ=255)
                                                             aTCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGaa                                                                                          >  2:4191207/1‑150 (MQ=255)
                                                             aTCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGaa                                                                                          <  2:3966650/150‑1 (MQ=255)
                                                              tCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTTGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAAt                                                                                         <  1:1540424/150‑1 (MQ=255)
                                                              tCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAAt                                                                                         <  1:4582989/150‑1 (MQ=255)
                                                               cACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATa                                                                                        >  1:2649716/1‑150 (MQ=255)
                                                                 ccAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACa                                                                                      >  1:3856200/1‑150 (MQ=255)
                                                                 ccAGTAACTCTGCAGGTGGTCAGCCGGGTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTGTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACa                                                                                      <  2:3343469/150‑1 (MQ=255)
                                                                    gTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAg                                                                                   >  1:976134/1‑150 (MQ=255)
                                                                       aCTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACt                                                                                <  2:1847554/150‑1 (MQ=255)
                                                                        ctctGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTg                                                                               >  1:491015/1‑150 (MQ=255)
                                                                         tctGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTgg                                                                              >  2:114743/1‑150 (MQ=255)
                                                                         tctGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTgg                                                                              <  1:5020556/150‑1 (MQ=255)
                                                                           tGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGAc                                                                            <  1:3305083/150‑1 (MQ=255)
                                                                           tGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGACTATATCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGAc                                                                            >  2:1654100/1‑150 (MQ=255)
                                                                             cAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCa                                                                          <  2:3303228/150‑1 (MQ=255)
                                                                             cAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTGTGCTGAAAACGCCTATTTCATCCAGCAGCGTGAAGCAGAATACAAAGACTGGAccc                                                                          >  2:4517393/1‑149 (MQ=255)
                                                                                gtggtCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCg                                                                       <  1:597288/150‑1 (MQ=255)
                                                                                 tggtCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGa                                                                      <  2:2483285/150‑1 (MQ=255)
                                                                                  ggtCAGCCGTTTTCACTGGCAAACTTCAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTGTGCGGAAAACGCCTATTTCATCGAGCAGCGTGAAGCTGTATACTAAGACTGGACCACCGAg                                                                     <  2:3649139/150‑1 (MQ=255)
                                                                                  ggtCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAg                                                                     <  1:544994/150‑1 (MQ=255)
                                                                                   gtCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGc                                                                    >  1:3955778/1‑150 (MQ=255)
                                                                                      aGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGa                                                                 >  2:4346582/1‑150 (MQ=255)
                                                                                       gCCGGTTTCACTTGCAAACTTAACAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGAt                                                                <  2:3439853/150‑1 (MQ=255)
                                                                                          ggTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCAc                                                             <  1:1223172/150‑1 (MQ=255)
                                                                                           gTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGACTACAAAGACTGGACACTCGAGCAGATCAcc                                                            >  2:2817007/1‑150 (MQ=255)
                                                                                           gTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCAcc                                                            <  1:761320/150‑1 (MQ=255)
                                                                                           gTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCAcc                                                            <  2:4901982/150‑1 (MQ=255)
                                                                                              tCACTGGCAAACTTAAAAGAGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCcgc                                                         <  2:2765849/150‑1 (MQ=255)
                                                                                                 cTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGaa                                                      <  1:3634949/150‑1 (MQ=255)
                                                                                                 cTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGaa                                                      <  1:1014623/150‑1 (MQ=255)
                                                                                                       aaCTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTAc                                                >  2:4398545/1‑150 (MQ=255)
                                                                                                           tAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAat                                            <  2:2932320/150‑1 (MQ=255)
                                                                                                             aaaGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAatat                                          >  2:3182827/1‑150 (MQ=255)
                                                                                                              aaGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATg                                         >  2:5090369/1‑150 (MQ=255)
                                                                                                              aaGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATg                                         <  1:2349945/150‑1 (MQ=255)
                                                                                                                  gATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGa                                     >  2:898863/1‑150 (MQ=255)
                                                                                                                  gATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGa                                     <  1:3914662/150‑1 (MQ=255)
                                                                                                                   aTGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGat                                    >  2:2310673/1‑150 (MQ=255)
                                                                                                                     gTACAGCATCGCGAAGAAATACGATATTCCGGTGGTCATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGatat                                  >  1:2475520/1‑150 (MQ=255)
                                                                                                                     gTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGatat                                  <  2:2723945/150‑1 (MQ=255)
                                                                                                                      tACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATg                                 >  1:1710132/1‑150 (MQ=255)
                                                                                                                      tACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATg                                 >  1:2492461/1‑150 (MQ=255)
                                                                                                                      tACAGCATCGCGAAGAAATACGATATTCCGGAGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATg                                 >  1:4149735/1‑150 (MQ=255)
                                                                                                                        cAGCATCGCGAAGAAATACGATATTCCGGAGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCt                               <  1:2882878/150‑1 (MQ=255)
                                                                                                                         aGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTg                              <  2:2410819/150‑1 (MQ=255)
                                                                                                                         aGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTg                              <  1:1881149/150‑1 (MQ=255)
                                                                                                                           cATCGCGAAGAAATAAGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGc                            <  2:1692262/150‑1 (MQ=255)
                                                                                                                            aTCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCg                           >  1:3648282/1‑150 (MQ=255)
                                                                                                                             tCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGa                          <  1:5200407/150‑1 (MQ=255)
                                                                                                                              cgcgAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGAt                         >  2:2880755/1‑150 (MQ=255)
                                                                                                                                cgAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGt                       <  2:3570016/150‑1 (MQ=255)
                                                                                                                                 gaagaaATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTc                      >  2:2646126/1‑150 (MQ=255)
                                                                                                                                 gaagaaATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGAAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTc                      >  2:525626/1‑150 (MQ=255)
                                                                                                                                     aaaTACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGcc                  >  2:1363085/1‑150 (MQ=255)
                                                                                                                                        taagATATACCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGAAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCaag               >  2:2524788/4‑150 (MQ=255)
                                                                                                                                        tACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCaag               <  2:180846/150‑1 (MQ=255)
                                                                                                                                             tatTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCaagaaaga          >  1:3410171/1‑150 (MQ=255)
                                                                                                                                              atTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGAt         >  2:825310/1‑150 (MQ=255)
                                                                                                                                               ttCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGATg        >  2:1267841/1‑150 (MQ=255)
                                                                                                                                                  cGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGATGCGa     <  2:721005/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGATTAGAACGCGGTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGATGCGATGG  >  CP009273/3882611‑3882909

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: