Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,001,919 C→A 9.2% H351N (CAT→AAT)  ygeX → 2,3‑diaminopropionate ammonia lyase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,001,9190CA 9.2% 1143.1 / 68.9 326H351N (CAT→AAT) ygeX2,3‑diaminopropionate ammonia lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (133/163);  new base A (14/16);  total (147/179)
Fisher's exact test for biased strand distribution p-value = 8.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 330 total aligned reads displayed.

TCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCA  >  CP009273/3001770‑3002068
                                                                                                                                                     |                                                                                                                                                     
 ccTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCa                                                                                                                                                      >  2:555830/1‑150 (MQ=255)
 ccTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGAAGCGGTTCa                                                                                                                                                      <  1:3034855/150‑1 (MQ=255)
    aCGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTa                                                                                                                                                   <  1:4272768/150‑1 (MQ=255)
      gTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATc                                                                                                                                                 >  2:5142518/1‑150 (MQ=255)
      gTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATc                                                                                                                                                 <  1:92336/150‑1 (MQ=255)
       tAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCa                                                                                                                                                <  1:2981872/150‑1 (MQ=255)
       tAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCa                                                                                                                                                <  1:4092176/150‑1 (MQ=255)
        aaCTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCAc                                                                                                                                               <  1:2640124/150‑1 (MQ=255)
          cTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCAccc                                                                                                                                             <  1:3418714/150‑1 (MQ=255)
              gcCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCaa                                                                                                                                         <  1:3761140/150‑1 (MQ=255)
               ccACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAAc                                                                                                                                        <  1:458494/150‑1 (MQ=255)
               ccACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAAc                                                                                                                                        <  1:3092378/150‑1 (MQ=255)
                cACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACg                                                                                                                                       >  2:4185518/1‑150 (MQ=255)
                  cccAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTc                                                                                                                                     >  1:61666/1‑150 (MQ=255)
                  cccAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTc                                                                                                                                     <  2:316601/150‑1 (MQ=255)
                  cccAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTc                                                                                                                                     <  1:1995393/150‑1 (MQ=255)
                  cccAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTc                                                                                                                                     <  2:4376658/150‑1 (MQ=255)
                  cccAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTAGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTc                                                                                                                                     >  2:2664282/1‑150 (MQ=255)
                   ccAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCa                                                                                                                                    <  1:1937923/150‑1 (MQ=255)
                      aTTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAg                                                                                                                                 >  2:2279128/1‑150 (MQ=255)
                      aTTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAg                                                                                                                                 <  2:4648441/150‑1 (MQ=255)
                            ctcCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGAt                                                                                                                           >  1:4364372/1‑150 (MQ=255)
                            ctcCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGAt                                                                                                                           >  1:4491865/1‑150 (MQ=255)
                              ccTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGGCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATgg                                                                                                                         >  2:4069810/1‑150 (MQ=255)
                               cTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGa                                                                                                                        <  1:4892457/150‑1 (MQ=255)
                               cTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGa                                                                                                                        >  1:3683262/1‑150 (MQ=255)
                                  ccAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCTGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGaaaa                                                                                                                     <  1:217380/150‑1 (MQ=255)
                                    aGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAAc                                                                                                                   >  2:4983521/1‑150 (MQ=255)
                                        cAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGgc                                                                                                               <  1:4761320/150‑1 (MQ=255)
                                         aGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGgcg                                                                                                              >  2:2318366/1‑150 (MQ=255)
                                            gTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTg                                                                                                           >  2:4723339/1‑150 (MQ=255)
                                            gTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTg                                                                                                           <  2:1146735/150‑1 (MQ=255)
                                             ttGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGa                                                                                                          >  1:4922398/1‑150 (MQ=255)
                                              tGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGaa                                                                                                         <  2:3966749/150‑1 (MQ=255)
                                               gCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAAc                                                                                                        <  2:2908374/150‑1 (MQ=255)
                                                ccGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACa                                                                                                       <  1:3482669/150‑1 (MQ=255)
                                                ccGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACa                                                                                                       <  2:1778950/150‑1 (MQ=255)
                                                  gCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACaaa                                                                                                     <  1:4412942/150‑1 (MQ=255)
                                                   cATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAg                                                                                                    >  1:865835/1‑150 (MQ=255)
                                                    aTTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGa                                                                                                   >  2:2469048/1‑150 (MQ=255)
                                                       aGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGTCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGc                                                                                                >  1:1242061/1‑150 (MQ=255)
                                                        ggTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGcc                                                                                               >  2:1052073/1‑150 (MQ=255)
                                                         gTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCg                                                                                              <  1:2391612/150‑1 (MQ=255)
                                                          tATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAATGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCgt                                                                                             >  2:4516535/1‑150 (MQ=255)
                                                              cgcgTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCgtggtg                                                                                         <  1:483594/150‑1 (MQ=255)
                                                               gcgTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGc                                                                                        <  2:49086/150‑1 (MQ=255)
                                                                cgTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGACAATCTGGCGCTGAACAAAGATGCCGTGGTGCt                                                                                       >  2:3618065/1‑150 (MQ=255)
                                                                  tGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTgg                                                                                     <  1:389389/150‑1 (MQ=255)
                                                                        gTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCa                                                                               <  1:5021541/150‑1 (MQ=255)
                                                                           aTCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGACAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCa                                                                            >  2:3389182/1‑150 (MQ=255)
                                                                            tCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCAc                                                                           >  1:1709629/1‑150 (MQ=255)
                                                                                tACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTACATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGaa                                                                       >  2:1651427/1‑150 (MQ=255)
                                                                                 aCGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAg                                                                      <  1:1001535/150‑1 (MQ=255)
                                                                                 aCGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAg                                                                      <  1:632299/150‑1 (MQ=255)
                                                                                   ggCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGc                                                                    >  2:2877257/1‑150 (MQ=255)
                                                                                   ggCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGc                                                                    >  1:4552241/1‑150 (MQ=255)
                                                                                   ggCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGc                                                                    >  2:1868541/1‑150 (MQ=255)
                                                                                       aCGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGaca                                                                <  1:1693884/150‑1 (MQ=255)
                                                                                        cGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGacac                                                               <  2:827554/150‑1 (MQ=255)
                                                                                         gACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACAcc                                                              >  2:496903/1‑150 (MQ=255)
                                                                                          aCCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCg                                                             >  2:1619489/1‑150 (MQ=255)
                                                                                           cccGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGa                                                            >  1:1541649/1‑150 (MQ=255)
                                                                                            ccGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGag                                                           <  2:3571290/150‑2 (MQ=255)
                                                                                                  atcatcTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACTCCGACGTGAAg                                                     <  2:2589556/150‑1 (MQ=255)
                                                                                                     atcTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCAc                                                  <  1:1076802/150‑1 (MQ=255)
                                                                                                     atcTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCAc                                                  <  2:441197/150‑1 (MQ=255)
                                                                                                     atcTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCAc                                                  <  1:159158/150‑1 (MQ=255)
                                                                                                      tctcCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACt                                                 <  2:3756007/150‑1 (MQ=255)
                                                                                                        tcCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTAc                                               >  1:3245911/1‑150 (MQ=255)
                                                                                                        tcCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTAc                                               >  1:1440255/1‑150 (MQ=255)
                                                                                                         ccGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTAcc                                              >  1:52031/1‑150 (MQ=255)
                                                                                                         ccGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTAcc                                              >  1:2705407/1‑150 (MQ=255)
                                                                                                         ccGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTAcc                                              <  2:4204526/150‑1 (MQ=255)
                                                                                                            gTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACcgcg                                           >  1:123592/1‑150 (MQ=255)
                                                                                                                 tCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGtt                                      >  2:2088032/1‑150 (MQ=255)
                                                                                                                      cgcgGTCGGTTTGGGCGGGCTCGCAGCGGTACATTATCACCCGCAACGTCACAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTg                                 <  2:3554039/150‑1 (MQ=255)
                                                                                                                      cgcTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTg                                 <  1:24949/150‑1 (MQ=255)
                                                                                                                       gcTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTgg                                <  2:2511721/150‑1 (MQ=255)
                                                                                                                       gcTGTCGGTTTGGGCGGTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGCAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTgg                                <  2:223053/150‑1 (MQ=255)
                                                                                                                        cTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTggg                               >  2:1409454/1‑150 (MQ=255)
                                                                                                                         tGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACCGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGa                              <  1:2520600/150‑1 (MQ=255)
                                                                                                                           tCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGACCACAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAg                            <  1:2851058/150‑1 (MQ=255)
                                                                                                                             ggTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGc                          >  1:3995717/1‑150 (MQ=255)
                                                                                                                                 tGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAac                      <  2:2793009/150‑1 (MQ=255)
                                                                                                                                   ggCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAacac                    <  2:3139908/150‑1 (MQ=255)
                                                                                                                                   ggCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAacac                    >  2:4842404/1‑150 (MQ=255)
                                                                                                                                   ggCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAacac                    <  2:2995866/150‑1 (MQ=255)
                                                                                                                                   ggCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAacac                    <  2:507742/150‑1 (MQ=255)
                                                                                                                                     cGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGc                  >  1:2072835/1‑150 (MQ=255)
                                                                                                                                      gTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCt                 >  2:3331443/1‑150 (MQ=255)
                                                                                                                                       ttCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTg                >  2:391355/1‑150 (MQ=255)
                                                                                                                                         ctcGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTa              >  1:3860038/1‑150 (MQ=255)
                                                                                                                                         ctcGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTa              <  1:2279128/150‑1 (MQ=255)
                                                                                                                                           cGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGc            <  2:212003/150‑1 (MQ=255)
                                                                                                                                           cGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGc            <  2:696736/150‑1 (MQ=255)
                                                                                                                                              aGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCAcc         >  1:817359/1‑150 (MQ=255)
                                                                                                                                                 ggTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTa      >  2:5179992/1‑150 (MQ=255)
                                                                                                                                                 ggTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTa      >  2:45903/1‑150 (MQ=255)
                                                                                                                                                  gTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTaa     <  1:4628165/150‑1 (MQ=255)
                                                                                                                                                    tCAGTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATc   <  1:3526120/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCA  >  CP009273/3001770‑3002068

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: