Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 757,353 A→T 11.5% N126Y (AAT→TAT)  sucB → dihydrolipoyltranssuccinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273757,3530AT11.5% 566.4 / 53.8 192N126Y (AAT→TAT) sucBdihydrolipoyltranssuccinase
Reads supporting (aligned to +/- strand):  ref base A (99/71);  new base T (16/6);  total (115/77)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 198 total aligned reads displayed.

CTGCGTGAAGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTA  >  CP009273/757206‑757500
                                                                                                                                                   |                                                                                                                                                   
cTGCGTGAAGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAAcaca                                                                                                                                                     <  2:2962033/148‑1 (MQ=255)
 tGCGTGAAGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACaa                                                                                                                                                    <  2:813278/148‑1 (MQ=255)
      gAAGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCg                                                                                                                                               <  2:3088223/148‑1 (MQ=255)
       aaGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACTATCTCGa                                                                                                                                              >  2:409436/1‑148 (MQ=255)
       aaGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGa                                                                                                                                              <  1:2168995/148‑1 (MQ=255)
         ggCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACg                                                                                                                                            >  1:886294/1‑148 (MQ=255)
          gCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGc                                                                                                                                           >  1:2704506/1‑148 (MQ=255)
                gcgcCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAgcgc                                                                                                                                     <  2:1395372/148‑1 (MQ=255)
                 cgcCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACTATCTCGACGCCAGCGcc                                                                                                                                    >  1:3267893/1‑148 (MQ=255)
                          aGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGc                                                                                                                           >  2:2068116/1‑148 (MQ=255)
                          aGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAAAACAATCTCGACGCCAGCGCCATTAAAGGc                                                                                                                           >  1:1899380/1‑148 (MQ=255)
                            aaaCCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCAc                                                                                                                         <  2:332214/148‑1 (MQ=255)
                              aCCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCg                                                                                                                       <  2:2036318/148‑1 (MQ=255)
                               ccAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCgg                                                                                                                      <  2:2514954/148‑1 (MQ=255)
                                  gcgcCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGgtgt                                                                                                                   >  2:472736/1‑148 (MQ=255)
                                   cgcCAAATCTGAAGACAACGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACTATCTCGACGCCAGCGCCATTAAAGGCACCGGtgtg                                                                                                                  <  1:852166/148‑1 (MQ=255)
                                    gcCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTgtgg                                                                                                                 <  1:1965375/148‑1 (MQ=255)
                                     ccAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTgtggg                                                                                                                >  1:337012/1‑148 (MQ=255)
                                     ccAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTgtggg                                                                                                                >  1:1292479/1‑148 (MQ=255)
                                       aaaTCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTgtgggtg                                                                                                              <  2:564922/148‑1 (MQ=255)
                                            tGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGgtcgt                                                                                                         <  1:3179914/148‑1 (MQ=255)
                                             gAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGgtcgtc                                                                                                        <  2:260523/148‑1 (MQ=255)
                                              aaGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCt                                                                                                       <  2:2497714/148‑1 (MQ=255)
                                               agagAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATGAAAGGCACCGGTGTGGGTGGTCGTCTg                                                                                                      >  2:2294994/1‑148 (MQ=255)
                                                   aaaGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTc                                                                                                  >  1:1931474/1‑148 (MQ=255)
                                                   aaaGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTc                                                                                                  >  1:1102367/1‑148 (MQ=255)
                                                     aGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACTATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGt                                                                                                >  2:571234/1‑148 (MQ=255)
                                                      gCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACTATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTg                                                                                               <  2:1690217/148‑1 (MQ=255)
                                                         tCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACTATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAg                                                                                            >  2:2865947/1‑148 (MQ=255)
                                                         tCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAg                                                                                            <  2:1806350/148‑1 (MQ=255)
                                                         tCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAg                                                                                            <  1:3209662/148‑1 (MQ=255)
                                                          ccACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGa                                                                                           >  1:2803152/1‑148 (MQ=255)
                                                          ccACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCAGTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGa                                                                                           >  1:2519548/1‑148 (MQ=255)
                                                                  gcgcAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGaaa                                                                                   >  1:3001912/1‑148 (MQ=255)
                                                                  gcgcAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGaaa                                                                                   >  1:3000854/1‑148 (MQ=255)
                                                                   cgcAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGaaaa                                                                                  >  2:2846333/1‑148 (MQ=255)
                                                                   cgcAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGaaaa                                                                                  <  2:3253395/148‑1 (MQ=255)
                                                                    gcAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGaaaaa                                                                                 <  2:1755583/148‑1 (MQ=255)
                                                                    gcAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGaaaaa                                                                                 <  1:1148024/148‑1 (MQ=255)
                                                                     cAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAAc                                                                                <  1:820792/148‑1 (MQ=255)
                                                                      aaCGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACa                                                                               <  2:534512/148‑1 (MQ=255)
                                                                      aaCGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACa                                                                               <  1:511250/148‑1 (MQ=255)
                                                                       aCGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACAt                                                                              >  1:753470/1‑148 (MQ=255)
                                                                        cGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACTATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATc                                                                             >  2:65779/1‑148 (MQ=255)
                                                                        cGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATc                                                                             >  1:1821020/1‑148 (MQ=255)
                                                                        cGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATc                                                                             <  2:2169674/148‑1 (MQ=255)
                                                                         gCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCt                                                                            <  1:1350672/148‑1 (MQ=255)
                                                                         gCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCt                                                                            <  1:1194007/148‑1 (MQ=255)
                                                                           cagcagGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTgg                                                                          <  2:1899380/148‑1 (MQ=255)
                                                                              cagGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACTATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGa                                                                       >  2:3013511/1‑148 (MQ=255)
                                                                               agGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGaa                                                                      <  2:2349207/148‑1 (MQ=255)
                                                                               agGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGaa                                                                      <  2:1064293/148‑1 (MQ=255)
                                                                                 gCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAg                                                                    >  2:3026949/1‑148 (MQ=255)
                                                                                 gCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAg                                                                    >  2:855843/1‑148 (MQ=255)
                                                                                   gTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGcc                                                                  <  1:642959/148‑1 (MQ=255)
                                                                                   gTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGcc                                                                  <  2:1764147/148‑1 (MQ=255)
                                                                                    tctcTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGccc                                                                 >  1:3102920/1‑148 (MQ=255)
                                                                                    tctcTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGccc                                                                 >  2:2942385/1‑148 (MQ=255)
                                                                                    tctcTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGccc                                                                 >  1:2127865/1‑148 (MQ=255)
                                                                                     ctctGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACTATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGcccc                                                                >  2:1851605/1‑148 (MQ=255)
                                                                                       ctGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCgg                                                              >  2:759613/1‑148 (MQ=255)
                                                                                         ggAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGGGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCg                                                            >  2:373706/1‑148 (MQ=255)
                                                                                          gAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGa                                                           >  1:1085418/1‑148 (MQ=255)
                                                                                               gCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTCAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAagag                                                      >  1:628847/1‑148 (MQ=255)
                                                                                                  aaaCAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCt                                                   >  2:1866432/1‑148 (MQ=255)
                                                                                                  aaaCAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCt                                                   >  1:2296245/1‑148 (MQ=255)
                                                                                                    acaacGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGc                                                 >  1:1541260/1‑148 (MQ=255)
                                                                                                     caacGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCt                                                >  1:536487/1‑148 (MQ=255)
                                                                                                     caacGATGCGTTAAGCCCGGCGATCCGTAGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCt                                                >  2:614260/1‑148 (MQ=255)
                                                                                                      aacGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACTATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTc                                               <  1:3013515/148‑1 (MQ=255)
                                                                                                        cGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACTATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCg                                             >  2:859027/1‑148 (MQ=255)
                                                                                                          aTGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGc                                           <  1:2464801/148‑1 (MQ=255)
                                                                                                            gCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCCCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAg                                         >  1:464907/1‑148 (MQ=255)
                                                                                                                tAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGgc                                     >  1:2553728/1‑148 (MQ=255)
                                                                                                                  aGCCCGGCGATCCGTCGCCTGCTGGCTGAACACACTCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGgctg                                   <  1:2074494/148‑1 (MQ=255)
                                                                                                                  aGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGgctg                                   >  2:1285629/1‑148 (MQ=255)
                                                                                                                  aGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGgctg                                   >  1:1790247/1‑148 (MQ=255)
                                                                                                                   gCCCGGCGATCCGTCGCCTGCTGGCTGAACACTATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGgctgc                                  >  2:2705821/1‑148 (MQ=255)
                                                                                                                   gCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTCAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGgctgc                                  >  2:2489499/1‑148 (MQ=255)
                                                                                                                     ccGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTc                                <  2:2801594/148‑1 (MQ=255)
                                                                                                                        gCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCcgg                             >  2:3226063/1‑148 (MQ=255)
                                                                                                                         cGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCcggc                            >  2:3292408/1‑148 (MQ=255)
                                                                                                                           aTCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCcggcgg                          >  1:1547462/1‑148 (MQ=255)
                                                                                                                            tCCGTCGCCTGCTGGCTGAACACTATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCggcggc                         >  2:273987/1‑148 (MQ=255)
                                                                                                                            tCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCggcggc                         >  2:1974312/1‑148 (MQ=255)
                                                                                                                              cGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGgcgc                       >  2:824528/1‑148 (MQ=255)
                                                                                                                                tCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCaa                     >  1:1635466/1‑148 (MQ=255)
                                                                                                                                tCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCaa                     <  1:2362689/148‑1 (MQ=255)
                                                                                                                                   ccTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCg                  >  2:314786/1‑148 (MQ=255)
                                                                                                                                   ccTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCg                  >  2:1960396/1‑148 (MQ=255)
                                                                                                                                    ctgctgGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCgg                 <  2:2673522/148‑1 (MQ=255)
                                                                                                                                      gctggctgAACACAATCTCGACGCCCGCGCCTTTATAGGAACGGGTGTGGTTGTTCTTGTGACTCTTGAAGATGTGGAAAAACTTCTGGCGAAAGCCCCGGCGAAATTGTGTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGct               <  1:731335/148‑1 (MQ=255)
                                                                                                                                      gctggctgAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGGGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGct               <  2:2089334/148‑1 (MQ=255)
                                                                                                                                       ctggctgAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGctc              >  2:731335/1‑148 (MQ=255)
                                                                                                                                       ctggctgAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGct               >  2:2393932/1‑147 (MQ=255)
                                                                                                                                       ctggctgAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGct               <  1:2393932/147‑1 (MQ=255)
                                                                                                                                               acacAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCa      >  1:2638377/1‑148 (MQ=255)
                                                                                                                                               acacAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCa      >  1:2638400/1‑148 (MQ=255)
                                                                                                                                                   aaTCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTa  >  2:1789482/1‑148 (MQ=255)
                                                                                                                                                   aaTCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTa  >  2:116439/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTGCGTGAAGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTA  >  CP009273/757206‑757500

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: