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,5610GT96.2% 565.6 / 12.0 195D18Y (GAC→TAC) mhpB2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (4/4);  new base T (105/82);  total (109/86)
Fisher's exact test for biased strand distribution p-value = 7.33e-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 195 total aligned reads displayed.

CGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCT  >  CP009273/364417‑364705
                                                                                                                                                |                                                                                                                                                
  cAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCg                                                                                                                                             >  1:3299682/1‑148 (MQ=255)
   aaaCCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCgg                                                                                                                                            >  1:1683275/1‑148 (MQ=255)
     aCCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGgcg                                                                                                                                          >  2:2316515/1‑148 (MQ=255)
        cTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTAGCCGCTGGTGGGGTATGTCTACCCGGCGCaa                                                                                                                                       >  2:829840/1‑148 (MQ=255)
        cTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTAGCCGCTGGTGGGGTATGTCTACCCGGCGCaa                                                                                                                                       >  2:1030879/1‑148 (MQ=255)
          gggCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAaga                                                                                                                                     >  1:1431148/1‑148 (MQ=255)
           ggCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAagag                                                                                                                                    >  1:100862/1‑148 (MQ=255)
            gCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAgg                                                                                                                                   <  1:224964/148‑1 (MQ=255)
            gCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAgg                                                                                                                                   <  1:2612149/148‑1 (MQ=255)
            gCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAgg                                                                                                                                   <  1:2861407/148‑1 (MQ=255)
            gCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAgg                                                                                                                                   <  1:1616068/148‑1 (MQ=255)
             cAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGt                                                                                                                                  >  2:1372223/1‑148 (MQ=255)
                   ccTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGa                                                                                                                            <  2:1284132/148‑1 (MQ=255)
                    cTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGAt                                                                                                                           >  2:2721691/1‑148 (MQ=255)
                     tGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATg                                                                                                                          >  1:1010540/1‑148 (MQ=255)
                      gAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGa                                                                                                                         >  2:645539/1‑148 (MQ=255)
                       aaTAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAg                                                                                                                        <  2:1049160/148‑1 (MQ=255)
                        aTAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAgg                                                                                                                       >  1:2929164/1‑148 (MQ=255)
                         tAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGt                                                                                                                      >  1:2307351/1‑148 (MQ=255)
                         tAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGt                                                                                                                      >  1:382581/1‑148 (MQ=255)
                          aaaCTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTc                                                                                                                     >  2:2289783/1‑148 (MQ=255)
                          aaaCTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTc                                                                                                                     >  1:1732308/1‑148 (MQ=255)
                           aaCTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCa                                                                                                                    <  1:112225/148‑1 (MQ=255)
                            aCTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCaa                                                                                                                   >  1:3299445/1‑148 (MQ=255)
                            aCTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCaa                                                                                                                   >  1:2988045/1‑148 (MQ=255)
                             cTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAAt                                                                                                                  <  1:2205976/148‑1 (MQ=255)
                               ggCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATgg                                                                                                                <  1:1109699/148‑1 (MQ=255)
                               ggCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATgg                                                                                                                <  1:124210/148‑1 (MQ=255)
                                  gTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGt                                                                                                             <  2:3488855/148‑1 (MQ=255)
                                     ggTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGAt                                                                                                          <  1:18936/148‑1 (MQ=255)
                                       tgatgaCGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTg                                                                                                        >  1:2431042/1‑148 (MQ=255)
                                       tgatgaCGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTg                                                                                                        >  1:3232933/1‑148 (MQ=255)
                                        gatgaCGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGTGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGc                                                                                                       <  2:3252066/148‑1 (MQ=255)
                                        gatgaCGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAATAGGTGCTCGATGAGGTCAATGGCGTGATTGc                                                                                                       <  2:410029/148‑1 (MQ=255)
                                         atgaCGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGcc                                                                                                      <  1:3391190/148‑1 (MQ=255)
                                          tgaCGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTTATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCa                                                                                                     >  2:3121243/1‑148 (MQ=255)
                                                  aCCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCcgc                                                                                             <  1:95607/148‑1 (MQ=255)
                                                    ccGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGa                                                                                           >  1:2089946/1‑148 (MQ=255)
                                                          tGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTAt                                                                                     <  1:1893836/148‑1 (MQ=255)
                                                          tGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTAt                                                                                     >  1:3158467/1‑148 (MQ=255)
                                                          tGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTAt                                                                                     >  1:2377997/1‑148 (MQ=255)
                                                                   cGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCAtt                                                                            <  1:602539/148‑1 (MQ=255)
                                                                     tttCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATtct                                                                          >  1:2572921/1‑148 (MQ=255)
                                                                      ttCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATtctc                                                                         >  2:3439649/1‑148 (MQ=255)
                                                                          gTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCt                                                                     >  2:251030/1‑148 (MQ=255)
                                                                           tCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTg                                                                    >  2:892475/1‑148 (MQ=255)
                                                                            cGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGa                                                                   >  2:1201452/1‑148 (MQ=255)
                                                                                caGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtg                                                               <  2:2645105/147‑1 (MQ=255)
                                                                                aaGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtg                                                               >  1:1253892/1‑148 (MQ=255)
                                                                                  ggTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtggt                                                             >  2:1543361/1‑148 (MQ=255)
                                                                                         tGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTg                                                      >  1:10641/1‑148 (MQ=255)
                                                                                         tGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGACATTCTCCCCTGAACTGGTGGTGCTg                                                      >  1:183393/1‑148 (MQ=255)
                                                                                          gATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGt                                                     >  1:34391/1‑148 (MQ=255)
                                                                                           atatGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGtt                                                    <  2:738867/148‑1 (MQ=255)
                                                                                           atatGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGtt                                                    >  2:516694/1‑148 (MQ=255)
                                                                                            tatGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGttt                                                   >  1:3517996/1‑148 (MQ=255)
                                                                                                cACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTgcgc                                               >  2:23433/1‑148 (MQ=255)
                                                                                                  cGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCa                                             >  2:2681713/1‑148 (MQ=255)
                                                                                                  cGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCa                                             >  2:3705566/1‑148 (MQ=255)
                                                                                                    cTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGa                                           <  1:3584374/148‑1 (MQ=255)
                                                                                                     ttATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGAt                                          >  1:3471565/1‑148 (MQ=255)
                                                                                                     ttATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGAt                                          >  1:102648/1‑148 (MQ=255)
                                                                                                       aTCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCa                                        <  2:1732308/148‑1 (MQ=255)
                                                                                                        tCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCAc                                       >  2:1940016/1‑148 (MQ=255)
                                                                                                          ttCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTa                                     <  2:3299682/148‑1 (MQ=255)
                                                                                                          ttCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTa                                     <  2:3219646/148‑1 (MQ=255)
                                                                                                           tCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTAc                                    >  2:1131169/1‑148 (MQ=255)
                                                                                                           tCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTAc                                    >  1:1605611/1‑148 (MQ=255)
                                                                                                            cACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACa                                   <  2:374959/148‑1 (MQ=255)
                                                                                                            cACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACa                                   >  2:2754036/1‑148 (MQ=255)
                                                                                                             aCTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACaa                                  <  2:2089946/148‑1 (MQ=255)
                                                                                                             aCTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACaa                                  >  2:2121649/1‑148 (MQ=255)
                                                                                                               tGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACg                                >  1:2404435/1‑148 (MQ=255)
                                                                                                                gTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACgg                               >  1:834855/1‑148 (MQ=255)
                                                                                                                  cTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCt                             <  2:3299445/148‑1 (MQ=255)
                                                                                                                  cTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCt                             <  2:2377578/148‑1 (MQ=255)
                                                                                                                    ttCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCttt                           <  2:78506/148‑1 (MQ=255)
                                                                                                                    ttCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCttt                           >  1:111399/1‑148 (MQ=255)
                                                                                                                     tCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCtttt                          <  2:39094/148‑1 (MQ=255)
                                                                                                                      cccACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCttttt                         <  1:3200499/148‑1 (MQ=255)
                                                                                                                        cACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCt                       >  2:3116902/1‑148 (MQ=255)
                                                                                                                        cACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCt                       >  2:3117508/1‑148 (MQ=255)
                                                                                                                         aCTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTa                      >  2:2378455/1‑148 (MQ=255)
                                                                                                                         aCTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTa                      <  1:1166441/148‑1 (MQ=255)
                                                                                                                         aCTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTa                      <  1:3072424/148‑1 (MQ=255)
                                                                                                                         aCTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTa                      <  2:1773388/148‑1 (MQ=255)
                                                                                                                         aCTAGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTa                      <  1:1030879/148‑1 (MQ=255)
                                                                                                                          cTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTAt                     >  1:2506638/1‑148 (MQ=255)
                                                                                                                          cTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTAt                     >  1:2665346/1‑148 (MQ=255)
                                                                                                                          cTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTAt                     >  2:3186808/1‑148 (MQ=255)
                                                                                                                            cgccgcTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGa                   <  2:3126657/148‑1 (MQ=255)
                                                                                                                             gccgcTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGAc                  <  2:1239827/148‑1 (MQ=255)
                                                                                                                             gccgcTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGAc                  <  2:290747/148‑1 (MQ=255)
                                                                                                                              ccgcTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACg                 <  2:3108486/148‑1 (MQ=255)
                                                                                                                                 cTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGa              <  2:3566896/148‑1 (MQ=255)
                                                                                                                                 cTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGa              <  1:378466/148‑1 (MQ=255)
                                                                                                                                   ggtggGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATg            >  2:1413014/1‑148 (MQ=255)
                                                                                                                                       gggTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCAc        <  2:100862/148‑1 (MQ=255)
                                                                                                                                            tGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTc   <  2:2848663/148‑1 (MQ=255)
                                                                                                                                            tGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTc   <  1:5373/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCT  >  CP009273/364417‑364705

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: