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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273344,8000GT88.3% 709.1 / 98.8 282E136D (GAG→GATprpB2‑methylisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base G (21/12);  new base T (110/139);  total (131/151)
Fisher's exact test for biased strand distribution p-value = 4.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.99e-01

Only 100 of 284 total aligned reads displayed.

TATCGGTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTT  >  CP009273/344653‑344949
                                                                                                                                                   |                                                                                                                                                     
tatCGGTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGatat                                                                                                                                                     <  2:1010265/150‑1 (MQ=255)
 atCGGTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATAtg                                                                                                                                                    >  2:5005884/1‑150 (MQ=255)
  tCGGTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGAtgg                                                                                                                                                   <  2:3462785/150‑1 (MQ=255)
  tCGGTTTTGGGTCTTCGGCCTTTACCGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATAtgg                                                                                                                                                   <  2:2157796/150‑1 (MQ=255)
    ggTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGAtggtg                                                                                                                                                 >  2:1403881/1‑150 (MQ=255)
        ttGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggatc                                                                                                                                             <  2:4916465/150‑1 (MQ=255)
         tGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggatcg                                                                                                                                            <  1:1090939/150‑1 (MQ=255)
          ggTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggatcgg                                                                                                                                           >  2:4751929/1‑150 (MQ=255)
           gTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggatcgga                                                                                                                                          <  1:2506735/150‑1 (MQ=255)
           gTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggatcgga                                                                                                                                          <  1:3296210/150‑1 (MQ=255)
             tcttcGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggatcggatc                                                                                                                                        >  1:132690/1‑150 (MQ=255)
              cttcGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATcc                                                                                                                                       >  1:1260594/1‑150 (MQ=255)
               ttcGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCcg                                                                                                                                      <  1:4857522/150‑1 (MQ=255)
                tcGGCCTTTAACGTGGCGCGCACCGTGAAATCTATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAATCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCcgc                                                                                                                                     >  1:4706385/1‑150 (MQ=255)
                tcGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCcgc                                                                                                                                     >  2:41660/1‑150 (MQ=255)
                    ccTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCgcgg                                                                                                                                 <  1:811040/150‑1 (MQ=255)
                     cTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCgcggc                                                                                                                                <  2:450155/150‑1 (MQ=255)
                       ttAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGcggcgg                                                                                                                              <  1:3815246/150‑1 (MQ=255)
                       ttAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGcggcgg                                                                                                                              >  1:3144169/1‑150 (MQ=255)
                        tAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGt                                                                                                                             >  1:5008413/1‑150 (MQ=255)
                        tAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGt                                                                                                                             <  1:2167353/150‑1 (MQ=255)
                         aaCGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTg                                                                                                                            <  2:468672/150‑1 (MQ=255)
                          aCGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTgg                                                                                                                           <  2:2479592/150‑1 (MQ=255)
                           cGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGACTAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGa                                                                                                                          >  2:3058302/1‑150 (MQ=255)
                              ggCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGc                                                                                                                       <  1:906715/150‑1 (MQ=255)
                                cgcgcACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGa                                                                                                                     <  2:3281168/150‑1 (MQ=255)
                                  cgcACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGaaa                                                                                                                   >  1:882592/1‑150 (MQ=255)
                                   gcACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGaaaa                                                                                                                  <  2:3577132/150‑1 (MQ=255)
                                       cGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGa                                                                                                              >  2:2247158/1‑150 (MQ=255)
                                       cGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGa                                                                                                              >  2:4792452/1‑150 (MQ=255)
                                         tGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATc                                                                                                            <  2:3734259/150‑1 (MQ=255)
                                          gAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATcc                                                                                                           >  2:602989/1‑150 (MQ=255)
                                           aaaTCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCt                                                                                                          >  1:2756191/1‑150 (MQ=255)
                                            aaTCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTg                                                                                                         >  1:577438/1‑150 (MQ=255)
                                             aTCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGa                                                                                                        <  1:1900870/150‑1 (MQ=255)
                                              tCGATGATTAAAACCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGAt                                                                                                       <  1:1628223/150‑1 (MQ=255)
                                                   gatTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTtgt                                                                                                  <  1:531639/150‑1 (MQ=255)
                                                       aaaGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATc                                                                                              >  2:458517/1‑150 (MQ=255)
                                                        aaGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCa                                                                                             >  2:4028340/1‑150 (MQ=255)
                                                         aGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCAt                                                                                            >  1:2657149/1‑150 (MQ=255)
                                                             tgagCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACTCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGCATCTGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGgcg                                                                                        >  2:2742405/4‑150 (MQ=255)
                                                             ggTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGgca                                                                                        <  2:157162/150‑2 (MQ=255)
                                                              gTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAAGAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAACACCGATCCTGATTTTGTGATCATGgcgc                                                                                       >  2:83007/1‑150 (MQ=255)
                                                                  ggCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCAc                                                                                   >  2:1268119/1‑150 (MQ=255)
                                                                    cAGGATTGCATATTGCAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCg                                                                                 <  2:339017/150‑1 (MQ=255)
                                                                    cAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCg                                                                                 >  2:4314494/1‑150 (MQ=255)
                                                                      ggATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGAt                                                                               <  2:3986900/150‑1 (MQ=255)
                                                                      ggATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGAt                                                                               >  1:4643253/1‑150 (MQ=255)
                                                                       gATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATg                                                                              >  2:3701389/1‑150 (MQ=255)
                                                                        cttGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGc                                                                             >  1:3970870/2‑150 (MQ=255)
                                                                        aTTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGACTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGc                                                                             >  1:3393505/1‑150 (MQ=255)
                                                                         ttGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGct                                                                            <  1:2702001/150‑1 (MQ=255)
                                                                         ttGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGct                                                                            <  1:2553574/150‑1 (MQ=255)
                                                                         ttGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGct                                                                            >  1:3761959/1‑150 (MQ=255)
                                                                         ttGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGct                                                                            <  1:1020939/150‑1 (MQ=255)
                                                                         ttGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGct                                                                            <  1:689110/150‑1 (MQ=255)
                                                                           gCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGctct                                                                          <  2:2379091/150‑1 (MQ=255)
                                                                               atTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCg                                                                      <  2:4420922/150‑1 (MQ=255)
                                                                                ttGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCgg                                                                     >  1:4322467/1‑150 (MQ=255)
                                                                                   aaGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGAGTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTag                                                                  >  1:348317/1‑150 (MQ=255)
                                                                                       tCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAggg                                                              >  2:3478553/1‑150 (MQ=255)
                                                                                        cAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAgggg                                                             >  2:126094/1‑150 (MQ=255)
                                                                                           gTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGCCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTg                                                          <  1:4094710/150‑1 (MQ=255)
                                                                                            ttGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTgg                                                         >  1:4798932/1‑150 (MQ=255)
                                                                                               gTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATg                                                      >  2:521965/1‑150 (MQ=255)
                                                                                                 gCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATgcg                                                    <  2:2402078/150‑1 (MQ=255)
                                                                                                   gAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATgcggc                                                  <  1:2793190/150‑1 (MQ=255)
                                                                                                       cGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATc                                              <  2:2915954/150‑1 (MQ=255)
                                                                                                       cGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATc                                              <  2:2846521/150‑1 (MQ=255)
                                                                                                       cGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATc                                              <  2:93562/150‑1 (MQ=255)
                                                                                                        gCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCg                                             >  2:3176077/1‑150 (MQ=255)
                                                                                                         cTGCGGTCATCGTCCGAATAACGCGATCGTCTCGACAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGa                                            <  1:3527880/150‑1 (MQ=255)
                                                                                                           gCGGTCATCGTCCGAATAAAGCGATCGTCGCGACAGAAGATCGGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGc                                          <  2:4037293/150‑1 (MQ=255)
                                                                                                            cGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCg                                         >  1:4894533/1‑150 (MQ=255)
                                                                                                             ggTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGt                                        <  1:1490965/150‑1 (MQ=255)
                                                                                                             ggTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGt                                        <  1:1491114/150‑1 (MQ=255)
                                                                                                               tCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTgc                                      <  2:2282580/150‑1 (MQ=255)
                                                                                                               tCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTgc                                      <  1:318320/150‑1 (MQ=255)
                                                                                                               tCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTgc                                      >  1:813203/1‑150 (MQ=255)
                                                                                                                  tCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGTTTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCa                                   <  1:13128/150‑1 (MQ=255)
                                                                                                                   cGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAg                                  >  1:204773/1‑150 (MQ=255)
                                                                                                                      ccGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGcc                               <  1:3676936/150‑1 (MQ=255)
                                                                                                                       cGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCt                              >  1:3541745/1‑150 (MQ=255)
                                                                                                                         aaTAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTAt                            >  1:3099034/1‑150 (MQ=255)
                                                                                                                          aTAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATg                           >  2:4219555/1‑150 (MQ=255)
                                                                                                                          aTAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATg                           <  1:4387879/150‑1 (MQ=255)
                                                                                                                             aaGCGATCGTCTCGACAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGCTCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTg                        <  2:1954659/150‑1 (MQ=255)
                                                                                                                             aaGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTg                        >  2:1787417/1‑150 (MQ=255)
                                                                                                                                cGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAg                     <  2:1430548/150‑1 (MQ=255)
                                                                                                                                 gATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGc                    >  1:3231974/1‑150 (MQ=255)
                                                                                                                                   tCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCgg                  >  1:1820635/1‑150 (MQ=255)
                                                                                                                                    cGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCggg                 >  2:853209/1‑150 (MQ=255)
                                                                                                                                    cGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCggg                 <  1:3257903/150‑1 (MQ=255)
                                                                                                                                     gTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGt                >  1:1448337/1‑150 (MQ=255)
                                                                                                                                       ctcGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGc              <  1:4453343/150‑1 (MQ=255)
                                                                                                                                       ctcGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGc              <  1:5080592/150‑1 (MQ=255)
                                                                                                                                       ctcGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGc              <  1:2896110/150‑1 (MQ=255)
                                                                                                                                         cGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCg            <  2:2672640/150‑1 (MQ=255)
                                                                                                                                          gAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCga           <  1:41745/150‑1 (MQ=255)
                                                                                                                                          gAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCga           >  1:1775327/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TATCGGTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTT  >  CP009273/344653‑344949

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: