Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,373,985 C→A 88.0% intergenic (+37/‑117) yhcB → / → degQ DUF1043 family inner membrane‑anchored protein/serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,373,9850CA88.0% 706.7 / 99.7 284intergenic (+37/‑117)yhcB/degQDUF1043 family inner membrane‑anchored protein/serine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (10/24);  new base A (95/155);  total (105/179)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-02

Only 100 of 292 total aligned reads displayed.

AATCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGT  >  CP009273/3373838‑3374128
                                                                                                                                                   |                                                                                                                                               
aaTCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTtatc                                                                                                                                               >  2:3284168/1‑150 (MQ=255)
aaTCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTtatc                                                                                                                                               >  1:203000/1‑150 (MQ=255)
 aTCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTtctct                                                                                                                                              >  2:4353505/1‑150 (MQ=255)
 aTCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTtatct                                                                                                                                              >  2:3817398/1‑150 (MQ=255)
     tCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATctccct                                                                                                                                          <  1:3365900/150‑1 (MQ=255)
     tCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATctccct                                                                                                                                          <  2:122024/150‑1 (MQ=255)
         gCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTcccc                                                                                                                                      <  1:376828/150‑1 (MQ=255)
          cAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCg                                                                                                                                     >  2:2258959/1‑150 (MQ=255)
          cAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCg                                                                                                                                     >  1:4139288/1‑150 (MQ=255)
          cAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCACCTCTTATCTCCCTCCCCg                                                                                                                                     >  1:4139110/1‑150 (MQ=255)
           agagTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGa                                                                                                                                    >  1:4420316/1‑150 (MQ=255)
                cTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATc                                                                                                                               >  2:141529/1‑150 (MQ=255)
                 tGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCa                                                                                                                              <  1:3285827/150‑1 (MQ=255)
                 tGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCa                                                                                                                              <  2:1091931/150‑1 (MQ=255)
                     gCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTa                                                                                                                          <  2:2041592/150‑1 (MQ=255)
                       cAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAAt                                                                                                                        >  1:3475907/1‑150 (MQ=255)
                        aGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATc                                                                                                                       <  1:3929423/150‑1 (MQ=255)
                         gCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCt                                                                                                                      <  1:988146/150‑1 (MQ=255)
                         gCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCt                                                                                                                      <  1:3098256/150‑1 (MQ=255)
                         gCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCt                                                                                                                      <  1:2276185/150‑1 (MQ=255)
                          cAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTg                                                                                                                     <  1:268800/150‑1 (MQ=255)
                          cAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTg                                                                                                                     >  2:271216/1‑150 (MQ=255)
                             cGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGgtg                                                                                                                  <  1:3671060/150‑1 (MQ=255)
                                  aGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCa                                                                                                             >  1:329998/1‑150 (MQ=255)
                                    gCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCAtt                                                                                                           <  2:1038322/150‑1 (MQ=255)
                                      aCCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGt                                                                                                         <  2:2926329/150‑1 (MQ=255)
                                      aCCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGt                                                                                                         >  2:1708788/1‑150 (MQ=255)
                                      aCCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGt                                                                                                         <  2:1090390/150‑1 (MQ=255)
                                       ccGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGtt                                                                                                        >  1:906095/1‑150 (MQ=255)
                                       ccGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGtt                                                                                                        <  1:2918842/150‑1 (MQ=255)
                                        cGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTa                                                                                                       >  2:3744530/1‑150 (MQ=255)
                                          gTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGc                                                                                                     <  1:2192471/150‑1 (MQ=255)
                                          gTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGc                                                                                                     <  2:539146/150‑1 (MQ=255)
                                           tGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGcc                                                                                                    <  1:692971/150‑1 (MQ=255)
                                           tGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGcc                                                                                                    >  2:3644757/1‑150 (MQ=255)
                                            gCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCg                                                                                                   <  1:3149203/150‑1 (MQ=255)
                                                aTGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTg                                                                                               >  2:4327234/1‑150 (MQ=255)
                                                 tGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGa                                                                                              >  2:2227305/1‑150 (MQ=255)
                                                 tGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGa                                                                                              <  2:1487300/150‑1 (MQ=255)
                                                     tCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAAt                                                                                          <  2:4031512/150‑1 (MQ=255)
                                                       gcgACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTc                                                                                        <  1:1138663/150‑1 (MQ=255)
                                                         gACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCaa                                                                                      >  2:1942900/1‑150 (MQ=255)
                                                          aCTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAAt                                                                                     >  2:4112119/1‑150 (MQ=255)
                                                           cTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATa                                                                                    <  2:1619007/150‑1 (MQ=255)
                                                                cTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATc                                                                               <  2:800887/150‑1 (MQ=255)
                                                                cTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATc                                                                               <  2:803582/150‑1 (MQ=255)
                                                                cTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATc                                                                               <  1:1078375/150‑1 (MQ=255)
                                                                  gAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCaa                                                                             <  2:3434330/150‑1 (MQ=255)
                                                                     ggCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACt                                                                          <  2:929556/150‑1 (MQ=255)
                                                                     ggCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACt                                                                          <  1:3794462/150‑1 (MQ=255)
                                                                         cATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGttt                                                                      <  2:4404356/150‑1 (MQ=255)
                                                                          aTCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGtttt                                                                     >  2:2542970/1‑150 (MQ=255)
                                                                          aTCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGtttt                                                                     >  1:1379198/1‑150 (MQ=255)
                                                                                ccTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAtct                                                               >  2:2340425/1‑150 (MQ=255)
                                                                                ccTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAtct                                                               <  1:3236127/150‑1 (MQ=255)
                                                                                 ctgctgCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAtctc                                                              <  1:613239/150‑1 (MQ=255)
                                                                                    ctgcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCttt                                                           <  1:1999220/150‑1 (MQ=255)
                                                                                     tgcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCtttt                                                          <  2:1265859/150‑1 (MQ=255)
                                                                                      gcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTc                                                         <  1:4191501/150‑1 (MQ=255)
                                                                                       cGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCt                                                        >  1:2965790/1‑150 (MQ=255)
                                                                                         tACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAACCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTa                                                      <  1:1934917/150‑1 (MQ=255)
                                                                                            tGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCa                                                   >  2:1235505/1‑150 (MQ=255)
                                                                                             ggCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCAt                                                  <  2:3885013/150‑1 (MQ=255)
                                                                                               cgcgAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTc                                                <  2:3492670/150‑1 (MQ=255)
                                                                                               cgcgAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTc                                                <  2:1204491/150‑1 (MQ=255)
                                                                                                  gAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAgg                                             <  1:3026243/150‑1 (MQ=255)
                                                                                                   aaGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGt                                            <  1:4033589/150‑1 (MQ=255)
                                                                                                   aaGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGt                                            <  1:1980434/150‑1 (MQ=255)
                                                                                                    aGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTa                                           <  1:1796971/150‑1 (MQ=255)
                                                                                                    aGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTa                                           <  1:2353803/150‑1 (MQ=255)
                                                                                                         gACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACgaga                                      <  2:1463796/150‑1 (MQ=255)
                                                                                                           cTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGc                                    <  2:2692609/150‑1 (MQ=255)
                                                                                                             aaTTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAg                                  <  1:3504730/150‑1 (MQ=255)
                                                                                                              atttatttTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAgg                                 <  2:582575/150‑1 (MQ=255)
                                                                                                                    ttttCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGaataat                           <  2:3078154/150‑1 (MQ=255)
                                                                                                                    ttttCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGaataat                           <  1:1223220/150‑1 (MQ=255)
                                                                                                                     tttCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATg                          <  2:870646/150‑1 (MQ=255)
                                                                                                                     tttCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATg                          <  1:3974199/150‑1 (MQ=255)
                                                                                                                      ttCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGa                         <  1:451440/150‑1 (MQ=255)
                                                                                                                      ttCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGa                         <  1:1912302/150‑1 (MQ=255)
                                                                                                                      ttCGGGCGCAGACATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGa                         <  1:449733/150‑1 (MQ=255)
                                                                                                                       tCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGaa                        <  1:1004472/150‑1 (MQ=255)
                                                                                                                       tCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGaa                        >  2:1777592/1‑150 (MQ=255)
                                                                                                                         gggCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGaaaa                      >  1:4171219/1‑150 (MQ=255)
                                                                                                                          ggCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGaaaaa                     <  1:594553/150‑1 (MQ=255)
                                                                                                                            cgcAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAaaac                   <  2:4488391/150‑1 (MQ=255)
                                                                                                                             gcAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAaaaca                  >  1:165837/1‑150 (MQ=255)
                                                                                                                               aGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAaaacaaa                >  1:2545044/1‑150 (MQ=255)
                                                                                                                               aGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAaaacaaa                >  2:1507437/1‑150 (MQ=255)
                                                                                                                                 ccATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAAcc              >  1:686894/1‑150 (MQ=255)
                                                                                                                                   aTTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCa            >  1:3723554/1‑150 (MQ=255)
                                                                                                                                   aTTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCa            <  1:1280967/150‑1 (MQ=255)
                                                                                                                                      gcgcCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCt         <  1:2500028/150‑1 (MQ=255)
                                                                                                                                       cgcCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTg        <  1:3583337/150‑1 (MQ=255)
                                                                                                                                       cgcCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTg        <  2:3779263/150‑1 (MQ=255)
                                                                                                                                         cccTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGtt      >  1:4499356/1‑150 (MQ=255)
                                                                                                                                          cctcctCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTg     <  2:2965790/150‑1 (MQ=255)
                                                                                                                                           ctcctcTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGa    <  2:2373022/150‑1 (MQ=255)
                                                                                                                                            tcctcTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAg   <  2:4090009/150‑1 (MQ=255)
                                                                                                                                             cctcTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGt  >  1:790261/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AATCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGT  >  CP009273/3373838‑3374128

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: