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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273757,3530AT 5.1% 555.2 / 15.8 175N126Y (AAT→TAT) sucBdihydrolipoyltranssuccinase
Reads supporting (aligned to +/- strand):  ref base A (92/74);  new base T (6/3);  total (98/77)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 183 total aligned reads displayed.

TGCGTGAAGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGT  >  CP009273/757207‑757499
                                                                                                                                                  |                                                                                                                                                  
tGCGTGAAGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACaa                                                                                                                                                   <  1:1196659/148‑1 (MQ=255)
tGCGTGAAGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACaa                                                                                                                                                   <  1:255278/148‑1 (MQ=255)
  cGTGAAGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAAtc                                                                                                                                                 <  1:1200324/148‑1 (MQ=255)
     gAAGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCg                                                                                                                                              <  1:638981/148‑1 (MQ=255)
         gCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCTCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGc                                                                                                                                          >  1:262457/1‑148 (MQ=255)
         gCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGc                                                                                                                                          <  2:562164/148‑1 (MQ=255)
          cAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGcc                                                                                                                                         >  1:3051112/1‑148 (MQ=255)
            aCAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAg                                                                                                                                       <  2:1871883/148‑1 (MQ=255)
            aCAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAg                                                                                                                                       <  2:1871866/148‑1 (MQ=255)
             cAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAgc                                                                                                                                      >  2:215529/1‑148 (MQ=255)
                cgcCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGcc                                                                                                                                   <  1:1150573/148‑1 (MQ=255)
                cgcCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGcc                                                                                                                                   >  1:3456729/1‑148 (MQ=255)
                cgcCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGcc                                                                                                                                   <  1:2923732/148‑1 (MQ=255)
                    ggTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTa                                                                                                                               <  1:1659775/148‑1 (MQ=255)
                      tAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCAGTaaa                                                                                                                             >  1:349211/1‑148 (MQ=255)
                       aaaGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAg                                                                                                                            <  2:170763/148‑1 (MQ=255)
                        aaGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAgg                                                                                                                           >  2:938968/1‑148 (MQ=255)
                         aGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACTATCTCGACGCCAGCGCCATTAAAGGc                                                                                                                          >  1:3147912/1‑148 (MQ=255)
                         aGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGc                                                                                                                          >  2:221256/1‑148 (MQ=255)
                         aGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGc                                                                                                                          >  1:1363654/1‑148 (MQ=255)
                             aCCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCg                                                                                                                      >  2:651959/1‑148 (MQ=255)
                              ccAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCgg                                                                                                                     >  1:529937/1‑148 (MQ=255)
                                     cAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGgggtgggt                                                                                                              >  2:802932/1‑148 (MQ=255)
                                     cAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTgtgggt                                                                                                              >  2:802937/1‑148 (MQ=255)
                                      aaaTCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTgtgggtg                                                                                                             >  2:726406/1‑148 (MQ=255)
                                       aaTCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTACAGGCACCGGTgtgggtgg                                                                                                            >  2:283572/1‑148 (MQ=255)
                                           tGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGgtcgt                                                                                                        <  2:902620/148‑1 (MQ=255)
                                           tGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGgtcgt                                                                                                        <  2:2681484/148‑1 (MQ=255)
                                           tGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGgtcgt                                                                                                        <  1:2148621/148‑1 (MQ=255)
                                            gAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGgtcgtc                                                                                                       >  2:3054747/1‑148 (MQ=255)
                                            gAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGggcgtc                                                                                                       >  2:264834/1‑148 (MQ=255)
                                             aaGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCt                                                                                                      >  1:2958058/1‑148 (MQ=255)
                                                agaAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGAc                                                                                                   <  2:2510672/148‑1 (MQ=255)
                                                agaAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGAc                                                                                                   >  2:2318633/1‑148 (MQ=255)
                                                   aaGCGTCCACTCCGGCGCAACGCCCGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCg                                                                                                <  1:2198367/148‑1 (MQ=255)
                                                     gCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTg                                                                                              <  2:3051112/148‑1 (MQ=255)
                                                      cGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGa                                                                                             <  1:189689/148‑1 (MQ=255)
                                                      cGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGa                                                                                             >  2:3030451/1‑148 (MQ=255)
                                                        tCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAg                                                                                           >  2:1823037/1‑148 (MQ=255)
                                                         ccACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGa                                                                                          >  1:2518841/1‑148 (MQ=255)
                                                         ccACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGa                                                                                          >  1:2518704/1‑148 (MQ=255)
                                                         ccACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGa                                                                                          >  1:527418/1‑148 (MQ=255)
                                                          cACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGAt                                                                                         <  2:1407816/148‑1 (MQ=255)
                                                             tCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGAtgtg                                                                                      >  1:682782/1‑148 (MQ=255)
                                                              ccGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTgg                                                                                     <  1:2973418/148‑1 (MQ=255)
                                                                  cgcAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGaaaa                                                                                 <  1:1444959/148‑1 (MQ=255)
                                                                      aCGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACAt                                                                             >  2:2049111/1‑148 (MQ=255)
                                                                       cGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATc                                                                            >  1:2468792/1‑148 (MQ=255)
                                                                       cGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATc                                                                            <  2:3372981/148‑1 (MQ=255)
                                                                         ccAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAACAACATCTg                                                                          >  1:1776895/1‑148 (MQ=255)
                                                                           agcagGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGc                                                                        <  1:2241931/148‑1 (MQ=255)
                                                                             cagGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGa                                                                      >  1:6809/1‑148 (MQ=255)
                                                                                gCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAg                                                                   >  2:877452/1‑148 (MQ=255)
                                                                                 cGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGc                                                                  >  2:3411470/1‑148 (MQ=255)
                                                                                  gTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGcc                                                                 <  1:2497125/148‑1 (MQ=255)
                                                                                   tctcTGGAAGAGCAAAACACCGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGccc                                                                <  2:3254033/148‑1 (MQ=255)
                                                                                   tctcTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGccc                                                                <  1:181237/148‑1 (MQ=255)
                                                                                   tctcTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGccc                                                                <  2:162479/148‑1 (MQ=255)
                                                                                    ctctGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGcccc                                                               >  1:2716143/1‑148 (MQ=255)
                                                                                      ctGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCCTTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCgg                                                             >  2:2010037/1‑148 (MQ=255)
                                                                                      ctGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCgg                                                             >  1:2784167/1‑148 (MQ=255)
                                                                                        ggAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCg                                                           >  1:1606790/1‑148 (MQ=255)
                                                                                          aaGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGaa                                                         <  2:1642752/148‑1 (MQ=255)
                                                                                          aaGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGaa                                                         >  2:1145099/1‑148 (MQ=255)
                                                                                          aaGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGaa                                                         >  2:3080242/1‑148 (MQ=255)
                                                                                           agagCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGaaa                                                        >  1:2011808/1‑148 (MQ=255)
                                                                                            gagCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGCGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAag                                                       <  2:1680885/148‑1 (MQ=255)
                                                                                              gCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAagag                                                     <  2:3062019/148‑1 (MQ=255)
                                                                                               cAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGt                                                    <  1:1155248/148‑1 (MQ=255)
                                                                                               cAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGt                                                    <  1:3221977/148‑1 (MQ=255)
                                                                                                aaaaCAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTc                                                   <  2:3121251/148‑1 (MQ=255)
                                                                                                 aaaCAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCt                                                  >  1:919532/1‑148 (MQ=255)
                                                                                                 aaaCAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCt                                                  >  1:3488798/1‑148 (MQ=255)
                                                                                                   acaacGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGc                                                >  2:2855304/1‑148 (MQ=255)
                                                                                                   acaacGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCAGTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGc                                                >  2:1549251/1‑148 (MQ=255)
                                                                                                      acGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTcc                                             >  1:2378013/1‑148 (MQ=255)
                                                                                                           gCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAg                                        >  1:2875194/1‑148 (MQ=255)
                                                                                                               tAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGgc                                    >  2:2332571/1‑148 (MQ=255)
                                                                                                                 aGCCCGGCGATCCGTCGCCTGCTGGCTGAACACTATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGgctg                                  >  1:888839/1‑148 (MQ=255)
                                                                                                                   cccGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGgctgct                                >  1:1710402/1‑148 (MQ=255)
                                                                                                                   cccGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGgctgct                                <  2:529730/148‑1 (MQ=255)
                                                                                                                    ccGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTc                               <  1:1794813/148‑1 (MQ=255)
                                                                                                                           tCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTAGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCggcggc                        <  1:2030720/148‑1 (MQ=255)
                                                                                                                             cGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGgcgc                      >  2:2805288/1‑148 (MQ=255)
                                                                                                                             cGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGgcgc                      >  2:2944264/1‑148 (MQ=255)
                                                                                                                              gTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCa                     <  1:802941/148‑1 (MQ=255)
                                                                                                                              gTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCa                     <  1:802932/148‑1 (MQ=255)
                                                                                                                              gTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCa                     >  1:3069086/1‑148 (MQ=255)
                                                                                                                               tCGCCTGCTGGCTGAACACTATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCaa                    <  2:1720788/148‑1 (MQ=255)
                                                                                                                               tCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCaa                    <  1:936707/148‑1 (MQ=255)
                                                                                                                                    tgctgGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGc               >  1:1449371/1‑148 (MQ=255)
                                                                                                                                    tgctgGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGc               >  2:954423/1‑148 (MQ=255)
                                                                                                                                      ctggctgAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGctc             >  1:1833814/1‑148 (MQ=255)
                                                                                                                                        ggctgAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTg           >  1:2766818/1‑148 (MQ=255)
                                                                                                                                        ggctgAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTg           >  1:2669287/1‑148 (MQ=255)
                                                                                                                                            gAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTg       >  1:438490/1‑148 (MQ=255)
                                                                                                                                             aaCACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGc      >  2:2303985/1‑148 (MQ=255)
                                                                                                                                              acacAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCa     >  1:3041169/1‑148 (MQ=255)
                                                                                                                                                acaATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACg   >  2:1252654/1‑148 (MQ=255)
                                                                                                                                                 caATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGt  <  2:349211/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGCGTGAAGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGT  >  CP009273/757207‑757499

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: