Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,247,740 G→T 7.7% G407C (GGC→TGC)  dhaR → dhaKLM operon transcription activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,247,7400GT 7.7% 906.6 / 36.8 244G407C (GGC→TGC) dhaRdhaKLM operon transcription activator
Reads supporting (aligned to +/- strand):  ref base G (104/121);  new base T (12/7);  total (116/128)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 254 total aligned reads displayed.

AGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGATGCGCGGCGTTTAATACCAATT  >  CP009273/1247591‑1247889
                                                                                                                                                     |                                                                                                                                                     
agGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATg                                                                                                                                                       >  1:4253339/1‑150 (MQ=255)
agGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATg                                                                                                                                                       <  1:3213886/150‑1 (MQ=255)
agGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATg                                                                                                                                                       >  1:1356049/1‑150 (MQ=255)
      tGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTc                                                                                                                                                 >  1:4354353/1‑150 (MQ=255)
           aaGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGt                                                                                                                                            >  2:2949313/1‑150 (MQ=255)
             ggCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCg                                                                                                                                          <  1:2661204/150‑1 (MQ=255)
                   gCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGa                                                                                                                                    >  2:1810169/1‑150 (MQ=255)
                       aGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTg                                                                                                                                <  1:3059301/150‑1 (MQ=255)
                         tCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGc                                                                                                                              <  1:4453175/150‑1 (MQ=255)
                          cAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGca                                                                                                                             <  2:523863/150‑1 (MQ=255)
                          cAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGca                                                                                                                             <  2:2027039/150‑1 (MQ=255)
                           aGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGGATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGcac                                                                                                                            <  2:802428/150‑1 (MQ=255)
                            ggCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGcaca                                                                                                                           <  1:3685030/150‑1 (MQ=255)
                                aTTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACAcggc                                                                                                                       >  2:506315/1‑150 (MQ=255)
                                 ttCCTAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACAcggcg                                                                                                                      <  2:4487105/150‑1 (MQ=255)
                                 ttCATAATGAAAGCTAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACAcggcg                                                                                                                      >  1:4826141/1‑150 (MQ=255)
                                 ttCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGGGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACAcggcg                                                                                                                      >  2:2350137/1‑150 (MQ=255)
                                    aTAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCa                                                                                                                   >  1:2105533/1‑150 (MQ=255)
                                     tAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCAc                                                                                                                  <  2:2643018/150‑1 (MQ=255)
                                      aaTGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACg                                                                                                                 <  1:8317/150‑1 (MQ=255)
                                      aaTGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACg                                                                                                                 <  2:963948/150‑1 (MQ=255)
                                        tGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCt                                                                                                               <  2:60959/150‑1 (MQ=255)
                                           aaGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGtt                                                                                                            <  1:1125009/150‑1 (MQ=255)
                                              cgagcgTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCt                                                                                                         <  2:1000331/150‑1 (MQ=255)
                                              cgagcgTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCt                                                                                                         <  1:45183/150‑1 (MQ=255)
                                               gagcgTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCtt                                                                                                        <  2:3127433/150‑1 (MQ=255)
                                                  cgTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGaa                                                                                                     <  2:4711624/150‑1 (MQ=255)
                                                     gctgcAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAg                                                                                                  >  1:2873460/1‑150 (MQ=255)
                                                      ctgcAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGa                                                                                                 >  2:2191660/1‑150 (MQ=255)
                                                        gcAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGAtt                                                                                               <  2:1626225/150‑1 (MQ=255)
                                                        gcAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGAtt                                                                                               <  2:866423/150‑1 (MQ=255)
                                                            gTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAat                                                                                           <  2:3083258/150‑1 (MQ=255)
                                                             tCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAata                                                                                          >  1:556121/1‑150 (MQ=255)
                                                              ccTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAatat                                                                                         <  2:2745852/150‑1 (MQ=255)
                                                              ccTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAatat                                                                                         <  2:4170419/150‑1 (MQ=255)
                                                              ccTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAatat                                                                                         <  1:1299890/150‑1 (MQ=255)
                                                               cTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATc                                                                                        >  1:332980/1‑150 (MQ=255)
                                                                ttATATCGCCGTCAATTGTGAGTTATATGGTGATGCTTCGCTGGCGGAAGAATGTATTGGTGGCGATCGCACTGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTAGAAAAGAGTGAATATCt                                                                                       >  2:4729031/1‑150 (MQ=255)
                                                                ttATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCt                                                                                       >  1:3571027/1‑150 (MQ=255)
                                                                 tataTCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTg                                                                                      <  2:1404837/150‑1 (MQ=255)
                                                                  atatCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTgg                                                                                     >  1:3681410/1‑150 (MQ=255)
                                                                   tatCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGc                                                                                    <  1:2067606/150‑1 (MQ=255)
                                                                    atCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCg                                                                                   >  2:1648173/1‑150 (MQ=255)
                                                                       gCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATTGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTg                                                                                <  2:3907177/150‑1 (MQ=255)
                                                                        ccGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTgg                                                                               <  2:3317692/150‑1 (MQ=255)
                                                                         cGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATTGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGACAAGATTGAATCTCTGGCGGTGGa                                                                              >  1:2324342/1‑150 (MQ=255)
                                                                           tCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGt                                                                            <  1:3409217/150‑1 (MQ=255)
                                                                           tCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGt                                                                            >  2:3326283/1‑150 (MQ=255)
                                                                           tCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGt                                                                            >  2:3791313/1‑150 (MQ=255)
                                                                            cAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGtt                                                                           >  2:1107605/1‑150 (MQ=255)
                                                                             aaTTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTa                                                                          <  2:3859001/150‑1 (MQ=255)
                                                                              aTTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTAc                                                                         <  1:925200/150‑1 (MQ=255)
                                                                               ttGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACa                                                                        >  1:3052272/1‑150 (MQ=255)
                                                                                 gtgAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGt                                                                      <  2:371268/150‑1 (MQ=255)
                                                                                  tgAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTc                                                                     >  1:4206445/1‑150 (MQ=255)
                                                                                    aGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCtg                                                                   >  1:1184695/1‑150 (MQ=255)
                                                                                    aGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCtg                                                                   >  2:3231055/1‑150 (MQ=255)
                                                                                      ttATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCtgct                                                                 <  2:1742818/150‑1 (MQ=255)
                                                                                        atatGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCtgcttt                                                               <  2:1192179/150‑1 (MQ=255)
                                                                                           tGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCtgctttgct                                                            >  2:3844075/1‑150 (MQ=255)
                                                                                             gTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTc                                                          <  1:3334407/150‑1 (MQ=255)
                                                                                               gATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAg                                                        <  1:1567515/150‑1 (MQ=255)
                                                                                                 tgctgcGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATTGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGt                                                      >  1:808586/1‑150 (MQ=255)
                                                                                                  gctgcGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGtt                                                     >  2:3437305/1‑150 (MQ=255)
                                                                                                   ctgcGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTa                                                    >  2:2381288/1‑150 (MQ=255)
                                                                                                    tgcGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTAt                                                   <  2:1335432/150‑1 (MQ=255)
                                                                                                         tGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGc                                              >  1:4092181/1‑150 (MQ=255)
                                                                                                           gCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATTGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAg                                            >  2:1972739/1‑150 (MQ=255)
                                                                                                             ggAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAggg                                          >  2:3049954/1‑150 (MQ=255)
                                                                                                               aagaaTTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAggggg                                        >  2:4413878/1‑150 (MQ=255)
                                                                                                                 gaaTTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGtt                                      >  2:3985565/1‑150 (MQ=255)
                                                                                                                 gaaTTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGtt                                      >  1:4583469/1‑150 (MQ=255)
                                                                                                                  aaTTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTa                                     <  2:2453903/150‑1 (MQ=255)
                                                                                                                   aTTTATTGGTGGCGATCGCACGGACAATGAAAATTGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTAt                                    >  1:3960697/1‑150 (MQ=255)
                                                                                                                    tttATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATc                                   >  1:3667739/1‑150 (MQ=255)
                                                                                                                     ttATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCa                                  <  2:2813067/150‑1 (MQ=255)
                                                                                                                            tggCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACt                           <  1:3014377/150‑1 (MQ=255)
                                                                                                                             ggCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTg                          <  2:2593380/150‑1 (MQ=255)
                                                                                                                             ggCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTg                          <  1:1038716/150‑1 (MQ=255)
                                                                                                                                gATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGAt                       >  2:3882578/1‑150 (MQ=255)
                                                                                                                                gATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGAt                       <  1:3905459/150‑1 (MQ=255)
                                                                                                                                gATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGAt                       <  1:904448/150‑1 (MQ=255)
                                                                                                                                gATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGAt                       <  1:1666226/150‑1 (MQ=255)
                                                                                                                                gATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGAt                       >  1:3117901/1‑150 (MQ=255)
                                                                                                                                 aTCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGATg                      <  1:4635511/150‑1 (MQ=255)
                                                                                                                                 aTCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGATg                      <  1:1800117/150‑1 (MQ=255)
                                                                                                                                 aTCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGATg                      <  1:2748790/150‑1 (MQ=255)
                                                                                                                                     cACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGATGcgcg                  <  1:1478521/150‑1 (MQ=255)
                                                                                                                                     cACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGATGcgcg                  >  2:4310812/1‑150 (MQ=255)
                                                                                                                                     cACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGATGcgcg                  >  2:3236297/1‑150 (MQ=255)
                                                                                                                                          aCAATGAAAATTGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGATGCGCGGCGtt             >  2:3546699/1‑150 (MQ=255)
                                                                                                                                          aCAATGAAAAGGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAAGATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGATGCGCGGCGtt             <  2:4059979/150‑1 (MQ=255)
                                                                                                                                           cAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGATGCGCGGCGttt            >  1:1007580/1‑150 (MQ=255)
                                                                                                                                             aTGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGATGCGCGGCGTTTaa          <  1:4631861/150‑1 (MQ=255)
                                                                                                                                             aTGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGATGCGCGGCGTTTaa          >  2:4365988/1‑150 (MQ=255)
                                                                                                                                                 aaaTGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGATGCGCGGCGTTTAATAcc      >  2:3385504/1‑150 (MQ=255)
                                                                                                                                                  aaTGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGATGCGCGGCGTTTAATACCa     <  2:4763456/150‑1 (MQ=255)
                                                                                                                                                   aTGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGATGCGCGGCGTTTAATACCaa    <  1:2541328/150‑1 (MQ=255)
                                                                                                                                                    tGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGATGCGCGGCGTTTAATACCAAt   <  2:1922236/150‑1 (MQ=255)
                                                                                                                                                     ggCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGATGCGCGGCGTTTAATACCAAtt  <  1:2018576/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGAACTGGCACACGGCGGCACGCTGTTTCTTGAAAAGATTGAATATCTGGCGGTGGAGTTACAGTCTGCTTTGCTTCAGGTTATCAAGCAGGGGGTTATCACGCGACTGGATGCGCGGCGTTTAATACCAATT  >  CP009273/1247591‑1247889

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: