Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,882,760 T→A 5.6% 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 5.6% 1335.6 / 51.2 388V224E (GTG→GAG) tnaAtryptophanase/L‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (194/170);  new base A (7/15);  total (203/185)
Fisher's exact test for biased strand distribution p-value = 7.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 399 total aligned reads displayed.

GATTAGAACGCGGTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGATGCGATGG  >  CP009273/3882613‑3882909
                                                                                                                                                   |                                                                                                                                                     
   tAGAACGCGGTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTaa                                                                                                                                                  >  1:5222905/1‑150 (MQ=255)
   tAGAACGCGGTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTaa                                                                                                                                                  >  2:4215033/1‑150 (MQ=255)
   tAGAACGCGGTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGGTGCAACCATCACCAGTAACTATGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAACTACGATATTCCGGTGGTaa                                                                                                                                                  <  2:499442/150‑1 (MQ=255)
    aGAACGCGGTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAATCGATGTACAGCATCGCGAAGAAATAAGATATTCCGGTGGTCAt                                                                                                                                                 >  1:3936439/1‑150 (MQ=255)
    aGAACGCGGTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAAt                                                                                                                                                 >  2:2166757/1‑150 (MQ=255)
     gAACGCGGTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATg                                                                                                                                                >  2:5244669/1‑150 (MQ=255)
           ggTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTcc                                                                                                                                          >  1:383740/1‑150 (MQ=255)
           ggTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTcc                                                                                                                                          >  1:779217/1‑150 (MQ=255)
           ggTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTcc                                                                                                                                          >  2:3906804/1‑150 (MQ=255)
           ggTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTcc                                                                                                                                          >  2:2325785/1‑150 (MQ=255)
            gTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTTGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCcg                                                                                                                                         >  1:45077/1‑150 (MQ=255)
             tATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCcgc                                                                                                                                        >  1:2838968/1‑150 (MQ=255)
                  aagaagTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCt                                                                                                                                   >  2:3189693/1‑150 (MQ=255)
                   agaagTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCgcgcgcgt                                                                                                                                  <  2:2669153/150‑3 (MQ=255)
                    taagTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCttt                                                                                                                                 <  2:2205746/149‑1 (MQ=255)
                    gaagTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCttt                                                                                                                                 >  2:3150321/1‑150 (MQ=255)
                      agTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATAACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGc                                                                                                                               >  1:2375239/1‑150 (MQ=255)
                         tGGTCCGAATTACGTGCCGTATATCGTTGCAACCATCACAAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGa                                                                                                                            <  2:3800719/150‑1 (MQ=255)
                         tGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGa                                                                                                                            <  2:3457814/150‑1 (MQ=255)
                          ggTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGaa                                                                                                                           >  2:1934983/1‑150 (MQ=255)
                            tCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGaaaa                                                                                                                         >  1:108518/1‑150 (MQ=255)
                            tCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGaaaa                                                                                                                         >  1:852253/1‑150 (MQ=255)
                             ccGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAAc                                                                                                                        <  1:5627045/150‑1 (MQ=255)
                               gAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGc                                                                                                                      >  1:4956398/1‑150 (MQ=255)
                               gAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGc                                                                                                                      >  2:4317855/1‑150 (MQ=255)
                                 aTAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCt                                                                                                                    <  1:3111714/150‑1 (MQ=255)
                                   aaCGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTAt                                                                                                                  >  1:3905485/1‑150 (MQ=255)
                                             atatCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAagc                                                                                                        >  1:28614/1‑150 (MQ=255)
                                             atatCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAagc                                                                                                        >  2:1374238/1‑150 (MQ=255)
                                             atatCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAagc                                                                                                        >  2:4548685/1‑150 (MQ=255)
                                              tatCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAagca                                                                                                       <  1:3711527/150‑1 (MQ=255)
                                                   ttGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTg                                                                                                  >  2:4474200/1‑150 (MQ=255)
                                                      cAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAg                                                                                               >  1:1947755/1‑150 (MQ=255)
                                                       aaCCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGc                                                                                              >  1:1498681/1‑150 (MQ=255)
                                                       aaCCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGc                                                                                              >  1:4932698/1‑150 (MQ=255)
                                                        aCCATCACCAGTAACTCTGCAGGTTGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACTGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCa                                                                                             <  1:1665118/150‑1 (MQ=255)
                                                           aTCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACTATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGaa                                                                                          <  2:4759675/150‑1 (MQ=255)
                                                           aTCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGaa                                                                                          <  2:323019/150‑1 (MQ=255)
                                                           aTCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGaa                                                                                          >  1:3542647/1‑150 (MQ=255)
                                                              aCCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATCCGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATAc                                                                                       <  2:561410/150‑1 (MQ=255)
                                                               ccAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACa                                                                                      <  2:3709567/150‑1 (MQ=255)
                                                                cAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACaa                                                                                     >  1:4116102/1‑150 (MQ=255)
                                                                 aGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACaaa                                                                                    >  1:503289/1‑150 (MQ=255)
                                                                 aGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACaaa                                                                                    >  2:2763230/1‑150 (MQ=255)
                                                                  gTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAg                                                                                   <  2:3869361/150‑1 (MQ=255)
                                                                  gTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAg                                                                                   <  1:652226/150‑1 (MQ=255)
                                                                   tAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGa                                                                                  >  1:1859330/1‑150 (MQ=255)
                                                                    aaCTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGAc                                                                                 >  2:3532324/1‑150 (MQ=255)
                                                                      ctctGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAAGACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTg                                                                               <  2:4025614/150‑1 (MQ=255)
                                                                          gCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGAcc                                                                           >  2:4638559/1‑150 (MQ=255)
                                                                            aGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCAt                                                                         <  1:4022831/150‑1 (MQ=255)
                                                                             ggtggtCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGACTACAAAGACTGGACCATc                                                                        <  1:4874279/150‑1 (MQ=255)
                                                                             ggtggtCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATc                                                                        <  2:576092/150‑1 (MQ=255)
                                                                                ggtCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAg                                                                     <  1:2488320/150‑1 (MQ=255)
                                                                                 gtCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGc                                                                    >  1:4370941/1‑150 (MQ=255)
                                                                                   cAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGAGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAg                                                                  >  1:4241063/1‑150 (MQ=255)
                                                                                    aGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGa                                                                 >  1:4700119/1‑150 (MQ=255)
                                                                                     gCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGAt                                                                >  1:421361/1‑150 (MQ=255)
                                                                                      ccGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATc                                                               <  1:4722709/150‑1 (MQ=255)
                                                                                      ccGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATc                                                               <  1:4667788/150‑1 (MQ=255)
                                                                                        ggTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCAc                                                             >  2:1870902/1‑150 (MQ=255)
                                                                                          tttCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCAccc                                                           <  1:1004884/150‑1 (MQ=255)
                                                                                           ttCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCcg                                                          <  1:923373/150‑1 (MQ=255)
                                                                                             cACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCcgcg                                                        >  1:3393428/1‑150 (MQ=255)
                                                                                               cTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGaa                                                      <  2:1271325/150‑1 (MQ=255)
                                                                                                   cAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGAGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCt                                                  <  1:1620269/150‑1 (MQ=255)
                                                                                                    aaaCTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTa                                                 <  2:1411722/150‑1 (MQ=255)
                                                                                                     aaCTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTAc                                                <  2:1195515/150‑1 (MQ=255)
                                                                                                      aCTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGACTCCAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACa                                               >  1:2336695/1‑150 (MQ=255)
                                                                                                      aCTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACa                                               <  2:2390999/150‑1 (MQ=255)
                                                                                                      aCTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACa                                               <  1:219149/150‑1 (MQ=255)
                                                                                                          aaaaGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGAGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAata                                           <  2:3651611/150‑1 (MQ=255)
                                                                                                           aaaGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAatat                                          >  2:273789/1‑150 (MQ=255)
                                                                                                           aaaGCGATGGACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAatat                                          <  2:4376182/150‑1 (MQ=255)
                                                                                                               cGATGTACAGCATCGCGAAGAAATACGATATTCCGGAGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCg                                      <  1:2659472/150‑1 (MQ=255)
                                                                                                                 aTGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGat                                    <  2:209316/150‑1 (MQ=255)
                                                                                                                 aTGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGat                                    >  2:4267660/1‑150 (MQ=255)
                                                                                                                  tGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGata                                   <  2:1635394/150‑1 (MQ=255)
                                                                                                                       aGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTg                              <  1:428907/150‑1 (MQ=255)
                                                                                                                        gCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTgg                             <  2:2520610/150‑1 (MQ=255)
                                                                                                                         cATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATTCCGATATGCTGGc                            >  1:3475132/1‑150 (MQ=255)
                                                                                                                         cATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGc                            <  1:2017376/150‑1 (MQ=255)
                                                                                                                             gcgAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATg                        <  2:3050696/150‑1 (MQ=255)
                                                                                                                               gaagaaATACGATATTACGGTGGTAATGGACTCCGCGCTCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTc                      <  2:5229998/150‑1 (MQ=255)
                                                                                                                                aagaaATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTcc                     <  1:2634319/150‑1 (MQ=255)
                                                                                                                                   aaaTACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGcc                  >  2:910321/1‑150 (MQ=255)
                                                                                                                                      tACGATATTCCGGTGGTAATGGTCTCCGCGCGCTTTGCTTAAATCGCCTAGTTCATCAAGCAGCGTGACGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATCTGCCGATATGCTGGCGATGTCCGCCaag               <  1:3830643/150‑1 (MQ=255)
                                                                                                                                       aCGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCaaga              >  1:2747068/1‑150 (MQ=255)
                                                                                                                                       aCGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCaaga              >  2:806303/1‑150 (MQ=255)
                                                                                                                                         gATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCaagaaa            >  2:102630/1‑150 (MQ=255)
                                                                                                                                         gATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCaagaaa            >  2:5336583/1‑150 (MQ=255)
                                                                                                                                          atatTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCaagaaag           >  2:2022130/1‑150 (MQ=255)
                                                                                                                                           tatTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCaagaaaga          >  1:5147397/1‑150 (MQ=255)
                                                                                                                                             ttCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGATg        >  1:5539070/1‑150 (MQ=255)
                                                                                                                                              tCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGATGc       <  1:5426039/150‑1 (MQ=255)
                                                                                                                                              tCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGATGc       <  1:1849015/150‑1 (MQ=255)
                                                                                                                                               ccGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGATGCg      <  2:4894824/150‑1 (MQ=255)
                                                                                                                                                 ggtggtAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGATGCGAt    >  2:4180003/1‑150 (MQ=255)
                                                                                                                                                 ggtggtAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGATGCGAt    >  2:3838721/1‑150 (MQ=255)
                                                                                                                                                  gtggtAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGATGCGATg   >  2:4821340/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GATTAGAACGCGGTATTGAAGAAGTTGGTCCGAATAACGTGCCGTATATCGTTGCAACCATCACCAGTAACTCTGCAGGTGGTCAGCCGGTTTCACTGGCAAACTTAAAAGCGATGTACAGCATCGCGAAGAAATACGATATTCCGGTGGTAATGGACTCCGCGCGCTTTGCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGATGCGATGG  >  CP009273/3882613‑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: