Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,551,620 C→A 70.0% D226Y (GAC→TAC)  ddpF ← D,D‑dipeptide permease system, ATP‑binding component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,551,6200CA70.0% 327.3 / 303.6 310D226Y (GAC→TAC) ddpFD,D‑dipeptide permease system, ATP‑binding component
Reads supporting (aligned to +/- strand):  ref base C (52/41);  new base A (103/113);  total (156/154)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 314 total aligned reads displayed.

AAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGATCTGCGCCTGAACG  >  CP009273/1551471‑1551768
                                                                                                                                                     |                                                                                                                                                    
aaaaCAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTc                                                                                                                                                      <  1:4181532/150‑1 (MQ=255)
  aaCAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGc                                                                                                                                                    >  1:3871223/1‑150 (MQ=255)
    cAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGccc                                                                                                                                                  >  1:4525018/1‑150 (MQ=255)
     aGCATTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTACTGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGATTCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCa                                                                                                                                                 <  1:1938949/150‑1 (MQ=255)
      gCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAg                                                                                                                                                <  1:3357818/150‑1 (MQ=255)
      gCCTTGCGGCAACGTGCGGTTTCCCGGCAGAGCCGTTTTACGTAATGCCCCTTCTTCCTCCAGCGGTTTATCACTGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAg                                                                                                                                                <  2:2081395/150‑1 (MQ=255)
        cTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGtt                                                                                                                                              <  2:4735402/150‑1 (MQ=255)
         ttGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTc                                                                                                                                             <  2:3355305/150‑1 (MQ=255)
          tGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCt                                                                                                                                            >  1:3422318/1‑150 (MQ=255)
          tGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCt                                                                                                                                            <  2:1596797/150‑1 (MQ=255)
           gCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTa                                                                                                                                           <  1:3631682/150‑1 (MQ=255)
           gCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTa                                                                                                                                           <  1:480671/150‑1 (MQ=255)
               cAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAAt                                                                                                                                       <  2:1396643/150‑1 (MQ=255)
               cAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAAt                                                                                                                                       <  1:270274/150‑1 (MQ=255)
                aaCGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATc                                                                                                                                      >  2:1428210/1‑150 (MQ=255)
                  cGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTg                                                                                                                                    >  2:3277782/1‑150 (MQ=255)
                  cGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTg                                                                                                                                    <  1:2832518/150‑1 (MQ=255)
                   gTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGc                                                                                                                                   >  1:358575/1‑150 (MQ=255)
                          ttCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGAt                                                                                                                            <  2:4598687/150‑1 (MQ=255)
                           tCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATa                                                                                                                           >  2:1277593/1‑150 (MQ=255)
                            cccGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATAc                                                                                                                          <  2:2024428/150‑1 (MQ=255)
                              cGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACAt                                                                                                                        >  2:3983492/1‑150 (MQ=255)
                               ggCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATc                                                                                                                       <  1:5028339/150‑1 (MQ=255)
                                  aGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAAAAATAGGGTGTATGGATGTGCAGGTGCCGTCAGCACCAGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCAcc                                                                                                                    >  1:1950406/1‑150 (MQ=255)
                                   gATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCg                                                                                                                   >  2:2807889/1‑150 (MQ=255)
                                    aTCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGc                                                                                                                  >  2:82854/1‑150 (MQ=255)
                                      ccGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCa                                                                                                                <  2:2317829/150‑1 (MQ=255)
                                        gTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTGGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCAcc                                                                                                              >  1:3017844/1‑150 (MQ=255)
                                            tACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGAt                                                                                                          <  2:4189979/150‑1 (MQ=255)
                                             aCGTAATGCCCATTCTTCCTCAAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATc                                                                                                         <  1:4163287/150‑1 (MQ=255)
                                              cGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCg                                                                                                        <  1:2637212/150‑1 (MQ=255)
                                               gTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGc                                                                                                       >  1:5302529/1‑150 (MQ=255)
                                                  aTGCCCATTCTTCCTCCAGCGGTTTATCAAAGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGAGGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGATCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCa                                                                                                    >  2:963018/1‑150 (MQ=255)
                                                     cccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCatat                                                                                                 <  2:3456373/150‑1 (MQ=255)
                                                     cccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCatat                                                                                                 >  1:811408/1‑150 (MQ=255)
                                                      ccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATg                                                                                                >  1:1505129/1‑150 (MQ=255)
                                                       cATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGa                                                                                               >  1:2198572/1‑150 (MQ=255)
                                                       cATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGa                                                                                               >  1:3492867/1‑150 (MQ=255)
                                                       cATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGa                                                                                               >  2:2460740/1‑150 (MQ=255)
                                                        aTTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGAc                                                                                              >  2:3229278/1‑150 (MQ=255)
                                                          tcttcCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGt                                                                                            <  2:2198359/150‑1 (MQ=255)
                                                           cttcCTCCAGCGGTTTATCAATGGCGGTGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGACCGAGATACATCACCGCCACCCGATCGCTCATATTACGTa                                                                                           >  1:3145229/1‑150 (MQ=255)
                                                              cctccAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCa                                                                                        >  2:2371564/1‑150 (MQ=255)
                                                                 ccAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCg                                                                                     >  2:5299095/1‑150 (MQ=255)
                                                                  cAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCga                                                                                    >  1:2104977/1‑150 (MQ=255)
                                                                  cAGCGGTCTATCCCTGGCGGGGCGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCGACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCga                                                                                    <  2:3736357/150‑1 (MQ=255)
                                                                   aGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCgag                                                                                   <  2:4493529/150‑1 (MQ=255)
                                                                         ttATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtg                                                                             <  1:3329804/150‑1 (MQ=255)
                                                                          tATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtgt                                                                            >  1:4485523/1‑150 (MQ=255)
                                                                           aTCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtgtg                                                                           <  1:2152739/150‑1 (MQ=255)
                                                                            tCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtgtgt                                                                          <  1:1230695/150‑1 (MQ=255)
                                                                             cAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtgtgtg                                                                         >  1:4759121/1‑150 (MQ=255)
                                                                              aaTGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGa                                                                        <  1:1322997/150‑1 (MQ=255)
                                                                               aTGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGaa                                                                       >  2:485202/1‑150 (MQ=255)
                                                                                tGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGaaa                                                                      <  2:567713/150‑1 (MQ=255)
                                                                                 ggCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAAt                                                                     <  2:1563740/150‑1 (MQ=255)
                                                                                 ggCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAAt                                                                     >  2:3497789/1‑150 (MQ=255)
                                                                                  gCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATc                                                                    <  1:4243506/150‑1 (MQ=255)
                                                                                   cGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCa                                                                   <  2:1377490/150‑1 (MQ=255)
                                                                                   cGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCa                                                                   >  2:4432028/1‑150 (MQ=255)
                                                                                        agggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGcaca                                                              >  1:980400/1‑150 (MQ=255)
                                                                                         gggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACAt                                                             <  2:2091805/150‑1 (MQ=255)
                                                                                           gaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAg                                                           <  2:3141704/150‑1 (MQ=255)
                                                                                            aaTCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAgg                                                          >  2:5001802/1‑150 (MQ=255)
                                                                                               ccAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCa                                                       >  1:2289497/1‑150 (MQ=255)
                                                                                                cAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAg                                                      <  1:3983492/150‑1 (MQ=255)
                                                                                                    aataatCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCg                                                  >  1:1305556/1‑150 (MQ=255)
                                                                                                    aataatCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCg                                                  >  1:1534788/1‑150 (MQ=255)
                                                                                                     ataatCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGt                                                 <  1:2574758/150‑1 (MQ=255)
                                                                                                      taatCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTg                                                >  1:105556/1‑150 (MQ=255)
                                                                                                        atCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGAt                                              >  2:430882/1‑150 (MQ=255)
                                                                                                         tCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGAtt                                             <  2:5047997/150‑1 (MQ=255)
                                                                                                              tgtATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCt                                        >  1:1781670/1‑150 (MQ=255)
                                                                                                                tATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGt                                      <  2:3430504/150‑1 (MQ=255)
                                                                                                                tATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGACCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGt                                      <  2:2646649/150‑1 (MQ=255)
                                                                                                                 aTGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTa                                     >  2:3845425/1‑150 (MQ=255)
                                                                                                                  tGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAg                                    <  2:2584596/150‑1 (MQ=255)
                                                                                                                  tGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAg                                    <  1:632044/150‑1 (MQ=255)
                                                                                                                   ggATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGc                                   <  1:4411978/150‑1 (MQ=255)
                                                                                                                      tgtgCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGtt                                >  1:2427390/1‑150 (MQ=255)
                                                                                                                         gCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTAcc                             <  2:1938017/150‑1 (MQ=255)
                                                                                                                             gTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTa                         >  2:2194107/1‑150 (MQ=255)
                                                                                                                             gTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTa                         >  2:1026120/1‑150 (MQ=255)
                                                                                                                              tGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTaa                        >  1:1027331/1‑150 (MQ=255)
                                                                                                                                 cGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAAtt                     <  1:4648512/150‑1 (MQ=255)
                                                                                                                                 cGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAAtt                     <  1:929347/150‑1 (MQ=255)
                                                                                                                                     aGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAgg                 >  1:2945726/1‑150 (MQ=255)
                                                                                                                                     aGCACCTGCTGCGCGACGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAgg                 >  1:1655519/1‑150 (MQ=255)
                                                                                                                                      gCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGa                <  1:684820/150‑1 (MQ=255)
                                                                                                                                       cACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGAt               <  2:260898/150‑1 (MQ=255)
                                                                                                                                       cACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGAt               <  1:4068033/150‑1 (MQ=255)
                                                                                                                                        aCCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGATc              <  2:4106162/150‑1 (MQ=255)
                                                                                                                                        aCCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGATc              <  2:4131090/150‑1 (MQ=255)
                                                                                                                                        aCCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGATc              >  2:3089982/1‑150 (MQ=255)
                                                                                                                                         ccTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGATCt             <  2:172292/150‑1 (MQ=255)
                                                                                                                                             ctgcGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGATCTgcgc         >  2:894767/1‑150 (MQ=255)
                                                                                                                                              tgcGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGATCTGCGcc        >  1:4957972/1‑150 (MQ=255)
                                                                                                                                              tgcGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGATCTGCGcc        <  2:319694/150‑1 (MQ=255)
                                                                                                                                               gcgcGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGATCTGCGCCt       <  2:3610355/150‑1 (MQ=255)
                                                                                                                                                    tCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGATCTGCGCCTGAACg  <  2:4386669/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGATCTGCGCCTGAACG  >  CP009273/1551471‑1551768

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: