Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 364,561 G→T 100% D18Y (GAC→TAC)  mhpB → 2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273364,5610GT95.7% 637.2 / 16.1 222D18Y (GAC→TAC) mhpB2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (4/6);  new base T (110/102);  total (114/108)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 222 total aligned reads displayed.

CGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTC  >  CP009273/364417‑364704
                                                                                                                                                |                                                                                                                                               
    aaCCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGgc                                                                                                                                          >  2:3091826/1‑148 (MQ=255)
     aCCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGgcg                                                                                                                                         >  1:502695/1‑148 (MQ=255)
     aCCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGgcg                                                                                                                                         <  2:3211387/148‑1 (MQ=255)
       ccTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCa                                                                                                                                       <  1:1221479/148‑1 (MQ=255)
       ccTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCa                                                                                                                                       <  2:2907184/148‑1 (MQ=255)
         tGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCacg                                                                                                                                     >  1:2398855/1‑146 (MQ=255)
         tGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAag                                                                                                                                     <  2:654552/148‑1 (MQ=255)
            gCAAGACCCTGAATACACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAgg                                                                                                                                  <  2:3334745/148‑1 (MQ=255)
            gCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAgg                                                                                                                                  <  1:220129/148‑1 (MQ=255)
             cAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGt                                                                                                                                 <  2:2148109/148‑1 (MQ=255)
              aaGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTg                                                                                                                                <  2:3440130/148‑1 (MQ=255)
               aGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGc                                                                                                                               <  1:335267/148‑1 (MQ=255)
               aGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGc                                                                                                                               <  1:1638840/148‑1 (MQ=255)
                gACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCt                                                                                                                              >  2:2308190/1‑148 (MQ=255)
                 aCCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTc                                                                                                                             >  1:3237686/1‑148 (MQ=255)
                  cccTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCg                                                                                                                            <  1:1931013/148‑1 (MQ=255)
                  cccTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCg                                                                                                                            <  2:2245233/148‑1 (MQ=255)
                       aaTAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAg                                                                                                                       >  2:2424155/1‑148 (MQ=255)
                            aCTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCaa                                                                                                                  >  2:461740/1‑148 (MQ=255)
                                 cGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCg                                                                                                             <  2:336556/148‑1 (MQ=255)
                                  gTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTNACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGt                                                                                                            <  1:1806822/148‑1 (MQ=255)
                                    cGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGa                                                                                                          <  1:3127203/148‑1 (MQ=255)
                                       tgatgaCGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTg                                                                                                       >  2:1638934/1‑148 (MQ=255)
                                       tgatgaCGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTg                                                                                                       >  1:91219/1‑148 (MQ=255)
                                         atgaCACTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGcc                                                                                                     >  1:834829/1‑148 (MQ=255)
                                             cGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCTCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAgcg                                                                                                 <  1:1943128/148‑1 (MQ=255)
                                                tGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGccc                                                                                              >  2:1441946/1‑148 (MQ=255)
                                                  aCCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCcgc                                                                                            >  2:3313024/1‑148 (MQ=255)
                                                     cGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAg                                                                                         >  2:1627528/1‑148 (MQ=255)
                                                         cTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTa                                                                                     >  1:3432008/1‑148 (MQ=255)
                                                         cTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTa                                                                                     <  1:106417/148‑1 (MQ=255)
                                                         cTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTa                                                                                     <  2:1489555/148‑1 (MQ=255)
                                                            aTGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTg                                                                                  >  1:3347360/1‑148 (MQ=255)
                                                              gCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCg                                                                                >  2:2857632/1‑148 (MQ=255)
                                                                cGACGTGTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGc                                                                              <  2:2312831/148‑1 (MQ=255)
                                                                 gACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCa                                                                             <  1:57832/148‑1 (MQ=255)
                                                                   cGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCAtt                                                                           <  2:704526/148‑1 (MQ=255)
                                                                     tttCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATtct                                                                         <  2:3158454/148‑1 (MQ=255)
                                                                     tttCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATtct                                                                         >  1:530439/1‑148 (MQ=255)
                                                                        cTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTccc                                                                      <  2:3010260/148‑1 (MQ=255)
                                                                         tGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTcccc                                                                     <  1:3308060/148‑1 (MQ=255)
                                                                         tGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTcccc                                                                     <  1:3308163/148‑1 (MQ=255)
                                                                          gTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCAGTCTCCCCt                                                                    >  1:66864/1‑148 (MQ=255)
                                                                              aaaaGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAAc                                                                >  1:3411866/1‑148 (MQ=255)
                                                                              aaaaGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAAc                                                                >  1:1652332/1‑148 (MQ=255)
                                                                               aaaGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACt                                                               >  1:1583586/1‑148 (MQ=255)
                                                                                aaGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtg                                                              >  1:1995663/1‑148 (MQ=255)
                                                                                 aGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtgg                                                             <  2:1290555/148‑1 (MQ=255)
                                                                                  ggTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCAGTCTCCCCTGAACtggt                                                            >  2:3060473/1‑148 (MQ=255)
                                                                                   gTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtggtg                                                           >  2:946020/1‑148 (MQ=255)
                                                                                   gTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtggtg                                                           <  1:3042370/148‑1 (MQ=255)
                                                                                    tAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGTGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtggtgg                                                          <  2:2944913/148‑1 (MQ=255)
                                                                                       ccTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGc                                                       <  2:596399/148‑1 (MQ=255)
                                                                                       ccTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGc                                                       <  1:2119801/148‑1 (MQ=255)
                                                                                       ccTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGc                                                       <  2:836097/148‑1 (MQ=255)
                                                                                         tGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTg                                                     <  1:1738759/148‑1 (MQ=255)
                                                                                         tGATATGCACGCTTATCTTCACAGTCTATCCCACTCGCCGCTGGTGGGGTATGACTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTg                                                     >  2:1105750/1‑148 (MQ=255)
                                                                                          gATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGt                                                    <  2:2054002/148‑1 (MQ=255)
                                                                                            tatGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGttt                                                  <  2:2186067/148‑1 (MQ=255)
                                                                                            tatGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGttt                                                  >  2:2837360/1‑148 (MQ=255)
                                                                                               gCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGATCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTgcg                                               >  1:249438/1‑148 (MQ=255)
                                                                                                cACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTgcgc                                              >  2:454747/1‑148 (MQ=255)
                                                                                                 aCGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGcc                                             <  1:2677070/148‑1 (MQ=255)
                                                                                                  cGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCa                                            >  1:2448763/1‑148 (MQ=255)
                                                                                                   gCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTTAACTGGTGGTGCTGTTTGCGCCAg                                           >  1:2373152/1‑148 (MQ=255)
                                                                                                       aTCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAAGGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCa                                       <  2:1914195/148‑1 (MQ=255)
                                                                                                         cTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACt                                     <  2:2763146/148‑1 (MQ=255)
                                                                                                             aCTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACaa                                 >  2:1222119/1‑148 (MQ=255)
                                                                                                              cTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAAc                                >  1:1866169/1‑148 (MQ=255)
                                                                                                              cTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAAc                                <  2:1571623/148‑1 (MQ=255)
                                                                                                              cTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAAc                                >  1:2753306/1‑148 (MQ=255)
                                                                                                              cTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAAc                                <  2:1268148/148‑1 (MQ=255)
                                                                                                               tGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACg                               <  1:115885/148‑1 (MQ=255)
                                                                                                               tGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACg                               >  2:1340075/1‑148 (MQ=255)
                                                                                                                 tCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGc                             >  1:1049274/1‑148 (MQ=255)
                                                                                                                   tttCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGGGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCtt                           <  1:900194/148‑1 (MQ=255)
                                                                                                                   tttCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCtt                           <  2:1861712/148‑1 (MQ=255)
                                                                                                                   tttCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCtt                           <  1:2412896/148‑1 (MQ=255)
                                                                                                                    ttCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCttt                          >  2:1127694/1‑148 (MQ=255)
                                                                                                                    ttCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCttt                          >  1:1364379/1‑148 (MQ=255)
                                                                                                                       ccACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTc                       <  2:2571808/148‑1 (MQ=255)
                                                                                                                       ccACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTc                       >  1:2883486/1‑148 (MQ=255)
                                                                                                                       ccACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTc                       >  2:2299442/1‑148 (MQ=255)
                                                                                                                         aCTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTa                     <  1:2320464/148‑1 (MQ=255)
                                                                                                                         aCTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTa                     <  1:2304148/148‑1 (MQ=255)
                                                                                                                         aCTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTa                     <  1:2857632/148‑1 (MQ=255)
                                                                                                                         aCTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTa                     <  1:2857328/148‑1 (MQ=255)
                                                                                                                          cTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTAt                    >  1:2214223/1‑148 (MQ=255)
                                                                                                                          cTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTAt                    >  2:246496/1‑148 (MQ=255)
                                                                                                                             gccgcTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGAc                 >  1:1431375/1‑148 (MQ=255)
                                                                                                                             gccgcTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGAc                 >  1:1228988/1‑148 (MQ=255)
                                                                                                                              ccgcTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACg                <  2:2398855/148‑1 (MQ=255)
                                                                                                                                   ggtggGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATg           <  2:3301912/148‑1 (MQ=255)
                                                                                                                                   ggtggGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGAATGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATg           >  2:487405/1‑148 (MQ=255)
                                                                                                                                    gtggGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGc          >  2:1678237/1‑148 (MQ=255)
                                                                                                                                      ggggTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCa        >  1:369328/1‑148 (MQ=255)
                                                                                                                                          tATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGt    <  2:2577803/148‑1 (MQ=255)
                                                                                                                                           aTGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGtt   <  2:3268236/148‑1 (MQ=255)
                                                                                                                                            tGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTc  <  2:2505638/148‑1 (MQ=255)
                                                                                                                                            tGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTc  <  2:1065323/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
CGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTC  >  CP009273/364417‑364704

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: