Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 754,572 G→T 91.0% E137D (GAG→GAT sucA → 2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273754,5720GT91.0% 699.7 / 59.5 252E137D (GAG→GATsucA2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base G (7/16);  new base T (115/114);  total (122/130)
Fisher's exact test for biased strand distribution p-value = 8.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 252 total aligned reads displayed.

AGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACA  >  CP009273/754427‑754718
                                                                                                                                                 |                                                                                                                                                  
aGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCt                                                                                                                                                <  1:357267/150‑1 (MQ=255)
  gTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTc                                                                                                                                              <  2:3870952/150‑1 (MQ=255)
   ttAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCa                                                                                                                                             <  1:3016566/150‑1 (MQ=255)
    tAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCaa                                                                                                                                            <  2:315673/150‑1 (MQ=255)
     aaaGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAAc                                                                                                                                           <  1:477069/150‑1 (MQ=255)
      aaGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACg                                                                                                                                          >  1:3917112/1‑150 (MQ=255)
            tGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGtt                                                                                                                                    <  2:629273/150‑1 (MQ=255)
            tGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGtt                                                                                                                                    <  2:1334693/150‑1 (MQ=255)
            tGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGtt                                                                                                                                    <  1:1535589/150‑1 (MQ=255)
            tGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGtt                                                                                                                                    >  2:2239406/1‑150 (MQ=255)
             gCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTc                                                                                                                                   <  2:677403/150‑1 (MQ=255)
             gCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTc                                                                                                                                   >  1:4333730/1‑150 (MQ=255)
              cAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTcc                                                                                                                                  >  2:962859/1‑149 (MQ=255)
              cAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCa                                                                                                                                  >  2:2507699/1‑150 (MQ=255)
               aGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCAt                                                                                                                                 >  1:1625845/1‑150 (MQ=255)
                    aTTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGcc                                                                                                                            <  1:995979/150‑1 (MQ=255)
                    aTTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGcc                                                                                                                            >  1:808842/1‑150 (MQ=255)
                     ttAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCa                                                                                                                           >  1:2117703/1‑150 (MQ=255)
                      tAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAg                                                                                                                          >  1:97830/1‑150 (MQ=255)
                       aaCGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGc                                                                                                                         <  2:813270/150‑1 (MQ=255)
                       aaCGCATACCGCGTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGc                                                                                                                         <  2:813230/150‑1 (MQ=255)
                          gCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGc                                                                                                                      >  2:3980013/1‑150 (MQ=255)
                             tACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCaaa                                                                                                                   <  1:841198/150‑1 (MQ=255)
                               ccGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGa                                                                                                                 <  1:826431/150‑1 (MQ=255)
                               ccGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGa                                                                                                                 <  1:103861/150‑1 (MQ=255)
                               ccGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGa                                                                                                                 <  1:2407004/150‑1 (MQ=255)
                                cGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGAGAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGaa                                                                                                                >  2:4537360/1‑150 (MQ=255)
                                       gTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGa                                                                                                         >  1:1177868/1‑150 (MQ=255)
                                        tGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGaa                                                                                                        <  1:2825607/150‑1 (MQ=255)
                                        tGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGaa                                                                                                        <  2:539059/150‑1 (MQ=255)
                                         ggTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGaaa                                                                                                       <  2:952130/150‑1 (MQ=255)
                                          gTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGACAGCGGCAAAGAAACCATGAAAc                                                                                                      >  2:1222487/1‑150 (MQ=255)
                                           tCACCCGCATCCGAATCTCGATCCGCTGGGACTGTGGCAGCCAGATACCGTCGCCGATCTGGATCCGTCTTTCCACGATCTGACCTAAGCAGACTTCCAGGATCCCTTCAACGTCGTTTCATGTGCCAGCGGCAAAGAAACCATGAAACt                                                                                                     <  2:2688668/150‑1 (MQ=255)
                                           tCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACt                                                                                                     <  1:2504306/150‑1 (MQ=255)
                                              ccAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCgg                                                                                                  <  1:2244202/150‑1 (MQ=255)
                                                 gCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGa                                                                                               <  2:3334060/150‑1 (MQ=255)
                                                 gCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGa                                                                                               <  2:3610311/150‑1 (MQ=255)
                                                    tGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAgct                                                                                            <  1:3854313/150‑1 (MQ=255)
                                                      cGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAgctgc                                                                                          <  1:2006403/150‑1 (MQ=255)
                                                      cGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAgctgc                                                                                          >  1:491122/1‑150 (MQ=255)
                                                       gAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAgctgct                                                                                         <  1:166131/150‑1 (MQ=255)
                                                          tctcGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGa                                                                                      <  1:448932/150‑1 (MQ=255)
                                                          tctcGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGa                                                                                      <  1:838395/150‑1 (MQ=255)
                                                          tctcGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGa                                                                                      <  1:838238/150‑1 (MQ=255)
                                                           ctcGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGaa                                                                                     <  1:706560/150‑1 (MQ=255)
                                                            tcGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAg                                                                                    >  2:2998024/1‑150 (MQ=255)
                                                            tcGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAg                                                                                    >  2:2046223/1‑150 (MQ=255)
                                                            tcGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTAGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGATGCAGTAAg                                                                                    >  1:3226724/1‑150 (MQ=255)
                                                              gctCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTggtatcc                                                                                  <  2:1973430/148‑6 (MQ=255)
                                                              gATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGcc                                                                                  >  1:545631/1‑150 (MQ=255)
                                                              gATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGcc                                                                                  >  1:545591/1‑150 (MQ=255)
                                                              gATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGcc                                                                                  <  2:776647/150‑1 (MQ=255)
                                                                 ccGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTc                                                                               >  1:758971/1‑150 (MQ=255)
                                                                   gCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCaa                                                                             <  1:1139229/150‑1 (MQ=255)
                                                                    cTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAg                                                                            >  1:1886554/1‑150 (MQ=255)
                                                                        gACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCaaa                                                                        >  1:4137721/1‑150 (MQ=255)
                                                                        gACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCaaa                                                                        >  2:2894545/1‑150 (MQ=255)
                                                                          cTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAAcc                                                                      <  1:2187565/150‑1 (MQ=255)
                                                                              ggCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACt                                                                  >  1:3795046/1‑150 (MQ=255)
                                                                                       aTAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGa                                                         <  1:3670238/150‑1 (MQ=255)
                                                                                         aaaGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGAtt                                                       <  2:2536095/150‑1 (MQ=255)
                                                                                         aaaGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGAtt                                                       >  1:1308515/1‑150 (MQ=255)
                                                                                         aaaGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGAtt                                                       >  1:1307455/1‑150 (MQ=255)
                                                                                           aGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTgg                                                     <  1:1005911/150‑1 (MQ=255)
                                                                                               gCCGATCTGGATCCGTCTTTCCACGATCTGACCTAAGCAGAATTCCAGGAGAACTTCAACGTCGGTTCATTTGCCAGCGGCAAATAAACCATGAAACTCTGCGAGATGCTGGAAGCCCTCAAGCTAACCTACTGCGGCCCGATTGGTgac                                                 <  2:2088996/150‑3 (MQ=255)
                                                                                               gCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGcc                                                 <  2:214561/150‑1 (MQ=255)
                                                                                               gCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGcc                                                 <  2:1516482/150‑1 (MQ=255)
                                                                                               gCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGATAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGcc                                                 <  2:4339437/150‑1 (MQ=255)
                                                                                                ccGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCg                                                <  1:1733456/150‑1 (MQ=255)
                                                                                                 cGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGa                                               <  1:3944303/150‑1 (MQ=255)
                                                                                                  gATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAg                                              >  2:4196893/1‑150 (MQ=255)
                                                                                                   aTCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGt                                             <  1:566218/150‑1 (MQ=255)
                                                                                                       ggATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTatat                                         >  1:2331151/1‑150 (MQ=255)
                                                                                                         aTCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGc                                       >  2:2272715/1‑150 (MQ=255)
                                                                                                             gTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATAGCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACAt                                   <  2:1790340/150‑1 (MQ=255)
                                                                                                               cTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCTATTGGTGCCGAGTATATGCACATTa                                 <  1:3516986/150‑1 (MQ=255)
                                                                                                                tttCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTAc                                <  1:1485057/150‑1 (MQ=255)
                                                                                                                  tCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGTCCCGATTGGTGCCGAGTATATGCACATTACCa                              >  1:696894/1‑150 (MQ=255)
                                                                                                                  tCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCa                              >  2:3051579/1‑150 (MQ=255)
                                                                                                                  tCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCa                              >  2:3051619/1‑150 (MQ=255)
                                                                                                                  tCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCa                              >  1:4468174/1‑150 (MQ=255)
                                                                                                                      cGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCAc                          >  1:2959438/1‑150 (MQ=255)
                                                                                                                       gATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCAcc                         >  2:1527230/1‑150 (MQ=255)
                                                                                                                          cTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaa                      >  2:3869817/1‑150 (MQ=255)
                                                                                                                           tGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaag                     <  2:2900079/150‑1 (MQ=255)
                                                                                                                           tGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaag                     >  1:4200948/1‑150 (MQ=255)
                                                                                                                           tGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaag                     >  1:2635997/1‑150 (MQ=255)
                                                                                                                            gACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaaga                    >  1:288534/1‑150 (MQ=255)
                                                                                                                            gACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaaga                    <  2:4441390/150‑1 (MQ=255)
                                                                                                                            gACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATAGGTGCCGAGTATATGCACATTACCAGCACCgaaga                    >  1:4182886/1‑150 (MQ=255)
                                                                                                                             aCCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaagaa                   >  2:3240373/1‑150 (MQ=255)
                                                                                                                              ccGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaa                  >  2:3728476/1‑150 (MQ=255)
                                                                                                                              ccGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaa                  <  1:2996191/150‑1 (MQ=255)
                                                                                                                               cGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaaa                 <  2:1784627/150‑1 (MQ=255)
                                                                                                                                gAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaaaa                <  2:504451/150‑1 (MQ=255)
                                                                                                                                gAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaaaa                >  2:2645922/1‑150 (MQ=255)
                                                                                                                                  aGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACg              >  1:175703/1‑150 (MQ=255)
                                                                                                                                         ttCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATc       >  2:3641598/1‑150 (MQ=255)
                                                                                                                                             aGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAAc   <  1:2731341/150‑1 (MQ=255)
                                                                                                                                             aGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAAc   <  2:910739/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACA  >  CP009273/754427‑754718

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: