Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 344,800 G→T 93.0% E136D (GAG→GAT)  prpB → 2‑methylisocitrate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273344,8000GT93.0% 790.8 / 47.2 268E136D (GAG→GAT) prpB2‑methylisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base G (7/12);  new base T (131/118);  total (138/130)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

Only 100 of 273 total aligned reads displayed.

ATATCGGTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGTT  >  CP009273/344652‑344946
                                                                                                                                                    |                                                                                                                                                  
atatCGGTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAgaga                                                                                                                                                   >  1:1413108/1‑150 (MQ=255)
 tatCGGTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGatat                                                                                                                                                  >  1:2934892/1‑150 (MQ=255)
    cGGTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATAtggt                                                                                                                                               <  2:1174122/150‑1 (MQ=255)
      gTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATAtggtgg                                                                                                                                             <  2:2521610/150‑1 (MQ=255)
       ttttGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTgga                                                                                                                                            >  1:1708360/1‑150 (MQ=255)
        tttGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggat                                                                                                                                           <  2:1136309/150‑1 (MQ=255)
         ttGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggatc                                                                                                                                          <  2:3293498/150‑1 (MQ=255)
          tGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggatcg                                                                                                                                         <  1:2562195/150‑1 (MQ=255)
          tGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggatcg                                                                                                                                         <  2:465717/150‑1 (MQ=255)
            gTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggatcgga                                                                                                                                       >  1:3200173/1‑150 (MQ=255)
             ttcttcGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTggatcggat                                                                                                                                      >  1:3715664/1‑150 (MQ=255)
              tcttcGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggatcggatc                                                                                                                                     <  2:3987878/150‑1 (MQ=255)
              tcttcGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggatcggatc                                                                                                                                     <  2:2803408/150‑1 (MQ=255)
               cttcGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATcc                                                                                                                                    <  2:1615522/150‑1 (MQ=255)
                  cGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCcgcg                                                                                                                                 <  1:3121554/150‑1 (MQ=255)
                   ggCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCcgcgc                                                                                                                                >  2:634811/1‑150 (MQ=255)
                   ggCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCcgcgc                                                                                                                                <  1:3705507/150‑1 (MQ=255)
                   ggCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCcgcgc                                                                                                                                >  1:536761/1‑150 (MQ=255)
                   ggCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCcgcgc                                                                                                                                >  2:3745012/1‑150 (MQ=255)
                     ccTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCgcgg                                                                                                                              >  1:1195059/1‑150 (MQ=255)
                       tttAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCgcggcg                                                                                                                            >  2:3026854/1‑150 (MQ=255)
                         tAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGt                                                                                                                          >  2:4220754/1‑150 (MQ=255)
                          aaCGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTg                                                                                                                         <  2:4448793/150‑1 (MQ=255)
                           aCGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTgg                                                                                                                        <  1:411380/150‑1 (MQ=255)
                             gTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGAt                                                                                                                      <  2:667996/150‑1 (MQ=255)
                              tGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGCAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATg                                                                                                                     <  1:46457/150‑1 (MQ=255)
                                 cgcgcACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGa                                                                                                                  >  2:2883712/1‑150 (MQ=255)
                                 agcgcACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGa                                                                                                                  >  2:1704639/2‑150 (MQ=255)
                                   cgcACCGTGAAATCTATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGaaa                                                                                                                <  2:2325245/150‑1 (MQ=255)
                                   cgcACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGaaa                                                                                                                <  2:2993679/150‑1 (MQ=255)
                                    gcACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGaaaa                                                                                                               >  1:3630768/1‑150 (MQ=255)
                                    gcACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGaaaa                                                                                                               >  1:1217609/1‑150 (MQ=255)
                                     cACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAAc                                                                                                              >  2:2705724/1‑150 (MQ=255)
                                       ccGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCg                                                                                                            >  2:1802985/1‑150 (MQ=255)
                                        cGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGa                                                                                                           <  2:3430319/150‑1 (MQ=255)
                                          tGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATc                                                                                                         >  2:4784172/1‑150 (MQ=255)
                                           gAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATcc                                                                                                        >  2:2126899/1‑150 (MQ=255)
                                            aaaTCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCt                                                                                                       <  1:1623592/150‑1 (MQ=255)
                                               tCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATACAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGAt                                                                                                    <  1:1891657/150‑1 (MQ=255)
                                                  atgatTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGAtttt                                                                                                 <  1:785232/150‑1 (MQ=255)
                                                  atgatTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGAtttt                                                                                                 <  2:3943305/150‑1 (MQ=255)
                                                    gatTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTtgt                                                                                               <  1:3926344/150‑1 (MQ=255)
                                                     atTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTtgtg                                                                                              <  2:3394772/150‑1 (MQ=255)
                                                      ttAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGa                                                                                             <  1:3469977/150‑1 (MQ=255)
                                                        acaGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATc                                                                                           <  1:657203/148‑1 (MQ=255)
                                                        aaaGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATc                                                                                           >  2:3455896/1‑150 (MQ=255)
                                                        aaaGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATc                                                                                           >  1:1270222/1‑150 (MQ=255)
                                                         aaGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCa                                                                                          >  1:272050/1‑150 (MQ=255)
                                                          aGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCAt                                                                                         <  1:327489/150‑1 (MQ=255)
                                                             cGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGgc                                                                                      >  1:486803/1‑150 (MQ=255)
                                                             cGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGgc                                                                                      <  2:1597918/150‑1 (MQ=255)
                                                              ggTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGgcg                                                                                     >  1:4761920/1‑150 (MQ=255)
                                                                tGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGgcgcg                                                                                   >  2:814886/1‑150 (MQ=255)
                                                                 gCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGgcgcgc                                                                                  <  1:4165664/150‑1 (MQ=255)
                                                                    gCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCAcc                                                                               <  1:1100675/150‑1 (MQ=255)
                                                                     cAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTTATTTTGTGATCATGGCGCGCACCg                                                                              >  2:4528876/1‑150 (MQ=255)
                                                                      aGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGa                                                                             >  2:456749/1‑150 (MQ=255)
                                                                       ggATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGAt                                                                            <  1:2381952/150‑1 (MQ=255)
                                                                       ggATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGAt                                                                            <  1:1601591/150‑1 (MQ=255)
                                                                             cATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTg                                                                      <  1:3020123/150‑1 (MQ=255)
                                                                              atatTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTgg                                                                     <  2:3080737/150‑1 (MQ=255)
                                                                                  tGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGt                                                                 <  2:1566466/150‑1 (MQ=255)
                                                                                  tGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGt                                                                 <  1:2630720/150‑1 (MQ=255)
                                                                                     aGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTaga                                                              >  1:860019/1‑150 (MQ=255)
                                                                                       aTCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAgg                                                            >  1:2137838/1‑150 (MQ=255)
                                                                                         cAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATTAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAgggg                                                          <  2:2204841/150‑1 (MQ=255)
                                                                                          aGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGc                                                         >  2:875024/1‑150 (MQ=255)
                                                                                           ggTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCt                                                        >  2:578678/1‑150 (MQ=255)
                                                                                           ggTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCt                                                        >  2:2969782/1‑150 (MQ=255)
                                                                                               ggTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATAGGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGAt                                                    >  1:2097668/1‑150 (MQ=255)
                                                                                                gTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGTGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATg                                                   <  1:4496837/150‑1 (MQ=255)
                                                                                                 tGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATgc                                                  <  2:2057959/150‑1 (MQ=255)
                                                                                                    gAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATgcggc                                               >  2:4491156/1‑150 (MQ=255)
                                                                                                    gAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATgcggc                                               >  1:1908916/1‑150 (MQ=255)
                                                                                                    gAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATgcggc                                               >  1:3790747/1‑150 (MQ=255)
                                                                                                    gAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATgcggc                                               >  1:4631592/1‑150 (MQ=255)
                                                                                                     aaaCGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATGTTGTGATCATGGCGCGCACCGAAGCTCTGGCGGTCGAGGGGCTGGATgcggcg                                              >  2:350746/1‑150 (MQ=255)
                                                                                                       aCGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGAt                                            >  1:998644/1‑150 (MQ=255)
                                                                                                        cGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATc                                           <  1:4309586/150‑1 (MQ=255)
                                                                                                                  aTCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTgcgc                                 >  2:1126709/1‑150 (MQ=255)
                                                                                                                        cGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCt                           <  2:3518675/150‑1 (MQ=255)
                                                                                                                         gAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTa                          <  2:195796/150‑1 (MQ=255)
                                                                                                                          aagaaaGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTAt                         <  2:1692954/147‑1 (MQ=255)
                                                                                                                          aaTAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTAt                         >  1:2819872/1‑150 (MQ=255)
                                                                                                                           aTAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATg                        <  2:1709352/150‑1 (MQ=255)
                                                                                                                           aTAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATg                        <  2:1916234/150‑1 (MQ=255)
                                                                                                                             aaaGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGtt                      <  1:4505687/150‑1 (MQ=255)
                                                                                                                              aaGCGATCGTCTCGAAAGAAGATATGGTGGTTCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTg                     <  2:101165/150‑1 (MQ=255)
                                                                                                                                 cGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAg                  >  2:1397191/1‑150 (MQ=255)
                                                                                                                                  gATCGTCTCGAAAGAAGATAAGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTATGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGc                 >  1:2962196/1‑150 (MQ=255)
                                                                                                                                  gATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGc                 >  1:4683379/1‑150 (MQ=255)
                                                                                                                                   aTCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCg                >  2:4652393/1‑150 (MQ=255)
                                                                                                                                      gTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGt             >  1:1145974/1‑150 (MQ=255)
                                                                                                                                       tctcGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTg            <  1:1005563/150‑1 (MQ=255)
                                                                                                                                          cGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCg         <  2:2444178/150‑1 (MQ=255)
                                                                                                                                           gAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCga        >  1:2201810/1‑150 (MQ=255)
                                                                                                                                           gAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCga        <  1:1768675/150‑1 (MQ=255)
                                                                                                                                            aaaGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCgag       >  1:3001879/1‑150 (MQ=255)
                                                                                                                                                aagaTATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGAtgt   <  2:1226791/150‑1 (MQ=255)
                                                                                                                                                 agatATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGAtgtt  <  2:3588594/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATATCGGTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGTT  >  CP009273/344652‑344946

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: ‑