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,5610GT97.4% 590.4 / 6.1 195D18Y (GAC→TAC) mhpB2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base T (103/87);  total (106/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 195 total aligned reads displayed.

CCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTC  >  CP009273/364416‑364704
                                                                                                                                                 |                                                                                                                                               
 cGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTAcc                                                                                                                                              <  2:884400/148‑1 (MQ=255)
    aaaCCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCgg                                                                                                                                           >  1:1898720/1‑148 (MQ=255)
     aaCCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGgc                                                                                                                                          >  2:2384734/1‑148 (MQ=255)
      aCCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGgcg                                                                                                                                         <  2:1040128/148‑1 (MQ=255)
       cccTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGgcgc                                                                                                                                        >  2:12035/1‑148 (MQ=255)
        ccTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCa                                                                                                                                       >  2:2033638/1‑148 (MQ=255)
        ccTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTAACTGTCTTTCCCACTCGCCGCTGGTGGGGGATGTCTACCCCGCGCa                                                                                                                                       >  2:2339024/1‑148 (MQ=255)
         cTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCACCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCaa                                                                                                                                      >  1:2280575/1‑148 (MQ=255)
          tGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAag                                                                                                                                     >  1:2063876/1‑148 (MQ=255)
          tGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAag                                                                                                                                     >  1:2290839/1‑148 (MQ=255)
          tGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAag                                                                                                                                     >  1:2176288/1‑148 (MQ=255)
            ggCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCATGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAagag                                                                                                                                   <  1:2326579/148‑1 (MQ=255)
               aaGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTg                                                                                                                                <  1:2389396/148‑1 (MQ=255)
                aGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGc                                                                                                                               <  2:1515681/148‑1 (MQ=255)
                 gACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCt                                                                                                                              >  2:3033055/1‑148 (MQ=255)
                  aCCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGAAGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTAACTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTAAGTCTACCCGGCGCAAGAGGTGCTc                                                                                                                             >  2:3080513/1‑148 (MQ=255)
                    ccTGAATAACCTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGa                                                                                                                           <  2:2438868/148‑1 (MQ=255)
                     cTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGAt                                                                                                                          >  2:2108380/1‑148 (MQ=255)
                      tGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATg                                                                                                                         <  1:2967000/148‑1 (MQ=255)
                      tGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATg                                                                                                                         <  1:1215795/148‑1 (MQ=255)
                         aTAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAgg                                                                                                                      <  2:86856/148‑1 (MQ=255)
                          tAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGt                                                                                                                     >  1:3200997/1‑148 (MQ=255)
                           aaaCTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTc                                                                                                                    >  2:1395436/1‑148 (MQ=255)
                                  cGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCg                                                                                                             >  2:922374/1‑148 (MQ=255)
                                    tCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTg                                                                                                           <  1:913353/148‑1 (MQ=255)
                                     cGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGa                                                                                                          >  1:2692164/1‑148 (MQ=255)
                                     cGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGa                                                                                                          >  1:938789/1‑148 (MQ=255)
                                         gatgaCGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGc                                                                                                      <  2:1968892/148‑1 (MQ=255)
                                          atgaCGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGcc                                                                                                     >  2:2938481/1‑148 (MQ=255)
                                            gaCGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAg                                                                                                   >  2:1992019/1‑148 (MQ=255)
                                            gaCGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAg                                                                                                   >  2:1122075/1‑148 (MQ=255)
                                                cTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGcc                                                                                               <  1:2489268/148‑1 (MQ=255)
                                                 tGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGccc                                                                                              >  1:3418518/1‑148 (MQ=255)
                                                  gACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCcg                                                                                             >  2:81362/1‑148 (MQ=255)
                                                   aCCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCcgc                                                                                            >  2:2782276/1‑148 (MQ=255)
                                                    cccGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCcgcg                                                                                           >  2:1775143/1‑148 (MQ=255)
                                                     ccGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGa                                                                                          >  1:1129230/1‑148 (MQ=255)
                                                             aTGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTAGTg                                                                                  <  2:1829036/148‑1 (MQ=255)
                                                              tGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGc                                                                                 >  2:1481334/1‑148 (MQ=255)
                                                              tGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGc                                                                                 >  2:2884972/1‑148 (MQ=255)
                                                                 cGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGc                                                                              <  1:2774110/148‑1 (MQ=255)
                                                                   aCGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCAt                                                                            >  1:260021/1‑148 (MQ=255)
                                                                    cGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCAtt                                                                           >  1:1891357/1‑148 (MQ=255)
                                                                      tttCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATtct                                                                         >  1:969549/1‑148 (MQ=255)
                                                                       ttCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATtctc                                                                        <  2:1255780/148‑1 (MQ=255)
                                                                       ttCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATtctc                                                                        >  1:196416/1‑148 (MQ=255)
                                                                       ttCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATtctc                                                                        >  2:3299049/1‑148 (MQ=255)
                                                                        tCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTcc                                                                       >  2:1385558/1‑148 (MQ=255)
                                                                        tCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTcc                                                                       >  2:1375233/1‑148 (MQ=255)
                                                                         cTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTccc                                                                      <  1:2865402/148‑1 (MQ=255)
                                                                           gTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGGATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCt                                                                    >  1:2664579/1‑148 (MQ=255)
                                                                              gAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGaa                                                                 <  2:427400/148‑1 (MQ=255)
                                                                              gAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGaa                                                                 <  1:2641034/148‑1 (MQ=255)
                                                                               aaaaGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAAc                                                                <  2:2061509/148‑1 (MQ=255)
                                                                                 aaGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtg                                                              <  1:33917/148‑1 (MQ=255)
                                                                                   ggTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtggt                                                            >  2:971620/1‑148 (MQ=255)
                                                                                     tAGCCTGATAGGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtggtgg                                                          <  1:410900/148‑1 (MQ=255)
                                                                                       gCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGGCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGgtggtg                                                        <  2:767040/148‑1 (MQ=255)
                                                                                       gCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGgtggtg                                                        <  2:2169084/148‑1 (MQ=255)
                                                                                        ccTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGc                                                       >  2:1576565/1‑148 (MQ=255)
                                                                                         cTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCt                                                      >  2:2469490/1‑148 (MQ=255)
                                                                                           gATATGCACGCTTATCTTCTCTGTCTTTCCAACTCGCCGCTGGTGGGGTATGTCTACCAGGCTCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGTCATTCTCCCCTGAACTGGTGGTGCTGt                                                    <  1:358857/148‑1 (MQ=255)
                                                                                               tGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTgc                                                >  1:1673524/1‑148 (MQ=255)
                                                                                                  aCGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGcc                                             <  2:3200997/148‑1 (MQ=255)
                                                                                                    gCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAg                                           <  1:430193/148‑1 (MQ=255)
                                                                                                     cTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCTGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGa                                          <  1:977342/148‑1 (MQ=255)
                                                                                                     cTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGa                                          <  2:2511453/148‑1 (MQ=255)
                                                                                                     cTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGa                                          <  2:2510900/148‑1 (MQ=255)
                                                                                                     cTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGa                                          <  2:3418518/148‑1 (MQ=255)
                                                                                                         tCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCAc                                      >  2:425236/1‑148 (MQ=255)
                                                                                                         tCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCAc                                      >  2:1167881/1‑148 (MQ=255)
                                                                                                          cTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACt                                     >  1:2992353/1‑148 (MQ=255)
                                                                                                           ttCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTa                                    <  1:524586/148‑1 (MQ=255)
                                                                                                           ttCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTa                                    >  2:160031/1‑148 (MQ=255)
                                                                                                             cACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACa                                  <  2:2744614/148‑1 (MQ=255)
                                                                                                             cACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACa                                  <  2:2744562/148‑1 (MQ=255)
                                                                                                             cACTGTCTGTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACa                                  <  2:2613244/148‑1 (MQ=255)
                                                                                                              aCTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACaa                                 <  1:2384734/148‑1 (MQ=255)
                                                                                                              aCTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACaa                                 <  1:621357/148‑1 (MQ=255)
                                                                                                                tGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACg                               >  1:1293449/1‑148 (MQ=255)
                                                                                                                 gTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACgg                              >  1:1060429/1‑148 (MQ=255)
                                                                                                                  tCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTGTGCGCCAGATCACTACAACGGc                             >  2:460183/1‑148 (MQ=255)
                                                                                                                   cTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCt                            <  1:236541/148‑1 (MQ=255)
                                                                                                                   cTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCt                            >  1:207524/1‑148 (MQ=255)
                                                                                                                    tttCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCtt                           <  2:2686234/148‑1 (MQ=255)
                                                                                                                     ttCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCttt                          <  1:1801182/148‑1 (MQ=255)
                                                                                                                      tCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCtttt                         >  1:3303761/1‑148 (MQ=255)
                                                                                                                       cccACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCttttt                        >  1:975707/1‑148 (MQ=255)
                                                                                                                        ccACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTc                       <  1:2624417/148‑1 (MQ=255)
                                                                                                                          aCTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTa                     >  2:1558187/1‑148 (MQ=255)
                                                                                                                           cTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTAt                    >  1:1605043/1‑148 (MQ=255)
                                                                                                                             cgccgcTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGa                  >  1:274050/1‑148 (MQ=255)
                                                                                                                             cgccgcTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGa                  >  1:77784/1‑148 (MQ=255)
                                                                                                                             cgccgcTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGa                  <  1:2033638/148‑1 (MQ=255)
                                                                                                                               ccgcTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACg                <  1:1339117/148‑1 (MQ=255)
                                                                                                                                 gcTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTg              <  2:2254038/148‑1 (MQ=255)
                                                                                                                                   tggtggGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGAt            <  1:2782276/148‑1 (MQ=255)
                                                                                                                                     gtggGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGc          >  1:106000/1‑148 (MQ=255)
                                                                                                                                       ggggTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCa        <  1:2634464/148‑1 (MQ=255)
                                                                                                                                           tATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGt    >  2:3322059/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTC  >  CP009273/364416‑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: