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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273344,8000GT89.1% 629.7 / 72.1 246E136D (GAG→GAT) prpB2‑methylisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base G (5/22);  new base T (100/119);  total (105/141)
Fisher's exact test for biased strand distribution p-value = 7.17e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 250 total aligned reads displayed.

TCGGTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGTT  >  CP009273/344655‑344946
                                                                                                                                                 |                                                                                                                                                  
tCGGTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATAtgg                                                                                                                                                <  2:2048294/150‑1 (MQ=255)
  ggTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATAtggtg                                                                                                                                              >  1:3233329/1‑150 (MQ=255)
   gTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATAtggtgg                                                                                                                                             <  1:3192210/150‑1 (MQ=255)
    ttttGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTgga                                                                                                                                            <  2:38743/150‑1 (MQ=255)
    ttttGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTgga                                                                                                                                            <  2:3024027/150‑1 (MQ=255)
     tttGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggat                                                                                                                                           <  1:1187890/150‑1 (MQ=255)
       tGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggatcg                                                                                                                                         >  1:47607/1‑150 (MQ=255)
        ggTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggatcgg                                                                                                                                        <  2:1096805/150‑1 (MQ=255)
        ggTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTggatcgg                                                                                                                                        <  2:2833815/150‑1 (MQ=255)
          ttcttcGGCCTTGAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTggatcggat                                                                                                                                      <  2:1389187/150‑1 (MQ=255)
                 gCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCcgcgcg                                                                                                                               <  2:1456426/150‑1 (MQ=255)
                  ccTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCgcgg                                                                                                                              >  2:4356081/1‑150 (MQ=255)
                  ccTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCgcgg                                                                                                                              >  1:1072661/1‑150 (MQ=255)
                   cTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCgcggc                                                                                                                             >  2:8966/1‑150 (MQ=255)
                      tAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGt                                                                                                                          <  2:539632/150‑1 (MQ=255)
                       aaCGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTg                                                                                                                         >  1:4245791/1‑150 (MQ=255)
                       aaCGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTg                                                                                                                         >  2:3206118/1‑150 (MQ=255)
                        aCGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGAGTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCTATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTgg                                                                                                                        <  2:2759164/150‑1 (MQ=255)
                        aCGTGGAGCGCACTGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTTCATATTGATTATCAGGTTGGTGCGAAACTCTGCGGTCATCGTCCGAATAAAGCTATCGTATCGAATGTAGATATGGTGGATCGGATCCGCGCGGCtgtgt                                                                                                                        <  2:1003302/150‑2 (MQ=255)
                         cGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGa                                                                                                                       >  1:580765/1‑150 (MQ=255)
                          gTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGAt                                                                                                                      <  1:2193424/150‑1 (MQ=255)
                            ggCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGc                                                                                                                    <  1:350409/150‑1 (MQ=255)
                             gcgcgcACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCg                                                                                                                   <  1:1647079/150‑1 (MQ=255)
                                cgcACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGaaa                                                                                                                <  2:3547128/150‑1 (MQ=255)
                                  cACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAAc                                                                                                              >  1:1686776/1‑150 (MQ=255)
                                    ccGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCg                                                                                                            <  1:87397/150‑1 (MQ=255)
                                    ccGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCg                                                                                                            <  1:2390668/150‑1 (MQ=255)
                                    ccGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCg                                                                                                            >  1:4259139/1‑150 (MQ=255)
                                        gAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATcc                                                                                                        <  1:993995/150‑1 (MQ=255)
                                          aaTCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCTCGGCGGTGGATGCGAAAACCGATCCTg                                                                                                      >  2:3796795/1‑150 (MQ=255)
                                          aaTCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTg                                                                                                      >  2:521656/1‑150 (MQ=255)
                                             cGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGAtt                                                                                                   >  2:1259280/1‑150 (MQ=255)
                                             cGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGAtt                                                                                                   >  2:4054578/1‑150 (MQ=255)
                                                  atTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTtgtg                                                                                              >  1:3756893/1‑150 (MQ=255)
                                                  atTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTtgtg                                                                                              >  1:1709356/1‑150 (MQ=255)
                                                  atTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTtgtg                                                                                              <  2:4305447/150‑1 (MQ=255)
                                                  atTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTtgtg                                                                                              <  2:4305168/150‑1 (MQ=255)
                                                    tAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGAt                                                                                            <  2:1838981/150‑1 (MQ=255)
                                                    tAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGAt                                                                                            <  2:2289287/150‑1 (MQ=255)
                                                    tAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGAt                                                                                            <  2:3533887/150‑1 (MQ=255)
                                                    tAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGAt                                                                                            <  2:1953735/150‑1 (MQ=255)
                                                     aaaGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATc                                                                                           >  2:4367397/1‑150 (MQ=255)
                                                     aaaGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATc                                                                                           >  2:2124580/1‑150 (MQ=255)
                                                           ggTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGgcg                                                                                     >  1:3247414/1‑150 (MQ=255)
                                                             tGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGgcgcg                                                                                   <  2:3250195/150‑1 (MQ=255)
                                                             tGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGgcgcg                                                                                   <  1:1308675/150‑1 (MQ=255)
                                                                  cAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCg                                                                              <  2:3261021/150‑1 (MQ=255)
                                                                  cAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCg                                                                              <  1:2187343/150‑1 (MQ=255)
                                                                   aGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGa                                                                             >  2:622912/1‑150 (MQ=255)
                                                                    ggATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGAt                                                                            <  1:1535545/150‑1 (MQ=255)
                                                                     gATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATg                                                                           <  1:1141049/150‑1 (MQ=255)
                                                                       ttGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGct                                                                         <  1:2162694/150‑1 (MQ=255)
                                                                        tGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGctc                                                                        >  1:2705806/1‑150 (MQ=255)
                                                                         gCAGATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGctct                                                                       <  2:3616101/150‑1 (MQ=255)
                                                                            tatTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGc                                                                    >  2:3031625/1‑150 (MQ=255)
                                                                            tatTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGc                                                                    >  2:753724/1‑150 (MQ=255)
                                                                              ttGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCgg                                                                  >  2:1686534/1‑150 (MQ=255)
                                                                               tGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGt                                                                 >  1:4014050/1‑150 (MQ=255)
                                                                                 aaGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTag                                                               <  2:1023530/150‑1 (MQ=255)
                                                                                  aGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTaga                                                              >  2:1932814/1‑150 (MQ=255)
                                                                                   gATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTagag                                                             >  1:2281478/1‑150 (MQ=255)
                                                                                    aTCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAgg                                                            >  2:3328422/1‑150 (MQ=255)
                                                                                    aTCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAgg                                                            >  2:2448934/1‑150 (MQ=255)
                                                                                     tCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAggg                                                           <  1:1746595/150‑1 (MQ=255)
                                                                                        ggTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCt                                                        <  1:1558912/150‑1 (MQ=255)
                                                                                              tGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATgc                                                  >  1:1193854/1‑150 (MQ=255)
                                                                                                 gAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATgcggc                                               >  1:3155736/1‑150 (MQ=255)
                                                                                                   aaCGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGa                                             <  2:398160/150‑1 (MQ=255)
                                                                                                       cTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGa                                         <  2:2018547/150‑1 (MQ=255)
                                                                                                                tCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCa                                >  1:4525617/1‑150 (MQ=255)
                                                                                                                tCGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCa                                >  1:1981189/1‑150 (MQ=255)
                                                                                                                 cGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAg                               <  1:597434/150‑1 (MQ=255)
                                                                                                                 cGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAg                               <  2:2254984/150‑1 (MQ=255)
                                                                                                                 cGTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAg                               >  1:4343308/1‑150 (MQ=255)
                                                                                                                  gTCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAgg                              >  1:392484/1‑150 (MQ=255)
                                                                                                                   tCCGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGc                             >  1:4090557/1‑150 (MQ=255)
                                                                                                                     cGAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCt                           <  1:1794223/150‑1 (MQ=255)
                                                                                                                      gAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGTTCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGAGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTa                          <  1:2312055/150‑1 (MQ=255)
                                                                                                                      gAATAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTa                          <  1:2589877/150‑1 (MQ=255)
                                                                                                                       aaTAAAGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTAt                         >  1:1419166/1‑150 (MQ=255)
                                                                                                                          aaaGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGtt                      <  2:3113097/150‑1 (MQ=255)
                                                                                                                          aaaGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGtt                      >  2:3224832/1‑150 (MQ=255)
                                                                                                                          aaaGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGtt                      <  1:3635853/150‑1 (MQ=255)
                                                                                                                          aaaGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGtt                      <  2:828958/150‑1 (MQ=255)
                                                                                                                            aGCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGa                    >  1:2951130/1‑150 (MQ=255)
                                                                                                                            aGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGa                    >  1:3458752/1‑150 (MQ=255)
                                                                                                                             gCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGaa                   <  2:1036161/150‑1 (MQ=255)
                                                                                                                             gCGATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGaa                   <  2:1056929/150‑1 (MQ=255)
                                                                                                                               gATCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGc                 >  2:1644785/1‑150 (MQ=255)
                                                                                                                                 tCGTCTCGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCgg               <  1:3085858/150‑1 (MQ=255)
                                                                                                                                    tctcGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTg            >  2:356930/1‑150 (MQ=255)
                                                                                                                                      tcGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGcc          <  2:1030344/150‑1 (MQ=255)
                                                                                                                                       cGAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCg         <  2:2578121/150‑1 (MQ=255)
                                                                                                                                        gAAAGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCga        <  1:4356241/150‑1 (MQ=255)
                                                                                                                                         aaaGAAGATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCgag       <  1:3618833/150‑1 (MQ=255)
                                                                                                                                          aagaagATATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCgaga      <  1:20836/150‑1 (MQ=255)
                                                                                                                                           agaagaTATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGAt     >  1:4266510/1‑150 (MQ=255)
                                                                                                                                              agatATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGAtgtt  >  1:1095293/1‑150 (MQ=255)
                                                                                                                                              agatATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGAtgtt  >  1:924302/1‑150 (MQ=255)
                                                                                                                                              agatATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGAtgtt  >  1:3376234/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TCGGTTTTGGTTCTTCGGCCTTTAACGTGGCGCGCACCGTGAAATCGATGATTAAAGCCGGTGCGGCAGGATTGCATATTGAAGATCAGGTTGGTGCGAAACGCTGCGGTCATCGTCCGAATAAAGCGATCGTCTCGAAAGAAGAGATGGTGGATCGGATCCGCGCGGCGGTGGATGCGAAAACCGATCCTGATTTTGTGATCATGGCGCGCACCGATGCTCTGGCGGTAGAGGGGCTGGATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGTT  >  CP009273/344655‑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: ‑