Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 364,561 G→T 93.3% 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,5610GT93.3% 502.5 / 27.4 187D18Y (GAC→TAC) mhpB2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (10/3);  new base T (81/93);  total (91/96)
Fisher's exact test for biased strand distribution p-value = 4.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 187 total aligned reads displayed.

CCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTG  >  CP009273/364423‑364706
                                                                                                                                          |                                                                                                                                                 
 ccTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCa                                                                                                                                         >  2:377113/1‑148 (MQ=255)
     ggCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAagag                                                                                                                                     >  2:2716699/1‑148 (MQ=255)
        aaGACCCTGAATACACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTg                                                                                                                                  <  2:149734/148‑1 (MQ=255)
            cccTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCg                                                                                                                              <  1:2404290/148‑1 (MQ=255)
            cccTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCg                                                                                                                              <  2:198280/148‑1 (MQ=255)
             ccTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGa                                                                                                                             <  2:2372454/148‑1 (MQ=255)
             ccTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGa                                                                                                                             >  1:2833055/1‑148 (MQ=255)
                gAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGa                                                                                                                          >  2:2762270/1‑148 (MQ=255)
                gAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGa                                                                                                                          <  1:274538/148‑1 (MQ=255)
                gAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGa                                                                                                                          >  2:188419/1‑148 (MQ=255)
                  aTAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAgg                                                                                                                        >  1:1484899/1‑148 (MQ=255)
                  aTAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAgg                                                                                                                        >  1:3040249/1‑148 (MQ=255)
                  aTAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAgg                                                                                                                        >  2:1218509/1‑148 (MQ=255)
                  aTAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAgg                                                                                                                        >  2:1009174/1‑148 (MQ=255)
                   tAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGt                                                                                                                       <  1:1450248/148‑1 (MQ=255)
                   tAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGt                                                                                                                       <  1:2204595/148‑1 (MQ=255)
                    aaaCTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTc                                                                                                                      >  1:507209/1‑148 (MQ=255)
                      aCTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCaa                                                                                                                    <  1:1392611/148‑1 (MQ=255)
                      aCTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCaa                                                                                                                    >  1:1780957/1‑148 (MQ=255)
                       cTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAAt                                                                                                                   >  2:3427215/1‑148 (MQ=255)
                        tGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATg                                                                                                                  >  1:3437938/1‑148 (MQ=255)
                           cGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCg                                                                                                               >  1:3344989/1‑148 (MQ=255)
                             tCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTg                                                                                                             <  2:3484279/148‑1 (MQ=255)
                                gTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGAtt                                                                                                          >  1:1308700/1‑148 (MQ=255)
                                 tgatgaCGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTg                                                                                                         <  1:3492336/148‑1 (MQ=255)
                                  gatgaCGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGc                                                                                                        <  2:2216459/148‑1 (MQ=255)
                                  gatgaCGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGc                                                                                                        >  1:1022662/1‑148 (MQ=255)
                                   atgaCGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGcc                                                                                                       <  1:864571/148‑1 (MQ=255)
                                   atgaCGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGcc                                                                                                       >  2:623468/1‑148 (MQ=255)
                                         cTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGcc                                                                                                 <  1:1553998/148‑1 (MQ=255)
                                            aCCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCcgc                                                                                              >  1:2805653/1‑148 (MQ=255)
                                             cccGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCcgcg                                                                                             <  2:658120/148‑1 (MQ=255)
                                               cGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAg                                                                                           <  1:154347/148‑1 (MQ=255)
                                                   cTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTa                                                                                       >  1:1119120/1‑148 (MQ=255)
                                                   cTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTa                                                                                       >  2:447518/1‑148 (MQ=255)
                                                    tGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTAt                                                                                      <  2:365376/148‑1 (MQ=255)
                                                    tGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTAt                                                                                      <  2:2340751/148‑1 (MQ=255)
                                                     gATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTAtt                                                                                     <  1:2456786/148‑1 (MQ=255)
                                                      aTGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTg                                                                                    >  2:3478182/1‑148 (MQ=255)
                                                      aTGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTg                                                                                    >  1:2961833/1‑148 (MQ=255)
                                                        gCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCg                                                                                  >  2:247605/1‑148 (MQ=255)
                                                        gCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCg                                                                                  >  2:298230/1‑148 (MQ=255)
                                                         ccGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCgg                                                                                 <  1:766180/148‑1 (MQ=255)
                                                            aCGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCAt                                                                              >  1:3382888/1‑148 (MQ=255)
                                                              gTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATtc                                                                            <  2:701833/148‑1 (MQ=255)
                                                               tttCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATtct                                                                           >  2:2079844/1‑148 (MQ=255)
                                                                ttCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATtctc                                                                          >  2:1170638/1‑148 (MQ=255)
                                                                 tCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTcc                                                                         <  2:1309592/148‑1 (MQ=255)
                                                                 tCTGTCGAAAAGGTAGCCTGATATGCACGCTGATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTcc                                                                         <  2:2585202/148‑1 (MQ=255)
                                                                  cTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTccc                                                                        <  2:2697816/148‑1 (MQ=255)
                                                                  cTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTccc                                                                        <  2:3238296/148‑1 (MQ=255)
                                                                  cTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTccc                                                                        >  1:2128947/1‑148 (MQ=255)
                                                                   tGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTcccc                                                                       <  2:3319537/148‑1 (MQ=255)
                                                                       gAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGaa                                                                   >  1:2102791/1‑148 (MQ=255)
                                                                       gAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGaa                                                                   >  1:2287023/1‑148 (MQ=255)
                                                                           aGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtgg                                                               <  2:900262/148‑1 (MQ=255)
                                                                           aGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtgg                                                               >  2:2293425/1‑148 (MQ=255)
                                                                           aGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtgg                                                               >  2:966655/1‑148 (MQ=255)
                                                                            ggTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtggt                                                              >  1:3065983/1‑148 (MQ=255)
                                                                             gTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtggtg                                                             <  2:2460405/148‑1 (MQ=255)
                                                                              tAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtggtgg                                                            >  2:1047/1‑148 (MQ=255)
                                                                                   tGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTg                                                       >  1:658554/1‑148 (MQ=255)
                                                                                   tGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTg                                                       <  1:1848774/148‑1 (MQ=255)
                                                                                   tGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTg                                                       <  2:1780957/148‑1 (MQ=255)
                                                                                     atatGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGtt                                                     <  2:2412981/148‑1 (MQ=255)
                                                                                     atatGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGtt                                                     >  1:2403905/1‑148 (MQ=255)
                                                                                     atatGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGtt                                                     <  2:585538/148‑1 (MQ=255)
                                                                                       atGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTg                                                   >  1:3017212/1‑148 (MQ=255)
                                                                                         gCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTgcg                                                 >  1:1868738/1‑148 (MQ=255)
                                                                                              cTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGa                                            <  1:1395569/148‑1 (MQ=255)
                                                                                                 aTCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCa                                         <  1:1719432/148‑1 (MQ=255)
                                                                                                       aCTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACaa                                   <  1:2578866/148‑1 (MQ=255)
                                                                                                       aCTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACaa                                   >  2:1070064/1‑148 (MQ=255)
                                                                                                       aCTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACaa                                   <  2:542641/148‑1 (MQ=255)
                                                                                                        cTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCCTTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAAc                                  >  1:860737/1‑148 (MQ=255)
                                                                                                        cTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAAc                                  <  2:3537112/148‑1 (MQ=255)
                                                                                                         tGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACg                                 >  2:1114236/1‑148 (MQ=255)
                                                                                                          gTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACgg                                >  2:3437527/1‑148 (MQ=255)
                                                                                                          gTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACgg                                <  1:3564242/148‑1 (MQ=255)
                                                                                                          gTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACgg                                <  1:2945007/148‑1 (MQ=255)
                                                                                                          gTCTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACgg                                <  2:3209651/148‑1 (MQ=255)
                                                                                                            cTTTCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCt                              >  1:1095048/1‑148 (MQ=255)
                                                                                                              ttCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCttt                            >  1:951698/1‑148 (MQ=255)
                                                                                                              ttCCCACTCGCCGCTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCttt                            <  1:208207/148‑1 (MQ=255)
                                                                                                                      cgccgcTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGa                    <  2:1556555/148‑1 (MQ=255)
                                                                                                                        ccgcTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACg                  <  1:1152732/148‑1 (MQ=255)
                                                                                                                          gcTGGTGGGGTATGTCTACCCGGCGCTAGAGGTGCTAGATGAGGTCAATGGCGTGATTGACAGTGCCCGCGAGCGTATTGCGGCTTTCTCCCCTTAACTGGTGGTGCTGTTTGCGCAAGTTCACTACAACGGCTTTTTCTATGAAGTg                <  2:1065246/148‑1 (MQ=255)
                                                                                                                          gcTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTg                <  1:1210022/148‑1 (MQ=255)
                                                                                                                           cTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGa               <  2:1696423/148‑1 (MQ=255)
                                                                                                                           cTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGa               <  2:1696246/148‑1 (MQ=255)
                                                                                                                           cTGGTGGGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGa               <  2:712370/148‑1 (MQ=255)
                                                                                                                               tggGGTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGTTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACTGCTTTTTCTATGACGTGATGcc           >  1:1918293/1‑148 (MQ=255)
                                                                                                                                ggggTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCa          >  2:1315650/1‑148 (MQ=255)
                                                                                                                                 gggTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCAc         >  1:961277/1‑148 (MQ=255)
                                                                                                                                 gggTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCAc         >  1:338756/1‑148 (MQ=255)
                                                                                                                                 gggTATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCAc         >  1:2647658/1‑148 (MQ=255)
                                                                                                                                    tATGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGt      >  2:2530693/1‑148 (MQ=255)
                                                                                                                                     aTGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGtt     <  1:2762270/148‑1 (MQ=255)
                                                                                                                                     aTGTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGtt     <  1:188419/148‑1 (MQ=255)
                                                                                                                                       gTCTACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCt   >  2:97218/1‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
CCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTG  >  CP009273/364423‑364706

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: