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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273344,8000GT95.2% 782.8 / 29.4 253E136D (GAG→GAT) prpB2‑methylisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base G (9/3);  new base T (131/110);  total (140/113)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 256 total aligned reads displayed.

TATCGGTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGT  >  CP009273/344653‑344945
                                                                                                                                                   |                                                                                                                                                 
tatCGGTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGatat                                                                                                                                                 >  2:4533897/1‑150 (MQ=255)
tatCGGTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGatat                                                                                                                                                 <  2:1043415/150‑1 (MQ=255)
 atCGGTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATAtg                                                                                                                                                <  1:2624960/150‑1 (MQ=255)
        ttGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggatc                                                                                                                                         <  1:3297575/150‑1 (MQ=255)
         tGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggatcg                                                                                                                                        <  2:502206/150‑1 (MQ=255)
           gTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggatcgga                                                                                                                                      >  2:1240815/1‑150 (MQ=255)
               ttcGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCcg                                                                                                                                  <  1:2497417/150‑1 (MQ=255)
                tcGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCcgc                                                                                                                                 <  1:863434/150‑1 (MQ=255)
                 cGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCcgcg                                                                                                                                <  2:2342854/150‑1 (MQ=255)
                  ggCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCcgcgc                                                                                                                               >  2:1991744/1‑150 (MQ=255)
                    ccTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCgcgg                                                                                                                             <  1:73433/150‑1 (MQ=255)
                       ttAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGcggcgg                                                                                                                          >  1:4028772/1‑150 (MQ=255)
                       ttAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGcggcgg                                                                                                                          <  1:581522/150‑1 (MQ=255)
                          aCGTGGCGCGCACCGAGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTgg                                                                                                                       <  2:3595167/150‑1 (MQ=255)
                           cGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGa                                                                                                                      >  2:112140/1‑150 (MQ=255)
                           cGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGa                                                                                                                      <  2:260926/150‑1 (MQ=255)
                            gTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGAt                                                                                                                     <  2:425882/150‑1 (MQ=255)
                              ggCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGc                                                                                                                   <  1:3120839/150‑1 (MQ=255)
                               gcgcgcACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCg                                                                                                                  <  2:3107471/150‑1 (MQ=255)
                                cgcgcACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGa                                                                                                                 >  2:735782/1‑150 (MQ=255)
                                cgcgcACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGa                                                                                                                 <  1:2155892/150‑1 (MQ=255)
                                   gcACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGaaaa                                                                                                              >  2:1422212/1‑150 (MQ=255)
                                   gcACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGaaaa                                                                                                              >  1:3574008/1‑150 (MQ=255)
                                    cACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAAc                                                                                                             <  2:2936532/150‑1 (MQ=255)
                                     aCCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAAcc                                                                                                            <  2:2703090/150‑1 (MQ=255)
                                      ccGTGAAATCGATGATTACAGCCGGTGCGGCAGGATGGCATATTGAAGATCAGGGTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCg                                                                                                           <  2:2809157/150‑1 (MQ=255)
                                       cGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGa                                                                                                          >  2:4165705/1‑150 (MQ=255)
                                        gTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGAt                                                                                                         >  1:2149597/1‑150 (MQ=255)
                                          gAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATcc                                                                                                       <  2:1380560/150‑1 (MQ=255)
                                           aaaTCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCt                                                                                                      <  2:3815237/150‑1 (MQ=255)
                                             aTCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGa                                                                                                    >  1:4625460/1‑150 (MQ=255)
                                               cGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGAtt                                                                                                  >  1:4163536/1‑150 (MQ=255)
                                               cGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGAtt                                                                                                  >  2:4068009/1‑150 (MQ=255)
                                                gatgatTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGAttt                                                                                                 >  1:2637441/1‑150 (MQ=255)
                                                gatgatTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGAttt                                                                                                 >  1:2259913/1‑150 (MQ=255)
                                                 atgatTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGAtttt                                                                                                >  2:1736561/1‑150 (MQ=255)
                                                   gatTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTtgt                                                                                              <  1:224377/150‑1 (MQ=255)
                                                        aaGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCa                                                                                         >  2:4302656/1‑150 (MQ=255)
                                                          gCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACTCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCTAAAGAAGATATGGTGTATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATg                                                                                       >  1:2241288/1‑150 (MQ=255)
                                                            cGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGgc                                                                                     >  2:3611577/1‑150 (MQ=255)
                                                            cGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGgc                                                                                     <  1:2446904/150‑1 (MQ=255)
                                                             ggTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGgcg                                                                                    >  1:1547512/1‑150 (MQ=255)
                                                             ggTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGgcg                                                                                    >  1:4442117/1‑150 (MQ=255)
                                                                 cGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCa                                                                                >  2:677702/1‑150 (MQ=255)
                                                                 cGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCa                                                                                <  1:3681738/150‑1 (MQ=255)
                                                                 cGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCa                                                                                >  2:1749540/1‑150 (MQ=255)
                                                                    cAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCg                                                                             <  1:113580/150‑1 (MQ=255)
                                                                     aGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGa                                                                            <  2:2112079/150‑1 (MQ=255)
                                                                       gATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATg                                                                          >  1:1845337/1‑150 (MQ=255)
                                                                       gATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATg                                                                          <  1:3544494/150‑1 (MQ=255)
                                                                        aTTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGc                                                                         <  1:3224535/150‑1 (MQ=255)
                                                                           gCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGctct                                                                      <  1:3002727/150‑1 (MQ=255)
                                                                           gCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGctct                                                                      <  2:3738491/150‑1 (MQ=255)
                                                                            cATATTTAAGAACAGGTTGGTGCGAAACGCTGCGGTCATCGACCGAATAAAGCGATCGTCTCGAAAGAATATATGGTGGATCGGATCCGAGCGGCGGTCGATGCGAAAACCGATCCTGATTTAGAGATCATGGCGCGCACCGATActcag                                                                     >  2:4465942/1‑148 (MQ=255)
                                                                            cATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTg                                                                     >  1:2067805/1‑150 (MQ=255)
                                                                            cATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTg                                                                     >  1:3967118/1‑150 (MQ=255)
                                                                              tatTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGc                                                                   >  2:2173350/1‑150 (MQ=255)
                                                                               atTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCg                                                                  >  1:4579663/1‑150 (MQ=255)
                                                                                 tGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGt                                                                >  2:849631/1‑150 (MQ=255)
                                                                                 tGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCAATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGt                                                                <  2:1276106/150‑1 (MQ=255)
                                                                                  gAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTa                                                               >  2:394109/1‑150 (MQ=255)
                                                                                  gAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTa                                                               <  1:3050240/150‑1 (MQ=255)
                                                                                  gAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTa                                                               >  2:4327004/1‑150 (MQ=255)
                                                                                  gAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTa                                                               >  2:4393811/1‑150 (MQ=255)
                                                                                    aGATCAGGTTGGTGCGAAACGCTGCGGTCANCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTaga                                                             >  1:4285170/1‑150 (MQ=255)
                                                                                      aTCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAgg                                                           <  2:3971083/150‑1 (MQ=255)
                                                                                      aTCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAgg                                                           >  1:3766883/1‑150 (MQ=255)
                                                                                       tCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAggg                                                          <  1:4309569/150‑1 (MQ=255)
                                                                                          ggTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCt                                                       >  1:3371720/1‑150 (MQ=255)
                                                                                             tGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGa                                                    <  2:1648923/150‑1 (MQ=255)
                                                                                                 gCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATgcg                                                <  1:3505142/150‑1 (MQ=255)
                                                                                                  cGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATgcgg                                               <  2:1033664/150‑1 (MQ=255)
                                                                                                  cGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATgcgg                                               >  2:2604983/1‑150 (MQ=255)
                                                                                                   gAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATgcggc                                              <  2:3293276/150‑1 (MQ=255)
                                                                                                      aCGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGAt                                           >  1:3434377/1‑150 (MQ=255)
                                                                                                          tGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAg                                       >  2:1285313/1‑150 (MQ=255)
                                                                                                           gCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGc                                      >  2:2769839/1‑150 (MQ=255)
                                                                                                             ggTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGt                                    <  1:2014076/150‑1 (MQ=255)
                                                                                                             ggTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGt                                    >  1:3668364/1‑150 (MQ=255)
                                                                                                             ggTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGt                                    >  2:887884/1‑150 (MQ=255)
                                                                                                              gTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTg                                   >  1:3621863/1‑150 (MQ=255)
                                                                                                               tCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTgc                                  <  2:870586/150‑1 (MQ=255)
                                                                                                                cATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTgcg                                 >  2:1332008/1‑150 (MQ=255)
                                                                                                                  tCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCa                               <  1:2937554/150‑1 (MQ=255)
                                                                                                                   cGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAg                              >  2:4626957/1‑150 (MQ=255)
                                                                                                                    gTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAgg                             >  2:3703746/1‑150 (MQ=255)
                                                                                                                      ccGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGcc                           <  2:1806829/150‑1 (MQ=255)
                                                                                                                       cGAATAAAGCGATCGTCTCGAAAGAATATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTAGTGATCATGGCGAGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGAGCGCAGGCCt                          >  1:2426619/1‑150 (MQ=255)
                                                                                                                        gAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTa                         <  2:4586864/150‑1 (MQ=255)
                                                                                                                         aaTAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTAt                        <  1:1091637/150‑1 (MQ=255)
                                                                                                                         aaTAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTAt                        >  2:3634389/1‑150 (MQ=255)
                                                                                                                             aaGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCAGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGAGGCGATCGAGCGTGCGCAGGCCTATGTTg                    >  1:3405631/1‑150 (MQ=255)
                                                                                                                               gCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGaa                  >  2:2618854/1‑150 (MQ=255)
                                                                                                                                cGATCGTCTCGAAAGAAGAGATGGTTGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGTCGCGCACCGATGCTCTGGCTGTAGAGTGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGTAg                 <  1:2791042/150‑1 (MQ=255)
                                                                                                                                  aTCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCg               <  2:752159/150‑1 (MQ=255)
                                                                                                                                      tctcGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTg           >  1:2555199/1‑150 (MQ=255)
                                                                                                                                           aaaGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGTGCAGGCCTATGTTGAAGCGGGTGCCgag      <  2:2375589/150‑1 (MQ=255)
                                                                                                                                            aagaagATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCgaga     >  2:3491746/1‑150 (MQ=255)
                                                                                                                                              gaagaTATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGAtg   >  1:4187176/1‑150 (MQ=255)
                                                                                                                                               aagaTATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGAtgt  >  2:1334500/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TATCGGTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGT  >  CP009273/344653‑344945

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